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Resumen 



En esta tesis, presentamos un estudio de la variacion de la masa del 
electron mg, y de la constante de estructura fina a, en distintas etapas de la 
evolucion del Universe. 

Estudiamos la epoca de formacion del hidrogeno neutro. Analizamos el 
escenario de recombinacion incluyendo los detalles en la recombinacion del 
helio, para encontrar las dependencias de las cantidades fisicas con las cons- 
tantes fundamentales. Modificamos las ccuacioncs para las fraccioncs de io- 
nizacion, para poder considerar valores arbitrarios de las constantcs. 

Utilizando los datos actualizados del fondo cosmico de radiacion, y el 
espectro de potencias de la distribucion de galaxias del catalogo 2dFGRS, 
buscamos Kmites a la posible variacion de las constantes en la epoca de re- 
combinacion. Estimamos el valor dc un conjunto de parametros cosmologicos, 
junto con la variacion de las constantes. Analizamos la variacion de a y rUe 
por separado, y el caso en que ambas varian simultaneamente. Presentamos 
resultados actualizados y comparamos con las cotas existentes. Estimamos 
tambien una relacion fenomenologica entre la variacion de a y la variacion 
de rrie entre el desacople y la actualidad. 

Analizamos el modelo fenomenologico de Barrow-Magueijo, que propone 
una variacion de la masa del electron inducida por cambios de un campo 
escalar en el espacio-tiempo. Mejoramos las soluciones originalmente presen- 
tadas por dichos autores. Estimamos los parametros del modelo usando los 
limites para la variacion de la masa del electron que se obtienen en esta tesis 
a partir de los datos del fondo cosmico de radiacion, y otros datos astronomi- 
cos y de laboratorio. Por otra parte, utilizando resultados de experimentos 
que testean la validez del Principio de Equivalencia Debil, encontramos una 
cota para el parametro del modelo, que viene de la variacion en la masa 
del electron inducida por cambios espaciales del campo escalar. Comparando 
ambos limites, se llega a la conclusion de que el modelo debe ser descartado. 

Por ultimo, rcanalizamos los cfcctos dc la variacion dc a o m^, en los radios 
planctarios. Explicamos por que la variacion dc la masa del electron no induce 
cambios en el radio factibles de ser observados. Presentamos una cota local 
e independiente para la variacion de a, obtenida a partir de principios muy 
generales. 
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Capitulo 1 
Introduccion 



La descripcion actual del Universo tiene su fundamento en dos teorias 
bien establecidas: la relatividad general, y el modelo estandar de fi'sica de 
particulas. En ambas existen parametros, Uamados "constantes fundamenta- 
les", cuyos valores no quedan determinados a partir de las teorias mismas, 
sino que deben ser hallados experimentalmente. Tal vez, uno de los objetivos 
mas ambiciosos de la fisica moderna sea formular una teoria del "todo" , una 
teon'a que unifique todas las interacciones de la naturaleza, donde los valores 
de diclios parametros puedan ser deducidos a partir de primeros principios. 
Hasta hoy, se han podido unificar las interacciones electromagneticas con la 
fuerza debil en las interacciones electro-debiles, y estas con las interacciones 
fuertes, en lo que se conoce como Teorias de Gran Unificacion. Sin embargo, 
aiin no se ha logrado encontrar una teon'a cuantica de la gravedad. 

Las teorias que intentan unificar la gravedad con las demas interacciones 

son de yariada naturaleza . Entre ellas se encuentran las teorias de super- 

cuerdas flWu fc W ang"l986': i MaedalllQSsI : iBarr fc MohaDatralllQsdjPamour fc Polvakov 



1993; Id amour et al.i.2002a .lj') , teorias de mundos br ana (|Younill2001a|jbl:|Palma et al] 



2003; Brax et al 



2003) y teorias de Kaluza-Klein (Kaluza 



Weinberglll983l : iGleiser fc Taylorlll985l : lOverduin fc Wesson 



19211: lKleinlll926: 



19971 ) . Las energfas 



involucradas en los fenomenos que permiten corroborar la unificacion son tan 
altas, que las posibles predicciones de las teorias serian muy diffciles de ve- 
rificar. Sin embargo, algunas consecuencias a bajas energfas son factibles de 
ser comprobadas a traves de experimentos u observaciones. Una de las pre- 
dicciones mas importantes de estas teorias es que el valor de las constantes 
fundamentals puede variar en distintos puntos del espacio-tiempo. La de- 
pendencia de las constantes con el tiempo es diferente para distintas teorias. 
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Es por ello que a traves de la comparacion con resultados de experimentos y 
observaciones es posible poner cotas a las mismas. 

Por otra parte, existen teon'as efectivas a bajas energias, formuladas para 
describir la variacion de determinadas constantes fundamentales. En este ca- 
so, se asocia el valor de las constantes a un campo en el espacio-tiempo, cuyo 
Lagrangiano se propone de manera tal de cumplir principios muy generales 
y las ecuaciones de movimiento tienen soluciones que pueden ser compara- 
das con datos observacionales. A modo de ejemplo, mencionamos la teoria 
de Bekenstein pa ra la variacion temporal de la constante de estructura fina 
(jBekensteirull982l ) y el modelo de Barrow fc M agueijo para la variacion de la 
masa del electron (IBarrow fc Magueijdl2005l ). 

Se ha estudiado tambien la variacion de otras constantes fundamentales. 



Albrecht fc Magueijo! (Il999l ) introdujeron un modelo cosmologico en el cual la 



velocidad de la luz es variable, con el fin de resolver los problemas basicos de 
la cosmologia (horizonte, planaridad, constante cosmologica, entropi'a y ho- 
mogeneidad) sin ac udir al paradigma inflacionario . Este r aodelo fue extendido 
por y arios autores fIClavton fc Moffatll200ll . l2000l . Il999l : iKehadas fc Kiritsis 
1999h . Asimismo, existen teorias que predicen la variacion de la constante de 



gravitacion universal G (IBrax fc David 1200 ll ). 

Una de las constantes fundamentales cuya variacion ha despertado ma- 
yor interes es la constante de estructura fina a, que esta relacionada con 
la constante de acoplamiento de las interacciones electromagneticas. Esta es 
una de las constantes estudiadas en esta tesis. Tambien ha suscitado interes 
la posible variacion de la constante de las interacciones debiles, parametri- 
zada por la constante de acopla miento de Fermi Gp. Sin embargo, como se 
enfatiza en iDixit fc Sherl (jl988l ). la constante de Fermi no es una constante 
fundamental en el mismo sentido que lo es la constante de estructura fina. 
Dado que Gp oc es mas apropiado examinar la variacion del valor de 

expectacion de vacio del campo Higgs (0). En particular, nos hemos concen- 
trado en la variacion de la masa en reposo del electron me, ya que es la linica 
consecuencia que tiene la variacion del Higgs en las etapas del Universo en 
las que esta tesis se focaliza. 

Desde el punto de vista observacional, el tema de la variacion de las cons- 
tantes es de gran interes. Son muchos los esfuerzos para poner cotas a la 
variacion de constantes como a y rrie, o el cociente /i = me/mp, entre la 
masa del electron y la masa del proton. Existen metodos locales, astronomi- 
cos y cosmologicos. Cada metodo considera el comportamiento temporal de 
las constantes a distintas escalas de energia y de distancia, y con distinta 
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precision, resultando todos de gran utilidad. Localmente, se pon en cotas a 



la tasa de variacion de l as constantes usando relojes atomicos (IBize et al 



2003; 



Fischer et al 



200G; 



Marion et al 



20041: iPeik et al.ll2004l : iPrestage et al.lll995l : ISortais et al 



20031 ). La variacion temporal de constantes a corrimien- 
tos al rojo (en adelante, redshifts z) del orden de z ~ 0,1 4 se puede acotar 



con datos del reactor n uclear natural que opero en Oklo fiFuiii et aljlioool : 



Damour fc Dysonlll996l ). Otras cotas se ob tienen a partir de la vida media de 



emisores de partfculas /3 (lOlive et al.ll2004l ). For otra parte, datos astronomi- 



2001 



variacion de a y u 


Webb et al. 


1999, 


Ivanchik et al. 


2005 


Tzanavaris et al. 



20071). Sin embargo, otros estudios de datos similares, arrojan resultados nu- 



los (iMartmez Fiorenzano et al.l l2003l: iQuast et al.ll2004j : iBahcall et all 12004 : 



Srianand et al.ll2004f ) 



Para estudiar el Universo temprano, y contrastar modelos alternativos a 
la cosmologia estandar, se pueden utilizar los datos de las abundancias pri- 
mitivas de elementos quimicos livianos, producidos durante la nucleosmtesis 
primordial (BBN), cuando el Universo tenia ~ 3 minutos de edad; y los datos 
del fondo cosmico de radiacion (CMB), que proveen informacion de la epo- 
ca de formacion del liidrogeno neutro, 400.000 afios despues del Big Bang. 
Dado que los procesos fisicos relevantes en ambas epocas son sensibles al 
valor de las constantes fundamentales, estos observables sirven como herra- 
mientas para poner limites a su variacion temporal. A pesar de que las cotas 
obtenidas a partir de BBN y de CMB son mas debiles que las mencionadas 
anteriormente, resultan sumamente importantes debido a que se refieren a 
las primeras etapas del Universo. 

En esta tesis se realizo un estudio a nivel cosmologico y local de la varia- 
cion temporal de dos constantes: la constante de estructura fina y la masa 
del electron. Uno de los objetivos mas importantes fue el estudio de la varia- 
cion de estas constantes durante la epoca en la cual el plasma de electrones 
e iones del Universo se recombino. Esta epoca se conoce con el nombre de 
recombinacion del Universo. 

El estudio de las dependencias con las constantes fundamentales de las 
cantidades fisicas relevantes durante la formacion del hidrogeno neutro se de- 
talla en el capitulo[2l AUi, ademas, se analiza por primera vez la infiuencia las 
constantes fundamentales en el escenario de recombinacion actualizado, que 
tiene en cuenta la recombinacion detallada del helio, el efecto del liidrogeno 
neutro en la recombinacion de este y la contribucion de las transiciones desde 
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estados tripletes del hidrogeno para la recombinacion. 

Utilizando el conjunto mas completo y actualizado de datos del fondo 
cosmico de radiacion, junto con datos del catalogo de galaxias 2dFGRS para 
el espectro de potencias de la distribucion de la materia, se ban obtenido 
cotas para la variacion de a y me en la epoca de formacion del hidrogeno 
neutro. Se considero la variacion de dichas constantes de manera separa- 
da, y tambien se realizo un analisis de la variacion conjunta de ambas, de 
manera totalmente fenomenologica, sin suponer ninguna relacion entre las 
variaciones, de manera tal que las conclusiones obtenidas en esta tesis resul- 
tan independientes de cualquier modelo. El linico trabajo an terior donde se 



consid ero la variacion conjunta de dichas constantes es el de llchikawa et al. 



(120061 ) ■ en el cual las variaciones estan relacionadas en el contexto de modelos 
dilatonicos de teorfas de cuerdas. En cambio, el resultado mostrado en esta 
tesis se obtiene suponiendo que las constantes varian en forma independiente 
una de la otra. A partir del ajuste, obtenemos una relacion fenomenologica 
entre las variacione s de a y mp en la epo ca de recombinacion. Comparamos 
con el resultado de llchikawa et al.l (120061 ) . y conclufmos que la propuesta de 
estos autores deben'a ser descartada. 

Los datos mas recientes del fondo cosmico de radiacion que habfan sido 
utilizados en esta tesis fueron los recolectados en los tres primeros afios de 
funcionamiento del satelite WMAP (WMAP3). En marzo de 2008, el satelite 
libero datos recolectados durante los primeros 5 anos (WMAP5). De mane- 
ra simultanea, el escenario de recombinacion fue actualizado, considerando 
procesos fisicos que habfan sido despreciados en el calculo estandar, y que 
afectan de manera sutil, aunque no despreciable, la historia de recombina- 
cion del Universo. For esta razon, se repitio el analisis estadistico para la 
variacion de a y m-e, incluyendo los datos de WMAP5, y el nuevo escenario 
de recombinacion, y se comparo con los resultados obtenidos con los datos 
de WMAP3, y el escenario de recombinacion estandar. 

En el Capitulo [3] se muestran los distintos ajustes realizados. Dado que 
se han llevado a cabo analisis estadisticos variando no solo las constantes 
fundamentals, sino tambien el conjunto basico de parametros cosmologicos, 
en esta tesis se ha abarcado el estudio y adaptacion de un programa numeri- 
co que utiliza cadenas de Markov Monte Carlo para muestrear el espacio 
de parametros, y hacer un analisis estadistico desde un enfoque bayesiano. 
Este metodo permite la exploracion del espacio de parametros en un tiempo 
computacionalmente mas accesible que el metodo frecuentista del calculo de 
la Likelihood en cada punto de una grilla en el espacio de parametros. En 
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dicho capi'tulo tambien se muestran y comparan los resultados obtenidos. 

Las cotas observacionales para la variacion de las constantes fundamenta- 
les son una herramienta importante para evaluar teon'as, por ejemplo, aque- 
llas que intentan unificar las cuatro interacciones de la naturaleza. En el 
Capitulo m se analiza el modelo de Barrow & Magueijo, que predice la varia- 
cion de la masa del electron como funcion del tiempo, a partir de la existencia 
de un campo escalar. Se obtuvieron mejores soluciones aproximadas del mo- 
delo, lo cual se logro teniendo en cuenta de manera mas rigurosa la evolucion 
de la expansion del Universo, y sin despreciar terminos en las ecuaciones di- 
ferenciales que describen la evolucion del campo escalar. Luego, se utilizaron 
las cotas para la variacion de rrie obtenidas en esta tesis a partir de los da- 
tos del fondo cosmico de radiacion, junto con otras cotas dadas a diferentes 
tiempos cosmologicos, para poner Ifmites a los parametros del modelo. Adi- 
cionalmente, se derivaron cotas a partir de resultados de experimentos que 
testean la validez del Principio de Equivalencia Debil, y se compararon con 
las anteriores. A la luz de los resultados, se Uego a la conclusion de que el 
modelo de Barrow & Magueijo debe ser descartado. 

Con el objetivo de estudiar la variacion de las constantes fundamentales 
a nivel local, tanto en el espacio como en el tiempo, se realizo un estudio 
en escalas correspondientes al sistema solar. Se analizo las consecuencias 
que un valor diferente de a y/o mg acarrean para el radio de la Tierra, la 
Luna y el planeta Mercurio. A partir de consideraciones generales, se deduce 
una expresion para la variacion de este parametro planetario. Utilizando los 
modelos mas actuales de la estructura interna de diclios cuerpos celestes y 
limites existentes para la variacion del radio de estos planetas, se dan cotas 
a la variacion de las constantes fundamentales entre la actualidad y la epoca 
de formacion del sistema solar. Asimismo, se analiza cuales de las constantes 
fundamentales estudiadas en esta tesis resultan relevantes en este problema. 
Los resultados se presentan en el Capitulo O 

Finalmente, en el Capitulo [6], se resumen los resultados de esta tesis, y se 
destacan las conclusiones mas relevantes. 

Parte de los resultados de esta tesis fueron incluidos en las siguientes 
publicaciones en revistas internacionales con referato: 

■ "Time variation of the fine structure constant and the Bekenstein mo- 
del". Mosquera, M. E., Scoccola, C. C, Landau, S. J., and Vucetich, 
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H. Astronomy and Astrophysics 478, 675 (2008). 

■ "Time variation of the electron mass and the B arrow- Magueijo model" . 
Scoccola, C. G., Mosquera, M. E., Landau, S. J., and Vucetich, H. 
Astrophysical Journal 681, 737 (2008). 

■ "Early Universe constraints on time variation of fundamental cons- 
tants". Landau, S. J., Mosquera. M. E., Scoccola, C. G., and Vucetich, 
H. Physical Review D 78 (8), 083527 (2008). 

■ "WMAP 5-year constraints on time variation of a and mg in a detailed 
recombination scenario". Scoccola, C. G., Landau, S. J., and Vucetich, 
H. Physics Letters B 669, 212 (2008). 
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Capitulo 2 

Fisica de Recombinacion 



En el modelo cosmologico estandar, la epoca de recombinacion es aque- 
Ua cuando el Universo se vuelve lo suficientemente frio como para que los 
protones, y niicleos de elementos qmmicos livianos, sean capaces de rete- 
ner electrones y se formen los atomos neutros. Los electrones, capturados en 
distintos niveles de energi'a, decaen al nivel fundamental. Sin embargo, este 
proceso se ve impedido por las ionizaciones y excitaciones que se producen 
rapidamente gracias al enorme reservorio de fotones de bajas energfas que 
hay en el Universo. Ademas, los decaimientos del nivel n = 2 al nivel n = 1 
y las recombinaciones directas desde el continuo al nivel fundamental se ven 
fuertemente disminufdos debido a que en ambos casos el foton emitido tiene 
una energfa en el rango donde existen muy pocos fotones de fondo, y por lo 
tanto sera inmediatamente absorbido por un atomo cercano que este en el 
nivel fundamental. Estos fenomenos traen como consecuencia que la recom- 
binacion no sea instantanea. Eventualmente, los atomos alcanzan el estado 
fundamental por el corrimiento al rojo de los fotones Lya (2p —>■ Is) causado 
por la expansion, o por la transicion 2s — *■ Is via dos fotones. Sin embargo, el 
Universo se expande y enfria antes de que la recombinacion termine, y queda 
una pequeiia fraccion de electrones y nucleos libres remanente. 

El elemento mas abundante en el Universo es el hidrogeno, y en menor 
proporcion el helio, siendo las abundancias de los demas elementos qufmicos 
formados en la nucleosfntesis primordial despreciables para el proceso de la 
recombinacion . L os primeros estudios e n este campo fueron realizados por 



Peebles! ( 119681 ) y IZeldovich et al.l (Il968l ). La metodologfa estandar era con- 
siderar un atomo de tres niveles (fundamental, n = 2 j continuo) con los 
niveles con n > 2 representados por un coeficiente de recombinacion. Se de- 
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rivaba una ecuacion diferencial ordinaria para la fraccion de ionizacion, pero 
debian hacerse muchas suposiciones para llegar a esta: i) el unico elemento 
en el Universe es el hidrogeno, ii) sus estados excitados estan en equilibrio 
con la radiacion, iii) se puede despreciar la emision estimulada para la transi- 
cion 2p —>■ Is (Imea Lya), iv) se puede usar un coeficiente de recombinacion 
simple, v) cada recombinacion neta resulta en un atomo en el nivel funda- 
mental, vi) el corrimiento al rojo de los fotones Lja se puede tener en cuenta 
usando una probabilidad de escape simple, y vii) los procesos colisionales 
son despreciables. Este tratamiento de la recombinacion fue usado duran- 
te treinta aiios con pocas modificaciones. Motivados por la expectativa de 
medir con gran precision las anisotropi'as de la radiacion cosmica de fondo 
con satelites y otro s experimentos dise fiados a tal fin, y gracias al poder de 
computo moderno, ISeager et al.l (I2OOOI ) presentaron un calculo detallado de 
la recombinacion del hidrogeno y del helio. El tratamiento consiste en consi- 
derar atomos con multiples niveles, y tener en cuenta todas las transiciones 
entre ellos, y tambien con el continuo, y tratar la recombinacion del hidrogeno 
y del helio en forma simultanea, teniendo en cuenta la evolucion de la po- 
blacion de el e ctrone s libres y de la temperatura de la materia. Ademas, en 
Seager et al.l (Il999al ). estos autores presentaron un conjunto de ecuaciones 
que permiten reproducir los resultados de los calculos detallados via para- 
metrizaciones en los coeficientes de recombinacion. En esta tesis Uamaremos 
escenario de recombinacion estandar a este tratamiento de la misma. 



Mas recientemente, la recombinacion del helio f ue tratada con mas preci- 
sion. Switzer y Hirata en una serie de tres trabaj os (ISwitzer fc Hiratall2008a : 
Hirata &: Switzerll2008l : ISwitzer fc Hiratall2008bl ). discutieron varios aspectos 
de la recombinacion y estudiaron los efectos de la retro-alimentacion produ- 
cida por las distorsiones espectrales, las lineas semi-prohibidas y prohibidas, 
el transporte radiativo en las lineas del He I con redistribucion parcial, y 
la opacidad del continuo producida por la fotoionizacion del hidrogeno. El 
efecto del hidrogeno neutro sobre la recombinacion del helio es acelerar la 
recombinacion de este respecto de los resultados obtenidos en Seager et al, 
debido a que fotones resonantes del He I ionizan atomos neutros de hidrogeno 
antes de ionizar otro atomo de helio, resultando en una recombinacion neta de 



helio. Este ultimo efecto fue tratado tambien por iKholupenko et al.l (120071 ) . 
quienes propusieron ademas una formula aproximada litil para ser inclufda 
en codigos como Recfast, que resuelve numericamente las ecuaciones de 
recombinacion. Nos referiremos a este tratamiento como escenario detallado. 
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2.1. Ecuaciones de recombinacion 



El conjunto de ecuaciones de recombinacion, conformado por dos ecua- 
ciones diferenciales ordinarias para la fraccion de ionizacion de H y de He I y 
una ecuacion para la temperatura de la materia Tm, se resuelven simultanea- 
mente en el codigo numerico Recfast. Los coeficientes de recombinacion en 
cada caso se parametrizan para reproducir los resultados de los analisis de 
multiples niveles realizados por Seager et al (para el hidrogeno) y Switzer 
& Hirata (para el helio), este ultimo teniendo en cuenta la aproximacion de 
Kholupenko et al para la opacidad del continuo del H sobre el He I. La ecua- 
cion de equilibrio de Saha es una buena aproximacion para la recombinacion 
del He 11. El conjunto de ecuaciones es el siguiente: 

H{z){l + z)-^ = [x^x.nnan - /3h(1 - x,)e-'^^H-/^^-) , (2.1) 



+ fxHeIlXenH«Ll " ^^^/?Hel(/He - XHeIl)e-'^He:,.3./fcT.A ^t^^^ ^ (2.2) 
V 5'HeI,lls / 

donde 

^ l + irH(AH + /5H)riH(l-Xp)' ^ • ^ 

1 + KHelAHenH(/He " XhcII ) 6^'^''- /"^^^ 
1 + irHel(AHe + Mu^^f^, - Xn,l{)e^^vs I^^M ' ^ " ^ 

^''^^ ^ 1 + nel/3Ll^H(/He - XHel^e^./'^^M " ^^.5) 

Las tres variables independientes son la fraccion de protones Xp = np/nn, 
la fraccion de helio simplemente ionizado xneii = ^Heii/^H, y la temperatura 
de la materia Tm- La variable dependiente es la fraccion de electrones libres 
Xe = ne/rin = Xp + xneii- Aqui, n se refiere a la densidad numerica, y nn 
es la densidad numerica total de hidrogeno. Los factores C tienen en cuenta 
el efecto cuello de hotella de la transicion 2p — > Is, que hace mas lenta la 
recombinacion. El segundo termino de la ecuacion (12. 2p se debe a la consi- 
der acion de las transiciones semi-prohibidas 2^p — > l^s del Hel como canal 
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adicional para la recombinacion al nivel fundamental. El supramdice t indica 
que se trata de los estados tripletes. Esta es una de las actualizaciones en el 
escenario de recombinacion. La ecuacion para la temperatura de la materia 
es 

HizKl + = '"f^^^^^;- {Tm - T^) + 2T,,H{z) . (2.6) 

dZ SrUeC 1 + /He + Xe 

Se utiliza esta ecuacion durante todo el calculo, debido al efecto que resulta 
de la diferencia entre Tm y Tr a bajo redsliift. 

A continuacion se listan todos los parametros que entran en las ecuaciones 
de mas arriba. Estos son la constante de Boltzmann k, la constante de Plan- 
ck h, la velocidad de la luz c, la seccion eficaz de Thomson ctt, la masa del 
electron rrie, y la constante de radiacion an. Con respecto a los datos atomi- 
cos, la longitud de onda en reposo de la Imea H Ly a es Ah2p = 121,5682 nm. 
La frecuencia de la Imea H 2s-ls, = c/\u2p, esta lo suficientemente cerca 
de Lya como para tomar el mismo valor de longitud de onda. La longitud 
de onda de Hel 2^p-l^s es AHei2ip = 58,4334 nm. Notar que la frecuencia 
de Hel 2^s-l^s es ^'He72s = c/60,1404nm. A diferencia de lo que pasa con 
el hidrogeno, la separacion de los niveles Hel 2^p y 2^s es suficientemente 
grande como para distinguir entre Ahc21p y ^Ue2^s, por lo tanto aparece un 
factor exponencial extra en las ecuaciones (12.41) y (12. 5p . que no aparecia en 
la ecuacion (12. 3p . con i^nei2^p2^s = ^Uei2^p ~ ^Uei2^s = ^ps, J equivalentemente 
para el termino que considera la transicion via tripletes. L a tasa de la tran - 
sicion 2s-ls via dos fotones para el H es Ah = 8.22 458 s"^ ffGoldmanl [l989|), 



mientras para el Hel es Ahc = 51,3 s"^ flPrake et al. l ll969f ). 



En las ecuaciones (12.11) y (12. 2p aparece el coeficiente de recombinacion 
Caso B, que es la suma de los coeficientes de recombinacion a todos los 
niveles e xcepto el nivel f undamental. Para el H, el coeficiente an fue calcu- 



lado por iHummed (119941 ). Efectivamente, los terminos de dicha suma estan 



constituidos por las tasas de recombinacion al nivel n: 

an{Te, Z) = ^Z^X^^n-^ / de{l + nhfe-^'dn{Z, e) (2.7) 

Jo 

donde a es la constante de estructura fina, dn{Z,e) es la seccion eficaz de 
fotoionizacion para el nivel n, e es la energfa del electron eyectado en uni- 
dades de Z^Ryd y A = Z^Ryd/Tg. El ajuste de esta cantidad, realizado por 
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Pequignot et al, 



fll99l[ ). es 



an = FlO 



-1Q S —1 



1 + ct'^ 



f2.8) 



con a = 4,309, 6 = -0,6166, c = 0,6703, = 0,5300 y t = Tm/10^K. El 
factor F vale 1.14, y permite que la ecuacion (12. ip coincida con el calculo 
obtenido con el codigo de niveles multiples, acelerando la recombinacion. 

Para el helio, el coeficient e de recombinacion Caso B, aHe, fue obtenido 
por iHummer &: Storeyl (Il998l ) , siendo la expresion de los terminos individua- 
les de la suma: 



a(T; nlS) 



^"-(2/ + l)fA^/V' 

'71 4 



dE{l + v^EYe-''^a{nlS] E) (2.9) 



donde a es la constante de estructura fina, a{nlS; E) es la seccion eficaz de 
fotoionizacion para el estado nlS a una energia del electron eyectado de E 
(en Ryd), v es el numero cuantico efectivo del estado, T es la temperatura 
electronica, y 



A 



1,5789 X 10^ 
T 



(2.10) 



El ajuste de este coeficiente se debe a IVerner fc Ferlandl (119961 ) y tiene la 
siguiente expresion: 




i+p 



ttHel = q 



Los parametros son q = 10 



-1 



-16,744 



m^s ^ 



(2.11) 



p = 0,711, Ti = y T2 fijado 

arbitrariamente a 3K. Este ajuste es bueno hasta < 0,1% sobre el rango 
de temperaturas relevante (4,000-10,000 K), y sigue siendo bastante preciso 
sobre un rango de temperatura aiin mayor. Para los tripletes, el ajuste se 
realiza con la misma forma funcional, pero con distintos valores para los 
parametros: p = 0,761; q = lO-^^.soe. ^ lO^-i^^K; and T2 = 3K. Este 
ajuste es mejor que el 1 % para temperaturas de entre 10^'^ y 10^ K. 

Los coeficientes de fotoionizacion /3x, donde X representa a H o Hel, 
se calculan a partir de los coeficientes de recombinacion haciendo: /3x = 
axi^TTniek Tm/Ii'^Y^'^ exp(—hh'2s/kTM)- Aqui J a son diferentes para H y 
Hel. Notar ademas que se utilizan Tm y i^2s, J que el uso incorrecto de Tr 
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o z/2p causara una diferencia pequefia pero importante en modelos con gran 
cantidad de bariones. Por ultimo, f3^^i es el coeficiente de fotoionizacion para 
los tripletes, y se calcula a partir de a^ei de esta manera 



/3, 



Hel 



a 



Hel 



M 



/l2 



3/2 



5'HeI,23s 



/kTy_ 



(2.12) 



donde (7He+ y 5'Hei,23s son las degeneraciones de los atomos de He^ y de He I 
con electrones en el ni vel 2^s, y hv2'is,c es la energia de ionizacion del estado 
2^8 f lWong et al.ll2008h . 

Los parametros cosmologicos son el redshift el factor de Hubble H{z), 
y la temperatura de la radiacion Tr = To(l + -z). La ab undancia primordial 
de H e se toma como Yp = 0,24 por unidad de masa (jSchramm fc Turner 
19981 ). y el valor actual de la temperatura del fondo cosmico de radiacion 
To se to ma como 2,728 K (el valor central determinado por el experimento 
FIRAS fiFixsen et al.lll996h l 

El corrimiento al rojo cosmologico de los fotones H Lya es 



(2.13) 



Para el helio, relacionamos los factores K de corrimiento al rojo de las Imeas 
Hel 2^p-l^s y Hel 2^p-l^s con la probabilidad de escape de Sobolev ps de 
la siguiente manera: 



Hel 



1 



Hel 



5'HeI,2ip ^HeI,lis^2ip-lis^'S 
S'Hel.lis 1 
5'HeI,23p '^Hel.lis^fsp.ilsPS 



(2.14) 
(2.15) 



donde AHei,2ip-iis y ^Hei,23p-iis SOU los cocficientes A de Einstein de las tran- 
siciones Hel 2^p-l^s y Hel 2'^p-l^s, respectivamente. Para incluir el efecto 
de la opacidad en el c o ntinu o debido al H, basados en la formula aproximada 
de iKholupenko et al.l (120071 ). se reemplaza ps Por la nueva probabilidad de 

escape pesc = Ps+ Pcon,H con 



Pcon.H 



1 + anel^^" 



(2.16) 
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donde ane y &He son parametros de ajuste. El factor 7 tiene la siguiente 
expresion 

^ _ gHcI,2lp ^ P ^ ^ ,^ 

87r3/2o-H,ls(^'HeI,2ip)'^HeI,2lp^'^D.2lp(l " ^^p) 

donde o"H,is(^'Hei,2ip) es la seccion eficaz de ionizacion del H a la frecuencia 

i^Hei,2ip y ^^^0,21? = '^Hei,2ip \/2/cB^M/''^HeC^ Gs el anclio termico de la Imea 
He I 2^p-l^s. La seccion eficaz de fotoionizacion desde el nivel n es ( ISeaton 



19591) 

a„(Z, h.) = ^^1^(1 + nhr^gnin, e) (2.18) 

donde gii{n,e) ~ 1 es el factor de Gaunt-Kramers, a es la constante de 
estructura fina j ao = h/ {nieca) es el radio de Bohr. 



2.2. Dependencia con las const antes funda- 
mentales 

Las constantes fundamentales a y rUe estan involucradas en cantidades 
fisicas relevantes durante la recombinacion. En esta seccion estudiamos los 
detalles de las dependencias con las constantes en cada una de las cantidades 
fisicas de interes. Luego analizaremos como cambia la liistoria de recombina- 
cion por variaciones temporales en los valores de las constantes fundamentales 
que estamos estudiando. Los primeros trabajos donde se analizo la dependen- 
cia de las constantes fundamentales en el escenario de recombinacion estandar 



son lo s de iKaplinghat et al.l (119991 ) , iHannestadl (119991 ) , y iKujat &: Scherrer 
(l200d ). 



La dependencia de los coeficientes de recombinacion Caso B puede ver- 
se de las ecuaciones (12.71) y (12. 9p . La dependencia del factor A aparece a 
traves del Ryd, una unidad de energia cuya dependencia con las constan- 
tes fundamentales coincide con la de la energia de ligadura del hidrogeno 
Bi = ^a^meC^, es decir a^mg. De la seccion eficaz de fotoionizacion viene 
una dependencia de aa^ ~ a~^m~^, que se deduce de la expresion de la 
ecuacion (I2.18p . Por lo tanto, la dependencia total de los coeficientes de re- 
combinacion Caso B con las constantes fundamentales es a^me^^"^. Por su 
parte, los coeficientes de fotoionizacion P tienen una dependencia con las 
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constantes que se deduce de la dependencia de los coeficientes de recombina- 
cion mas el factor extra m^J"^ , y la dependencia en el exponencial que va como 
Q^mei como todos los terminos proporcionales a la energia que aparecen en 
las ecuaciones. 

La dependencia con las constantes fundamentals que viene del redsliift 
de los fotones Lya se deduce de la expresion = X^/8TrH{z), con Aq, = 
Snhc/SBi, es decir, a~^m~^. 

En esta tesis se estudio la influencia de las constantes en el escenario ac- 
tualizado de la recombinacion. Para ello, buscamos la dependencia que viene 
de los factores de corrimiento al rojo de las Ifneas del helio, analizando los 
factores de las ecuaciones fl2.14p y fl2.15l) . En primera instancia, buscamos 
como dependen los coeficientes A de Einstein de las constantes fundamen- 
tales. Las tasas de pro babilidad de transicion ^Hei,2ip-iis y ^Hci,23p-iis se 



pueden expresar como (jPrake &: Mortonll2007l ): 



^] = ^^imn+r2\^P,)\' (2.19) 

donde Uij es la frecuencia de la transicion, j hace referenda a los estados 
inicial (final) del atomo. Primero analizamos el braket. A primer orden en 
teon'a de perturbaciones, todas las funciones de onda se pueden aproximar a 
la correspondiente funcion de onda del hidrogeno. Estas ultimas se pueden 
analizar como P{r/ao) * exp(— gr/ao) en donde P es un polinomio, oq es el 
radio de Bohr y g es un mimero. Se puede mostrar que cualquier integral 
del tipo de la ecuacion (12.191) puede ser resuelta con un cambio de variables 
X = r/at). Si las funciones de onda estan apropiadamente normalizadas, la 
dependencia con las constantes fundamentales viene del operador, en este 
caso ri + r2. Por ende, la dependencia del braket va como Qq. Por otra parte, 
Uij es proporcional a la diferencia de los niveles de energia y por lo tanto, 
dependen de las constantes fundamentales como uJij ~ mgO;^. En consecuen- 
cia, la dependencia de las probabilidades de transicion de Hel con a y me se 
puede expresar como: 

Affj ~ a V (2.20) 

En segunda instancia, buscamos la dependencia que viene de la modifica- 
cion de la probabilidad de escape. En la expresion de 7 (ver ecuacion fl2.17p ) 
ademas de los coeficientes A de Einstein, aparece la seccion eficaz de fotoio- 
nizacion del hidrogeno, cuya dependencia con las constantes es a~^m~'^ como 
vimos anteriormente. Solo queda analizar la dependencia de los factores ane 
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y ^He que aparecen en Pcon,H- El calculo detallado de estos parametros de 
ajuste no esta disponible aiin, por lo que n o es posible determ inar como se 
ven afectados por variaciones de a y/o vde- Wong et al. ( 2008 1) han mostra- 



do que estos parametros deben conocerse con un 1 % de precision para los 
datos que se obtendran con el satelite Planck. Sin embargo, en esta tesis se 
ha trabajado con los datos del satelite WMAP, para el cual esta precision 
no es necesaria. Para Uegar a esta conclusion hemos calculado los espectros 
de potencia de la temperatura, la polarizacion y la correlacion cruzada del 
fondo cosmico de radiacion, permitiendo que los parametros ohg y &He van'en 
un 50 %. Hemos encontrado que para los espectros de temperatura y pola- 
rizacion la variacion es siempre menor que el error observacional ( 1 % para 
la temperatura y casi 40% para la polarizacion). Las mayores variaciones 
aparecen en el espectro de correlacion cruzada. Para este, hemos calculado 
los errores observacionales para todos los ^'s medidos y los hemos comparado 
con la variacion inducida cuando cambiamos one y ^He en un 50 %. En todos 
los casos, el error observacional es varios ordenes de magnitud mayor que 
la variacion de los C^'s por el cambio en los parametros. Por lo tanto, para 
analizar los datos de WMAP no es necesario la modificacion de los mismos. 
Todo este analisis y co nclusiones son contrib uciones originales de esta tesis 



y fueron publicadas en IScoccola et al.l (l2008al ) . 



Para encontrar la dependencia de las tasas de recombinacion via dos foto- 
nes, Ah y Ahc, con las constantes fundamentals, consideremos una transicion 
donde uno de los fotones tiene una frecuencia y J el otro tiene frecuencia 
(1 — ?/)z/i2, siendo la frecuenci a correspondiente a la diferencia de energia 
entre el nivel 2s y el nivel Is (iBreit fc Teller Il940l : ISpitzer fc GreensteinI 



La probabilidad de emision de un foton de frecuencia z/12 es 



My) = ^^(^) (2-21) 

donde a es la constante de estructura fina, R es el Rydberg (IRyd = ^meC^a^), 
y ipiy) es una combinacion lineal de funciones de onda normalizadas. De 
aqui se ve que dichas tasas de recombinacion dependen de las constantes 
fundamentals como a^rrie. 

En la Tabla (12. ip . se presenta un resumen de las dependencias con las 
constantes a y mg, para las cantidades fisicas relevantes en los procesos que 
ocurren durante la recombinacion. 

La variacion temporal de las constantes fundamentals afecta la fisica 
durante la recombinacion. El principal efecto es el corrimiento de la epoca de 
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Descripcion 


Magnitud ffsica 


Dependencia 


Energia de ligadura del hidrogeno 




2 

a me 


Seccion eficaz de fotoionizacion 




— 1 —2 

a Trig 


Seccion eficaz de Thomson 






Coeficientes de recombinacion Caso B 


«H, «HeI, «HcI 


3 -3/2 
a TTLe 


Coeficientes de fotoionizacion 


/5h, /^Hel 




Redsliift de los fotones Lya 


i^H, ^Hel, ^Hel 




Coeficiente A de Einstein 


/I Hel 




Tasa de recombinacion via 2 fotones 


Ah, Ahbi 





Tabla 2.1: Dependencias con a y mg de las magnitudes ffsicas relevantes en 
los procesos que tienen lugar durante la recombinacion. 



recombinacion a redsliifts mayores cuando aumenta a o mg. Esto se entiende 
facilmente dado que la energfa de ligadura Bn escalea como a^mg, de manera 
tal que los fotones deberian tener mayor energfa para poder ionizar atomos 
de hidrogeno. En las Figs. (12.11) se muestra como se ve afectada la historia de 
ionizacion por cambios en a y en mg, en un Universo piano con parametros 
cosmologicos {Vtbh^^QcDMh'^^h^T) = (0,0223; 0,1047; 0,73; 0,09). Cuando a 
y/o rrie tienen valores mayores que los actuales, la recombinacion ocurre mas 
temprano (a redshifts mas grandes). La historia de ionizacion es mas sensible 
a variaciones en a que a variaciones en mg por las dependencias en Bn- 

El mecanismo mas eficiente para alcanzar el equilibrio termico en el gas 
de fotones del Universo temprano es el scattering de Thomson por electrones 
libres, por lo tanto, otro efecto importante es el corrimiento en la seccion 
eficaz de Thomson ctt, definida segiin: 

= ^r^" (2.22) 

3 mgC"^ 

de donde se ve su dependencia con las constantes fundamentales. 

La funcion visibilidad, la cual mide la probabilidad diferencial de que un 
foton se disperse por ultima vez al tiempo conforme r], depende de a y de 
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Est a funcion se define como 



giv) = e 



dn 
drj 



donde 



dn 

— = XeUpaaT 



(2.23) 



es la profundidad optica diferencial de los fotones debido a la dispersion de 
Thomson. La densidad numerica total de protones es Up (tanto libres como 
ligados), Xe es la fraccion de electrones libres, y a es el factor de escala. 
El efecto mas fuerte de variar a y rrie sobre la funcion visibilidad ocurre 
debido a la alteracion de la historia de ionizacion Xe{ri). En la Fig. (12. 2p se 
muestra que si el valor de a o me fuera mas pequeiio (grande) en la epoca 
de recombinacion que su valor actual, el maximo en la funcion visibilidad se 
correria a valores de redshift mas pequenos (grandes), y su ancho creceria 
(decreceria) levemente. 
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2000 1500 1000 500 



Figura 2.1: Fraccion de ionizacion como funcion del redshift, para distintos 
valores de a (panel superior) y nie (panel inferior) en la epoca de recombina- 
cion. ao y ('Tie)o son los valores actuales de las constantes. Si a o me tenfan 
en la epoca de recombinacion un valor mayor que el actual, la recombinacion 
ocurre mas temprano. 
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Figura 2.2: Funcion visibilidad como funcion del tiempo conforme, para di- 
ferentes valores de a (panel superior) y me (panel inferior). Si a o rrie tenfan 
en la epoca de recombinacion un valor mayor que el actual, la recombinacion 
ocurre mas temprano, es decir, a 7] menor. 
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Capitulo 3 



Fondo cosmico de radiacion 



3.1. Consideraciones generales 



El Universo primordial esta compuesto por un plasma de fotones, electro- 
nes y micleos de elementos qufmicos livianos. A medida que el Universo se 
expande, la temperatura baja, y es posible la formacion de atomos neutros. 
La fraccion de electrones libres se hace practicamente despreciable, por lo 
que la dispersion de Thomson entre fotones y electrones deja de tener lugar. 
La opacidad del medio cae bruscamente, y los fotones comienzan a viajar 
libremente. Estos fotones constituyen el fondo cosmico de radiacion (CMB 
por cosmic microwave background) , y traen informacion directa de como era 
el Universo en el momento del desacople de la materia y la radiacion, apro- 
ximadamente 400.000 afios despues del Big Bang. 

Los fotones del fondo cosmico de radiacion fueron dete ctados por pri- 
mera vez de manera accidental por Wilson &: Penziai ( 1967 ). Mas adelante, 
el instrumento FIRAS a bordo del satelite COBE (1992) confirmo que di- 
cha radia cion tiene un espec tro de cuerpo negro, con una temperatura de 
2.728 K (IMather et al.lll994l ). Ademas, con el instrumento DMR (differen- 
tial microwave radiom eter) se obtuv o el primer mapa de las anisotropias en 
CMB en todo el cielo (1Smootlll993l ). Luego de COBE, distintos experimen- 
tos en tierra y en globos fueron lanzados para medir con naayor precision 



la distribuc ion de dichas anisotro pias: A RCHEOPS flBenoitM2004f) . BOO- 



MERANG fiNetterfield et al.ll2002h. DASI flHalverson et al.ll2002h . MAXIMA 
fiLee et al.lll999h. VSA flCrainge et al.ll2003h~ y con m ejor resolucion angular: 
CBI flMason et al.ll200lh . v ACBAR (IKuo et al.ll2002h . En Junio de 2001, NA- 
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SA lanzo el satelite WMAP, que en el ano 2003 libero un primer conjunto de 
datos que constituyero n un mapa completo d el cielo con una precision muclio 
mayor a la de COBE flHinshaw et al.ll2003h . En el 2005, WMAP libero un 



conjunto de datos mayor, que Uamamos WMAP3 (jSpergel et al.l 120071 ) . y 



en marzo de 2008 publico los datos recolectad os durante 5 anos, conjunto 
que Uamaremos WMAP5 (jPunkley et al.ll2009l ). La Agenda espacial euro- 
pea (ESA) planea lanzar el satelite Planck durante el ano 2009, para realizar 
el mapa definitivo de las anisotropias en la temperatura y polarizacion de 
CMB en todo el cielo , con una resolucion angular de un decimo de grado 
teurigana et all 12004 ). 

La distribucion de la temperatura de CMB en el cielo es una funcion 
definida sobre una esfera, por lo tanto, es natural analizarla con un desarrollo 
en armonicos esfericos 



T(e,0) = ^a,^y,„(^,0). (3.1) 

En esta combinacion lineal, el termino del monopolo da la temperatura 
media de CMB, T = 2,728K. El termino dipolar, con £ = 1, se interpreta 
como el resultado del corrimiento Doppler causado por el movimiento del 
sistema solar relativo a CMB. Una vez que se remueven estos dos terminos, 
lo que queda son las anisotropias intrinsecas de CMB, que son del orden de 
10~^ en todas las escalas angulares, y contienen la informacion de la fi'sica del 
Universo en la epoca del desacople de la materia y la radiacion. La mayor par- 
te de la informacion cosmologica esta contenida en la funcion de correlacion 
de dos puntos de la temperatura. Esta funcion se define como el promedio 
en el cielo del producto de la desviacion fraccional de la temperatura en las 
direcciones n y n', desarroUando el resultado en polinomios de Legendre 

C(e) J^^) .t?^c,p,i.o.e, (3.2) 

Los coeficientes del desarrollo, Q, son las varianzas de los coeficientes a^m, 
de manera tal que 

{O'imO'e'm') = Su'Smm'Ce. (3.3) 

Esta funcion, Uamada espectro angular de potencias de las anisotropias en el 
fondo cosmico de radiacion, es la herramienta mas importante al comparar 
teoria con observaciones. 
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Las anisotropias en la temperatura estan mtimamcntc rclacionadas con 
las fluctuaciones dc la densidad en la epoca de recombinacion. Los modelos 
inflacionarios predicen que el campo de inhomogeneidades en la densidad, 
S{x), es un campo aleatorio gaussiano. Estas inhomogeneidades son la se- 
milla para la formacion de estructura a gran escala, y se supone que fueron 
generadas a partir de fluctuaciones cuanticas en el campo del inflaton, que 
crccicron con un factor de expansion exponencial. Definiendo la desviacion 
media cuadratica de las inhomogeneidades como = {S{x)'^), la funcion de 
autocorrelacion se define como 

e(|x2-xi|)^e(^) = l(5(x,)5(x0). (3.4) 

En la aproximacion donde las fases son aleatorias, todos los modos de Fourier 
S{k) estan no correlacionados, lo que significa que 

(5*(k)5(k')) = (27r)'5'(k - k')P(k). (3.5) 

La funcion P(k) es el espectro de potencies de las fluctuaciones. En la mayon'a 
de los modelos, el espectro inicial toma la forma 

P(k)=A.[ff- (3.6) 

donde As es la amplitud y Us es el indice espectral. El valor = 1 corres- 
ponde a las Uamadas fluctuaciones invariantes de escala, dado que para tal 
espectro, las fluctuaciones S tienen la misma amplitud para todas las esca- 
las de longitud. La teoria inflacionaria predice perturbaciones que sean casi 
invariantes de escala. 

Los fotones que fueron dispersados por ultima vez en el desacople, nos 
Uegan de la epoca cuando se formaron los atomos neutros en el Universo. A 
esa superflcie en el cielo, que se encuentra a una distancia correspondiente 
al redshift del desacople, se la conoce como superflcie de ultimo scattering. 
Sin embargo, la recombinacion no fue instantanea, por lo que la superflcie 
de ultimo scattering tiene cierto espesor en redshift. La dispersion de Thom- 
son entre fotones y electrones en la superflcie de ultimo scattering genera 
polarizacion en los fotones de CMB. Por lo tanto, se pueden deflnir nuevas 
funciones de correlacion angular a partir de la polarizacion, y asi tener mas 
herramientas con las cuales extraer informacion acerca de las condiciones 
fisicas del Universo. La manera mas conveniente de describir la polarizacion 
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es a traves de sus modos E (escalar) y B (pseudo-escalar). La funcion de 
correlacion cruzada entre temperatura y polarizacion se define como 

{'^Jm'^f'm') = ^U'^mm'Cj^ (3.7) 

y la funcion de autocorrelacion de los modos E de polarizacion se define como 

{(^fm'^fim') = ^U'^mm'Cf^ (3.8) 

donde af^ son los coeficientes del desarroUo en armonicos esfericos de la pola- 
rizacion modo E. Dado que las perturbaciones iniciales escalares no generan 
modos B de polarizacion, no estudiamos las correlaciones que involucren es- 
tos modos. 

Tanto los parametros cosmologicos como el valor de las constantes fun- 
damentales a y rUe afectan la forma de las funciones de correlacion, y de 
esta manera sus valores pueden estimarse a partir de observaciones de CMB. 
Los efectos en el espectro angular de potencias de CMB debido a variar las 
constantes fundamentales son similares a los producidos por cambios en los 
parametros cosmologicos, es decir, cambios en las amplitudes relativas de los 
picos, y un corrimiento en sus posiciones. Efectivamente, el aumento de a o 
rrie aumenta el redshift de la superficie de ultimo scattering, lo que correspon- 
de a un horizonte de sonido mas pequefio. La posicion del primer pico (^i) 
es inversamente proporcional a este ultimo, por lo que resulta un £i mayor. 
Tambien se produce un efecto ISW (integrated Sachs- Wolfe) mayor, hacien- 
do que el primer pico Doppler sea mas alto. Mas aiin, un aumento en a o 
rrie decrece el damping de difusion para £'s grandes, el que se produce por el 
ancho finito de la superficie de liltimo scattering, y por lo tanto, aumenta la 
potencia en pequefias escalas. Todos estos efectos se ilustran en la Fig. fl3.1l) . 
donde ao y (me)o son los valores actuales de las constantes fundamentales a 
y rrie- En este ejemplo, los parametros cosmologicos fueron fijados a: 

l^b/i^ = 0.0223 VLcDMh'^ = 0.1047 Hq = 73 km/s/Mpc 
r^e = 0.09 = 0.951 A, = 2.3x10"^ 

donde Qbh'^ es la densidad de materia barionica, en unidades de la den- 
sidad critica, flcDMh"^ es la densidad de materia oscura, en las mismas uni- 
dades, Hq es la tasa actual de expansion de Hubble, r^e es la profundidad 
optica de reionizacion, Ug es el indice espectral de perturbaciones escalares, 
y As es la amplitud de las fluctuaciones escalares. 
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Figura 3.1: Espectro de anisotropias de la temperatura de CMB para dife- 
rentes valores de a (panel superior) y me (panel inferior). 
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La manera en la que cambia el espectro de correlacion cruzada entre tem- 
peratura y polarizacion para distintos valores de las constantes fundamenta- 
les, se ve en la Fig. (13.21) . El mismo efecto, para el espectro de potencias de 
la polarizacion, esta ilustrado en la Fig. (13. 3p . 



3.2. Estadisticas 

En general, es de interes en estadistica usar una muestra de datos para 
hacer inferencias sobre un modelo probabilfstico, por ejemplo para evaluar si 
el modelo es valido o bien para estimar los valores de sus parametros. 

Hay dos enfoques principales para la inferencia estadistica, que llamare- 
mos frecuentista y bayesiano. En la estadistica frecuentista, la probabilidad 
se interpret a como la frecuencia de aparicion de un resultado en un experi- 
mento que se repite. Las herramientas mas importantes en este enfoque son 
la estimacion de parametros, y los tests estadi'sticos. La estadistica frecuen- 
tista provee las herramientas usuales para informar objetivamente acerca del 
resultado de un experimento sin necesidad de incorporar creencias previas 
concernientes al parametro a ser medido o a la teoria a ser testeada. 

En la estadistica bayesiana, la interpret acion de probabilidad es mas ge- 
neral, e incluye el grado de creencia. Se puede entonces hablar de la funcion 
densidad de probabilidad para un parametro, lo cual expresa el grado de 
conocimiento que tenemos acerca de donde se encuentra su valor verdadero. 
Los metodos bayesianos permiten agregar de manera natural informacion ex- 
tra tal como limitaciones fisicas o informacion subjetiva. De hecho, requieren 
la imposicion de un prior para los parametros, es decir, el grado de creencia 
sobre los valores de los parametros antes de realizar la medida. Usando el 
Teorema de Bayes, el grado de creencia a priori se actualiza con los datos 
obtenidos en el experimento. 

Para muclios problemas de inferencia, los enfoques frecuentista y baye- 
siano dan los mismos resultados numericos, a pesar de estar basados en in- 
terpretaciones diferentes de la probabilidad. Sin embargo, para muestras de 
datos pequefias, o para medidas de un parametro cerca de un limite fisico, los 
distintos enfoques pueden dar resultados diferentes, y es necesario hacer una 
eleccion. Elegimos el enfoque bayesiano porque es el adoptado en la literatura 
y, por lo tanto, podremos comparar nuestros resultados con los obtenidos por 
otros autores. 
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Figura 3.2: Espectro de correlacion cruzada entre T y polarizacion modo E 
de CMB para diferentes valores de a (panel superior) y rrie (panel inferior). 
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Figura 3.3: Espectro de potencias de la polarizacion modo E de CMB para 
diferentes valores de a (panel superior) y rrie (panel inferior). 
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3.2.1. El enfoque bayesiano 



Desde la perspectiva bayesiana, no hay una distincion fundamental entre 
observables y parametros de un modelo estadistico: todos son considerados 
cantidades "desconocidas" , y se les puede asociar un concepto de probabili- 
dad. Sea D los datos observados, y sea 9 los parametros del modelo y los datos 
faltantes. Formalmente, se puede construir una distribucion de probabilidad 
conjunta P{D, 6) sobre todas las cantidades aleatorias y/o desconocidas. Es- 
ta distribucion conjunta consta de dos partes: una distribucion a priori P{6) 
y una funcion Likelihood P{D\9). La Likelihood es una funcion que expresa, 
luego de que el resultado de un experimento es conocido, que tan probable 
era obtener ese resultado. Es una funcion de los parametros y no debe con- 
fundirse con una densidad de probabilidad. Especificar P{6) y P{D\6) da un 
modelo probabilfstico completo, en el cual P(-D, 6) = P{D\6)P{6). Habiendo 
observado D, se usa el teorema de Bayes para determinar la distribucion de 
9 condicionada a D: 

^^^'^^ - j p{9)p{D\9)de ^^-^^ 

esto se llama la distribucion posteriori de 0, y es el objeto de la inferencia ba- 
yesiana. A partir de esta distribucion pueden calcularse los valores esperados 
a posteriori de cualquier funcion f{9): 

^f^^^)'^] - jPi9)PiD\9)d9 ■ ^^-^^^ 

Las integraciones en esta expresion eran la dificultad maxima de la inferencia 
bayesiana, y el metodo de las cadenas de Markov Monte Carlo (MCMC) 
surge como una manera viabl e de ey aluar este tipo de expresiones. Seguimos 



el tratamiento de iGilks et al.l (1l996l ) 



3.2.2. Cadenas de Markov Monte Carlo 

Las cadenas de Markov Monte Carlo surgen como un metodo para evaluar 
expresiones del tipo de la Ec. (13.101) y estan compuestas de dos partes: la 
integracion Monte Carlo y la cadena de Markov. 
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Integracion Monte Ccirlo 



La integracion Monte Carlo evalua E[f{X)] tomando muestras {Xt,t = 
1, . . . , n} a partir de la funcion P{9\D), y luego aproximando 



Por lo tanto, la media de la poblacion se estima con la media de la muestra. 
Cuando las muestras {Xt} son independientes, la Ley de los grandes mimeros 
de la estadistica ascgura que la aproximacion se hace cada vcz mas precisa 
al aumentar n. Es importante notar que n no es el mimero de datos, sino el 
tamafio de la muestra, que esta controlado por el analista. 

En general, tomar muestras {Xt} a partir de P{9\D) no es posible dado 
que esta funcion puede ser muy poco comun. Sin embargo, se pueden generar 
por cualquier proceso que tome muestras dentro del soporte de P{9\D), en las 
proporciones correctas. Una manera de hacer esto es a traves de una cadena 
de Markov, que tenga a P{9\D) como distribucion estacionaria. En ese caso, 
se trata de una cadena de Markov Monte Carlo. 

Cadenas de Markov 

Supongamos que generamos una secuencia de variables aleatorias 



de tal manera que, para cada t > 0, el siguiente estado Xt+i se muestree de 
una distribucion P{Xt+i\Xt) que dependa solo del estado actual de la cadena, 
Xt. Esta secuencia se llama cadena de Markov, y P{-\-) se llama el kernel de 
transicion de la cadena. La cadena ira perdiendo informacion de cual era el 
Xq inicial, y la distribucion de la cual se muestrean los valores tender a a una 
distribucion estacionaria </>(■). 

Luego de un tiempo de "burn-in" suficientemente largo, digamos de m 
iteraciones, los puntos {Xt; t = m + 1, . . . ,n} seran extraidos de la distribu- 
cion 0(-). Hay metodos para estimar el valor de m. Ademas, es necesaria una 
manera de asegurar que la cadena construida es tal que su distribucion es- 
tacionaria 0(-) coincide con la posteriori P{9\D). Una de las posibilidades es 
utilizar el algoritmo de Metropolis-Hastings. Luego de descartar las muestras 




(3.11) 



t=i 



{Xq, Xi, X2, . . .} 
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del burn-in, el estimador de E[f{X)] es: 



f = ^ E /™ (3.12) 



n — m 

t=m+l 



que se conoce como promedio ergodico. La convergencia al valor de expecta- 
cion requerido se asegura mediante el teorema ergodico. 

Algoritmo de Metropolis-Hastings 

Para este algoritmo, en cada paso t, el siguiente estado X^+i se elige 
primero muestreando un punto candidato Y, de la distribucion propuesta 
q{-\Xt), que por ejemplo, puede ser una distribucion normal de variables 
multiples, con valor medio X y una matriz de covarianza fija. Para el punto 
candidato Y se calcula una probabilidad a{Xt, Y) tal que 

A la vez, se elige un numero aleatorio entre y 1. Luego se compara dicho 
mimero con a(Xt,y); si es menor, se toma X^+i = Y, sino, se toma Xf+i = 
Xt. Finalmente se incrementa t. Cabe destacar que la distribucion propuesta 
g(-|-) puede tener cualquier forma, y la distribucion estacionaria de la cadena 
sera P{9\D). 

Convergencia 

El problema fundamental de la inferencia a partir de simulaciones de ca- 
denas de Markov es que siempre habra regiones del espacio de parametros 
que no seran cubiertas por una cadena de longitud finita. A medida que la 
simulacion avanza, la propiedad ergodica de la cadena de Markov hace que 
eventualmente se cubra toda la region, pero en el corto plazo, las simula- 
ciones no pueden decir nada acerca de regiones que la cadena aiin no haya 
explorado. Esto es un problema cuando la convergencia es lenta. La mejor 
manera de notar la convergencia lenta es examinando multiples simulaciones 
independientes. Tambien es importante que los puntos iniciale s de cada una 
de las cadenas esten suficientemente separados (lGelmanlll996l ). 

La manera recomendada para monitorear la convergencia se basa en de- 
tectar cuando las cadenas de Markov se ban "olvidado" de sus puntos de 
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inicio, comparando varias secuencias largadas desde diferentes puntos inicia- 
les y chequeando que sean indistinguibles. Una forma cuantitativa utiliza el 
analisis de las varianzas: se diagnostica convergencia aproximada cuando la 
varianza entre cadenas diferentes no es mayor que la varianza dentro de una 
secuencia individual. 

De manera independiente, se monitorea la convergencia de todos los 
"rcsumenes escalares" de intcrcs. Por cjcmplo, podemos estar interesados cn 
todos los parametros de la distribucion, y varias cantidades predichas. Cada 
uno de esos parametros se considera un resumen escalar ip. Consideraremos 
uno solo, para ilustrar el metodo. Supongamos que tenemos m simulaciones 
paralelas, cada una de longitud n. Queremos una manera cuantitativa de 
determinar si "las diferentes secuencias se apartan entre si mas de lo espe- 
rado, respecto de su variabilidad interna". Llamamos ipij con i = 1, . . . ,m 
y j = 1, ... ,n a los eslabones de cada cadena. Se define la varianza entre 
cadenas B, y \a varianza dentro de las cadenas W: 

B = — - y'(V'i - donde ^/Sj = - V ^p^, -0 = — V 

1=1 j=l i=l 



w 



^ m ^ n 

-J2s^ donde s^^ (3-14) 



A partir de estas dos cantidades, construimos dos estimadores de la va- 
rianza de ■0. El primero es 

n — 1 1 

vB^m^ W + -B (3.15) 

n n 

es un estimador de la varianza "no-sesgado" si los puntos iniciales de las 
simulaciones fueron verdaderamente extraidos de la distribucion P{6\D), pero 
es "sobre-estimado" si los puntos iniciales estan apropiadamente dispersados. 

Por otra parte, para cualquier n finito, la varianza dentro de la secuencia, 
W , debe subestimar la varianza de ip porque las secuencias individuales no 
han tenido tiempo de recorrer toda la distribucion, y como resultado tendran 
menos variabilidad. En el limite n — > oo tanto var{'ip) como W se acercan a 
var{i/j) pero desde direcciones opuestas. 

Existe un estimador R definido como 



R-Xl^ (3.16) 
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que a medida que la secuencia converge, va acercandose a 1, reflejando que 
las cadenas de Markov se estan superponiendo. Se considera que la cadena ha 
convergido cuando i? es 1,2 o 1,1. Una practica mas conservadora es tomar 
i?-l<0,03. 

3.2.3. Funcion Likelihood para CMB 

Las teorfas mas simples que explican la generacion de inhomogeneida- 
des en el Uni verso, tales como la inflacionaria, predicen que la serial de la 
temperatura en cierta direccion del cielo es una muestra de una distribucion 
gaussiana. A partir de allf se puede deducir la funcion Likelihood L para Np 
pixeles de datos de la anisotropia en la temperatura del fondo cosmico de 
radiacion como 

L = — ^ — exp|--AC-^Al (3.17) 

(27r)^-/'(detC)^/' 12 / ^ ^ 

donde A es el vector de datos que consiste en las Np medidas de la anisotropia 
de la temperatura en distintas direcciones en el cielo y C es la matriz de 
covarianza total, suma de la matriz de covarianza del ruido mas la matriz de 
covarianza de la serial. 



3.3. Ajustes 



Uno de los objetivos principales de esta tesis fue encontrar llmites a la 
variacion de las constantes fundamentales entre la epoca de la recombinacion 
del Universo y la actualidad. Hemos estudiado los casos donde solo una de 
las constantes fundamentales varia con el tiempo (ya sea a o rrie), y el caso 
donde la variacion es simultanea, aunque independiente. 

Para el analisis estadistico, hemos utilizado distintos conjuntos de datos y 
comparado como cambian las cotas obtenidas en cada caso. Asimismo, coteja- 
mos nuestros resultados con los disponibles en la literatura. Hemos utilizado 
el conjunto mas completo de datos de las anisotropi'as en la temperatura y 
polari zacion del Fondo Cosm ico de Ra diacion, que incluye los resultados de 
CBI flReadhead et al.l liool. ACBAR (IKuo et al.l liooi). y BOOMERANG 



(jPiacentini et al.ll2006l : I Jones et al.ll2006l ). Utilizamos tambien los datos del 



satelite WMAP. Para estos, tenemos dos conjuntos de datos, que fueron libe- 
rados en distintos momentos durante la realizacion de este trabajo de tesis: 
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los recolectados durante los primeros tres anos de funci onamiento del satel i- 
te, liberados en el afio 2005, que llamaremos WMAP3 fepergel et aDl2007l) 



y los publicados en el ano 2008, que llamaremos WMAP5 (iDunkley et al 



20091 ). y son el resultado de cinco anos de observacion. 



Dada la fuerte degeneracion que existe entre la variacion de las constan- 
tes y el factor de expansion de Hubble Hq, hemos realizado un ajuste con 
los datos de CMB y un prior gaussiano para Hp, el cual fue obtenid o del 
Proyecto Clave del Telescopio espacial Hubble (IFreedman et al.l l200ll ) . No 
utilizamos el ajuste final para Hq que presentan en dicho trabajo, dado que 
algunos de los objetos utilizados estan a alto redshift. En cambio, utilizamos 
los valores obtenidos a partir de los objetos celestes mas cercanos, de manera 
tal de poder despreciar cualquier posible diferencia entre el valor de a y me 
a ese redshift determinado y el valor actual. El valor medio de la gaussiana 
considerada en esta tesis es de 71 km/s/Mpc, y la desviacion estandar es 
de 10 km/s/Mpc, teniendo en cuenta errores a leatorios y sistematicos . Men- 
cionamos que el valor final que presentan en iFreedman et al.l (120011 ) es de 
i/o = 72 ± 8 km/s/Mpc. 

El espectro de potencias de las galaxias (o de la materia) P{k) es creciente 
para valores pequefios de k (grandes escalas) y decreciente para k grandes. La 
escala para la cual se presenta el quiebre corresponde a aquella que entra al 
horizonte en la epoca de igualdad de materia y radiacion. El factor de escala 
correspondiente es aeq = 4, 15 x 10~^ /Q^h"^ ■ Esta dependencia con Qmh"^ nos 
permite romper la degeneracion que existe con a y rrie, y obtener cotas mas 
restrictivas para sus variaciones. Por lo tanto, hemos utilizado los datos del 
espectro de pote ncias de la distribucion de galaxias del catalogo 2dFGRS 
(ICole et al.ll2005l ). junto con los de CMB para obtener cotas a la variacion 
de las constantes fundamentales. 

Clasificamos los ajustes realizados segiin el conjunto de datos que se haya 
utilizado: 



■ Todas las observaciones de CMB, incluyendo los datos de WMAP3. 

■ Todas las observaciones de CMB, incluyendo los datos de WMAP3, 
junto con el prior para Hq del HST Key Project, obtenido a partir de 
objetos cercanos. 

■ Todas las observaciones de CMB, incluyendo los datos de WMAP3, 
junto con el espectro de potencias P{k) del catalogo 2dFGRS. 
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■ Todas las observaciones de CMB, incluyendo los datos de WMAP5, 
junto con el espectro de potencias P{k) del catalogo 2dFGRS. 

En los analisis estadisticos, hemos considerado un modelo cosmologico 
piano con fluctuaciones en densidad adiabaticas, con el siguiente conjunto de 
parametros cosmologicos, que tambien son ajustados: 



(3.18) 



donde Qbh"^ es la densidad de materia barionica, en unidades de la densidad 
crftica, VtcDuh'^ es la densidad de materia oscura, en las mismas unidades, 9 
da el cociente entre el horizonte comovil de sonido en la epoca de desacople y 
la distancia angular a la superficie de ultimo scattering, r^e es la profundidad 
optica de reionizacion, es el mdice espectral de perturbaciones escalares, 
y As es la amplitud de las fluctuaciones escalares. 

Para llevar a cabo el analisis estadfstico, utilizamos el metodo de las ca- 
denas de Markov para explorar el espacio de parametros porque el uso de 
grillas cuando el numero de parametros es tan grande demanda un tiem- 
po de computo excesiva mente largo. Estas tecn icas estan implementadas en 
el codig o CosmoMC de iLewis fc Bridld (120021 ). El mismo, utiliza el codigo 
CAMB (ILewis et al.ll2000l ) par a calcular los espectr os de potencias angulares 
de CMB, y el codigo Recfast (ISeager et al.lll999bl ) para resolver las ecuacio- 
nes de recombinacion. Modificamos dichos programas para incluir variaciones 
en el valor de a y de me durante la recombinacion. En estos a nalisis, hemos 



considerado l a recombinacion estandar, segiin el tratamiento de ISeager et al. 



( 2000l . 1999a ). Hemos realizado el analisis variando las constantes por se- 
parado, y en forma conjunta, y los resultados se muestran en las secciones 
subsiguientes. Se ban considerado priors pianos para todos parametros (ex- 
cepto cuando se tiene en cuenta los datos del HST, en cuyo caso se toma 
un prior gaussiano para Hq, como se menciono anteriormente) . Para cada 
uno de los analisis llevados a cabo, se ban co rrido 8 cadenas d e Mar kov, y 
se ha utilizado el criterio de convergencia de iRaftery fc Lewis (Il992h para 
detener las cadenas cuando el estimador _R — 1 toma un valor bajo (que espe- 
cificaremos en cada caso). En general hemos trabajado con cadenas largas ya 
que, para asegurar la convergencia, no hemos cortado el proceso hasta que el 
valor del estimador R — 1 fuera mas bajo que los valores que se consideran 
habitualmente. 

Finalmente, comentamos que Planck sera, la primera mision en realizar un 
mapa de todo el cielo con una sensibilidad del orden del mJy y una resolucion 
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mejor que 10' ( Villa et al.l 2003 ). Tal resolucion permitira ver dentro de la 
cola de decaimiento (damping) del espectro de anisotropfas, cerca del tercer 
y cuarto pico, con un a precision limitada exclusivamente por la varianza 
cosmica ( Whitel 2006 ). Esto permitira una determinacion muy precisa del 
cociente barion-foton a partir de la altura relativa de los picos del espectro. 



3.3.1. Variacion de a 



Para poner cotas a la variacion temporal de a durante la recombinacion, 
al conjunto de parametros fl3.18p le agregamos Aa/ao, que tiene en cuenta 
posibles variaciones en el valor de la constante de estructura fina respecto de 
su valor actual a^. 

Existen algunos analisis previos de los datos de CMB incluyendo una po- 
sible variacion de la constante de estructura fina a. Los primeros trabajos que 
presentan cotas a partir de datos del fondo cosmico de radiacion son previos 
al satelite WMAP. Ajustes a la variacio n de a usand o dato s de BOOME- 
RANG y MAXIMA son presentados en lAvelino et aP fl2000h 
fl200lh . y usando BOOMERANG , COBEy DASI en [ Xvelino et alj fl200lf ) 



En los trabajos de Martins et al.1 ( 2002 ): Rocha et al. ( 20031 ) se utilizan los 
datos de WMAPl (recolectados durante el primer afio de funcionamiento del 
satelite) y se obtiene 

■ -0,06 < A«/ao < 0,01 



con un 95 % de confianza. Por otra parte, en llchikawa et al.l (120061) utilizan 
los datos de WMAPl y el prior gaussiano para i^o dado por el HST Key 
Project (i^o = 72 ± 8 km/s/Mpc). El metodo estadistico que utilizan es el 
calculo del mmimo del como funcion de a. Para minimizar sobre los otros 
parametros cosmologicos, utilizan iterativamente el metodo de Brent, de la 
interpolacion sucesiva por parabolas. Las cotas que obtienen son 

■ -0,107 < Aa/ao < 0,043, usando los datos de WMAPl, y 

■ —0,048 < Aa/ao < 0,032, cuando se agrega el prior de HST. 



Stefanescu ( IStefanescul 120071 ) encontro limites a la variacion de a utilizando 
los datos de WMAP3 junto con el prior para Hq dado por el HST Key Project. 
Por otra parte, en diclio trabajo se desprecia el efecto de la variacion de a 
sobre la recombinacion del lielio. El li'mite presentado por esta autora es 
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-0,039 < A«/ao < 0,010, 



con un 95 % de confianza. Por ultimo, iNakashima et al.l (120081 ) utilizan los 
datos de WMAP5 para obtener cotas sobre la variacion de a. Los resultados 
que obtienen son 

■ —0,028 < Aa/ao < 0,026, si se usa el prior del HST, y 

■ —0,050 < Aa/ao < 0,042, si no se lo utiliza. 

Resaltamos que estos autores tambien usan el resultado final para Hq, cal- 
culado a partir de objetos a distintos redsliifts. 

En nuestro analisis considerando solo datos de WMAP3 encontramos la 
siguiente cota 

- 0,065 < Aa/ao < 0,019 (3.19) 

con una confianza del 95 %. Para estas cadenas, el estimador de convergencia 
vale R — 1 = 0,0266. Cuando incorporamos el prior para Hq del HST, calcu- 
lado a partir de objetos cercanos {Hq = 71 ± 10 km/s/Mpc), el intervalo de 
confianza se hace mas pequeiio 

- 0,047 < Aa/tto < 0,014. (3.20) 
La diferencia entre nuestro resultado y el presentado por lStefanescu fl2007h se 



debe al prior elegido para Hq. En esta tesis preferimos considerar el prior obte- 
nido a partir de objetos cercanos. Si a las cadenas correspondientes al analisis 
realizado con los datos de WMAP3 se les realiza un post-procesamiento para 
agregar el prior para Hq dado por el valor final del HST Key Project, el 
intervalo del 95 % de confiabilidad que obtenemos es 

- 0,040 < Aa/ao < 0,014. (3.21) 

Tanto Stefanescu como Nakasliima et al, no ban considerado la posibilidad 
de que la cota sobre Ho pueda verse afectada si a tiene un valor diferente al 
actual en las epocas correspondientes a los redsliifts de los objetos utilizados 
en el HST Key Project. 

Las cotas mas fuertes se obtienen al combinar los datos de WMAP con 
los del P{k) del 2dFGRS (ver Figs. ([32D y fl33D l 

Cuando utilizamos los datos de WMAP3 y el espectro del 2dFGRS, corta- 
mos las cadenas cuando i? — 1 < 0,003. En este caso, el intervalo de confianza 
del 95 % es 

- 0,038 < Aa/ao < 0,009. (3.22) 
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Figura 3.4: Contornos de confianza del 68% y 95% para los parametros 
Aa/ao y flbh"^ obtenidos con datos de WMAP3, con datos de WMAP3 mas 
el espectro de potencias P{k) del 2dFGRS, y con WMAP3 y el prior para 
Hq calculado a partir de objetos cercanos, para el analisis con variacion de 
a. 



Estos resultados estan publicados en iMosquera et al.l (|2008[ ) . 

Finalmente, realizamos el ajuste con los datos de WMAP5 y 2dFGRS. 
Las cadenas se cortaron cuando R—l = 0,0213. En este caso, el intervalo de 
confianza del 95 % es 



0,019 < Aa/ao < 0,017. 



(3.23) 



Esta cota es mas fuerte que la presentada en iNakashima et al.l (120081 ). Efec- 
tivamente, en las Figs. (13.41) y (13.51) se puede ver que los contornos son mas 
pequenos cuando se agregan a CMB los datos del P{k) que cuando se agrega 
el prior para Hq. 

En la Fig. (13.61) se muestran las distribuciones a posteriori marginalizadas, 
para el caso de variacion de a obtenidas con los datos de CMB, incluyendo 
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CMB 




Figura 3.5: Contornos de confianza del 68 % y 95 % para los parametros 
Aa/tto y Hq obtenidos con datos de WMAP3, con y sin datos del espectro 
de potencias del 2dFGRS, para el analisis con variacion de a. 



los datos de WMAP5 y el espectro de potencias del 2dFGRS. En la diagonal 
se muestran las distribuciones de densidad de probabilidad a posteriori para 
los parametros individuales indicados en el eje x. Estas distribuciones estan 
marginalizadas sobre los otros parametros y normalizadas a 1 en el maximo. 
Los otros paneles muestran los contornos de confianza para pares de parame- 
tros. Los contornos representan 68 %, 95 % y 99 % de confiabilidad, desde 
el mas interno al mas externo, respectivamente. En dicha figura se puede 
ver una degeneracion fuerte entre a y G, esta ultima directamente relacio- 
nada con Hq, y tambien entre a y Qbh"^. Los valores obtenidos para Qbh"^, 
h, flcDAih'^, Tre, Y estau en acuerdo, dentro del intervalo a la, con los 
correspondientes valor es cuando no se inc luye variacion de a, obtenidos por 
el equipo de WMAP (jSpergel et al.l 120071 ). Los resultados aquf presentados 
son consistentes a 2cr con variacion nula de a en la epoca de recombinacion. 
Para el ajuste con los datos de WMAP3 se obtienen resultados similares, 
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aunque las cotas son mas precisas cuando se utilizan los datos de WMAP5. 

Los resultados de los ajustes, con WMAP3 y WMAP5, para los parame- 
tros cosmologicos y Aa/a^ se muestran en la Tabla fl3.ip . A modo de com- 



Parametro 


WMAP3 + P{k) 


WMAP5 ± P{k) 




0,0216 ±0,0009 


0,0226 ± 0,0007 


^CDIVlh'^ 


0,102 ±0,006 


0,110 ±0,005 


e 


1,021 ±0,017 


1,040 ±0,013 




0,092 ±0,014 


0,086:^o;ooo8 




-0,015 ±0,012 


-0,002 ± 0,009 


ris 


0,965 ±0,016 


0,964 ±0,014 


log (101%) 


o Qoq+0,064 

0,UOJ_g Qg5 


3,060 ±0,037 


Ho 




71 7+4,0 
' -L) ' -4,1 



Tabla 3.1: Valores medios y errores a la para los parametros incluyendo la 
variacion de a, para el analisis con WMAP3 y el analisis con WMAP5. Las 
unidades de Hq son km s"^ Mpc"^. 

paracion, en la Tabla (13.21) se muestran los ajustes para los parametros cos- 
mologicos, publicados por el equipo del satelite WMAP, en el modelo cos- 
mologico canonico donde no hay variacion de las constantes. Vemos que los 
valores medios de los parametros cosmologicos practicamente no cambian si 
se considera variacion de a, a excepcion de que en caso de variar a 

prefiere valores mas bajos que en el caso canonico, para ambos conjuntos de 
datos. Otra diferencia que se observa es en los valores de Tre y -f^o, para el 
ajuste con WMAP3. En este caso, cuando hay variacion de a, se prefiere un 
valor ligeramente mas alto para r^e, y un valor mas bajo para Hq. 

En la Fig. (13. 7p se compara como cambia la distribucion de densidad de 
probabilidad para el parametro Aa/ag cuando se consideran los datos de 
WMAP3 y los de WMAP5. En ambos casos, tambien se utilizaron los datos 
del espectro de potencias del 2dFGRS. 

En la Fig. (13. 8p se resume cuales son las cotas mas actuales para la va- 
riacion de a, entre la epoca de recombinacion y la epoca actual, existentes 
en la literatura. 
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Figura 3.6: Distribuciones a posteriori marginalizadas, para el caso de varia- 
cion de a obtenidas con los datos de CMB, incluyendo los datos de WMAP5 
y el espectro de potencias del 2dFGRS. En la diagonal se muestran las dis- 
tribuciones de densidad de probabilidad a posteriori para los parametros 
individuales indicados en el eje x. Estas distribuciones estan marginalizadas 
sobre los otros parametros y normalizadas a 1 en el maximo. Los otros paneles 
muestran los contornos de confianza para pares de parametros. Los contornos 
representan 68%, 95% y 99% de confiabilidad, desde el mas interno al mas 
externo, respectivamente. Analisis estadistico con variacion de a. 
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Parametro 


WMAP3 


WMAP5 


100 Qbh^ 


2,229 ± 0,073 


2,273 ± 0,062 




0,1054 ±0,0087 


0,1099 ±0,0062 




0,089 ± 0,030 


0,087 ±0,017 


Us 


0,958 ±0,016 


0,963+J°- 


Ho 




71,9+1? 



Tabla 3.2: Valores medios y errores a la para los parametros cosmologicos a 
partir de los datos de WMAP, sin variacion de las constantes. Las unidades 
de Hq son km s~^ Mpc~^. 




.04 0.06 



Figura 3.7: Distribucion a posteriori unidimensional para Aa/aQ, comparan- 
do los ajustes con WMAP3 y WMAP5. Analisis con variacion solo de a. 

3.3.2. Variacion de iJie 

Para encontrar cotas a la variacion de m-e entre la epoca de recombinacion 
y la actualidad, al conjunto de parametros (13.181) le agregamos la variable 
Ame/(me)o. 

Existen algunos trabajos previos q ue estudian la var i acion de la masa 
del electron usando datos de CMB. En Kujat fc Scherrer ( 2000l ) realizan un 
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A a / ttQ 



Figura 3.8: Comparacion de las cotas existentes para la variacion de a entre 
la recombinacion y la epoca actual. La Imea vertical indica variacion nula. 
Referencias: 1). Martins et al. 2002, Rocha et al. 2003 - 2). Ichikawa et al. 
2006 - 3). Stefanescu. 2007 - 4). Nakashima et al. 2008. - 5). Resultado de 
esta tesis (Ec. ^^). 



analisis de como afecta al espectro de potencias de las anisotropfas de CMB 
un valor distinto de la masa del electron, y estiman la sensibilidad de los 
experimentos a la variacion de mp, pero no se llevan a cabo estudios es- 
tadisticos. En IYoo fc Scherreii (120031 ) se realiza un estudio de los h'mites que 
imponen los datos de nucleosmtesis primordial y CMB sobre la variacion del 
valor medio de vacio del campo Higgs (que, como mencionamos en la Intro- 
duccion, durante la recombinacion solo afecta a la masa del electron). Estos 
autores utilizan datos de VSA, BOOMERANG, DASI, MAXIMA y CBI, to- 
dos previos a WMAP. Presentan graficos con los contornos de probabilidad, 
y mencionan que las cotas obtenidas con CMB son mas debiles que las obte- 
nidas con los datos de BBN, pero no publican explicitamente los valores que 
obtienen. 
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En esta tesis se presenta la linica cota que existe en la literatura para 
la variacion de la masa del electron encontrada a partir de datos del fondo 
cosmico de radiacion. En los analisis estadisticos que realizamos para este 
estudio, nuevamente hemos utilizado los datos de CMB, el prior de HST 
y el espectro de potencias del 2dFGRS para poner cotas a los parametros 
cosmologicos y al valor de la masa del electron cuando se permite variar esta 
ultima durante la epoca de recombinacion. A continuacion detallamos los 
resultados. 

Para el analisis utilizando todos los datos de CMB, incluyendo los datos 
de WMAP3, obtuvimos el siguiente intervalo de confianza de 95 % 

- 0,154 < Ame/(me)o < 0,091 (3.24) 

aunque la degeneracion entre rrie y Hq es tan grande que las cotas no resultan 
confiables, ya que las cadenas Uegan a los extremos del prior piano para Hq 
(40 < Hq < 100). Al agregar a este conjunto de datos el prior para Hq del 
HST Key Project obtenido a partir de objetos cercanos, se rompe la fuerte 
degeneracion, y se obtiene un intervalo mas estricto. El intervalo de confianza 
del 95 % es 

- 0,104 < Ame/(me)o < 0,060. (3.25) 

En este caso, las cadenas fueron cortadas cuando el valor del estimador de 
convergencia fue R — \ = 0,029. 

Cuando realizamos el analisis utilizando los mismos datos de CMB, es 
decir incluyendo WMAP3, mas el espectro P{k) del 2dFGRS, obtenemos 
cotas consistentes, aunque mas estrictas: 

- 0,093 < Ame/(me)o < 0,040. (3.26) 

Esto muestra tambien la consistencia del analisis. En las Figs. (13. 9p y (I3.10p 
se muestra como los contornos se hacen mas pequefios cuando se agrega los 
datos del 2dFGRS al analisis que cuando se agrega el prior para Hq. En este 
caso, las cadenas se cortaron cuando se alcanzo el valor R — 1 < 0,0044. Los 
resultados se muestran en la Tabla (13.31) . 

Finalmente, realizamos el ajuste con los datos de WMAP5 y 2dFGRS. 
Las cadenas se cortaron cuando R~l = 0,0169. En este caso, el intervalo de 
confianza del 95 % es 

- 0,071 < Ame/(me)o < 0,053. (3.27) 
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Figura 3.9: Contornos de confianza del 68 % y 95 % para los parametros 
Ame/(me)o y ^bh"^ obtenidos con datos de WMAP3, con datos de WMAP3 
mas el espectro de potencias P{k) del 2dFGRS, y con WMAP3 y el prior 
para Hq calculado a partir de objetos cercanos, para el analisis con variacion 
de rUe. 

En la Fig. (13.1 ip se muestran las distribuciones a posteriori marginali- 
zadas, para el caso de variacion de rUe obtenidas con los datos de CMB, 
incluyendo los datos de WMAP5 y el espectro de potencias del 2dFGRS. 
En la diagonal se muestran las distribuciones a posteriori para los parame- 
tros individuales. Los otros paneles muestran los contornos 2D para pares de 
parametros, marginalizando sobre los demas. En dicha figura se puede ver 
una fuerte degeneracion entre m-e y 0, y tambien entre rrie y Qbh"^, J entre rrie 
y QcDAih"^- En el caso de utilizar los datos de WMAP3, en lugar de WMAP5, 
se observan las mismas degeneraciones, solo que en el ultimo caso, los ajustes 
son mas precisos. 

Los resultados de los ajustes, con WMAP3 y WMAP5, para los parame- 
tros cosmologicos y Ame/(me)o se muestran en la Tabla fl3.3l) . Comparando 
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Figura 3.10: Contornos de confianza del 68% y 95% para los parametros 
Ame/(me)o y Hq obtenidos con datos de WMAP3, con y sin datos del es- 
pectro de potencias del 2dFGRS, para el analisis con variacion de rrie. 

con los valores de los parametros en el caso canonico, se observa que el valor 
de flbh"^ es menor cuando se permite que uie van'e, para ambos conjuntos de 
datos. Por otra parte, para el ajuste con WMAP3, el valor preferido para 
es menor cuando rrie varia que en el caso canonico. 

En la Fig. (13.121) se compara como cambia la distribucion de densidad de 
probabilidad para el parametro Ame/(me)o cuando se consideran los datos 
de WMAP3 y los de WMAP5. En ambos casos, tambien se utilizaron los 
datos del espectro de potencias del 2dFGRS. 

3.3.3. Variacion conjunta de a y rrie 

Para poner cotas a la variacion conjunta de a y me durante la recombi- 
nacion, al conjunto de parametros (13.181) le agregamos las variables Aa/a^ 
y Ame/(me)o. 
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Figura 3.11: Distribuciones a posteriori marginalizadas, para el caso de varia- 
cion de me obtenidas con los datos de CMB, incluyendo los datos de WMAP5 
y el espectro de potencias del 2dFGRS. En la diagonal se muestran las dis- 
tribuciones de densidad de probabilidad a posteriori para los parametros 
individuales indicados en el eje x. Estas distribuciones estan marginalizadas 
sobre los otros parametros y normalizadas a 1 en el maximo. Los otros paneles 
muestran los contornos de confianza para pares de parametros. Los contornos 
representan 68%, 95% y 99% de confiabilidad, desde el mas interno al mas 
externo, respectivamente. Analisis estadistico con variacion de mg. 
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Parametro 


ASriV/rAOO 1 ID/ u\ 

WiVlAro + ("^j 


WMAro ± -f (fcj 




U,Uzi ± U,UUiU 


0,0ZZ4 ± U,000o 




u,iui ± u,uuy 


U,iUo ± U,UU 1 


e 


1,020 ±0,025 


1,034 ±0,022 




u,uyi_Q Q^4 




Ame/(me)o 


-0,029 ± 0,034 


-0,010 ±0,030 




0,960 ±0,015 


0,963 ±0,013 


log(10iOA,) 


3,020 ± 0,064 


3,055 ±0,037 




68,1^^ 


70,6 ± 5,8 



Tabla 3.3: Valores medios y errores a Icr para los parametros incluyendo la 
variacion de rrie, para el analisis con WMAP3 y el analisis con WMAP5. Las 
unidades de Hq son km s~^ Mpc~^ 




0.15 



Figura 3.12: Distribucion a posteriori unidimensional para Amg/ (me)o5 com- 
parando los ajustes con WMAP3 y WMAP5. Analisis con variacion solo de 
me. 
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El linico trabajo en la literatura que presenta ba cotas para l a varia cion 



conjunta de a y me usando datos de CMB es el de llchikawa et al.l (120061 ). En 
dicho trabajo se propone una relacion funcional para las variaciones de a y 
rrie que estaria motivada en la teoria de cuerdas. Los resultados que presentan 
son los limites sobre la variacion temporal de la constante de estructura fina 
a utilizando los datos de WMAPl y el prior gaussiano para Hq dado por el 
HST Key Project. Como no se conoce la relacion entre las variaciones porque 
la teoria esta en construccion, suponen que me varfa con una ley de potencias 
de a. A partir de la accion a bajas energias derivada de la Teoria de cuerdas 
heterotica en el marco de Einstein, derivan la relacion 



me + Ame fa + Aa\^^'^ 



Tambien consideran otras posibilidades, fenomenologicamente, pero con for- 
ma de ley de potencias: 

m^^^/a + Aay ^^^^^^ 
(me)o V " / 

con p = 0, 2, 4. Notar que el caso p = corresponde a variacion solo de a. La 
cota que obtienen en este caso fue comparada en la seccion 13.3. 1[ 

El metodo estadistico que utilizan es el calculo del minimo del como 
funcion de a. Para minimizar sobre los otros parametros cosmologicos, uti- 
lizan iterativamente el metodo de Brent, de la interpolacion sucesiva por 
parabolas. La cota que obtienen es 

- -0,097 < Aa/a< 0,034 al 95 % usando solo los datos de WMAPl, y 

■ —0,042 < Aa/a < 0,026 al combinar estos datos con la cota para Hq 
que se obtiene del HST Key Project. 

En esta tesis se realizo el unico analisis que hay en la literatura sobre la 
variacion conjunta de a y rrie, permitiendo que las r aismas sean in d epend ien- 



tes. Parte de estos resultados estan publicados en iLandau et al.l (120081 ). Se 
utilizaron datos actualizados de CMB, y en distintos estudios se combine esta 
informacion con el prior para Hq calculado a partir de los objetos mas cer- 
canos, y con el espectro de potencias P{k) del catalogo 2dFGRS . Nuestros 



limites no pueden ser comparados con los de llchikawa et al.l (120061 ) dado que 
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estos autores realizan el ajuste fijando una relacion determinada entra las 
variaciones de ambas constantes. Los resultados se detallan a continuacion. 

En el analisis estadi'stico que realizamos con los datos de CMB incluyendo 
WMAP3, encontramos los siguientes intervalos de 95 % de confianza 

- 0,077 < Aa/ao < 0,024 - 0,135 < Ame/(me)o < 0,237 (3.30) 

y al combinar estos datos con el prior para Hq calculado a partir del HST 
con objetos cercanos, los intervalos son 

- 0,078 < Aa/ao < 0,023 - 0,096 < Ame/(me)o < 0,215. (3.31) 

Vemos que no se modifican demasiado los intervalos de confianza al agregar 
el prior para Hq. 

En el analisis con los datos de WMAP3 mas el P{k) del catalogo 2dFGRS, 
cortamos las cadenas cuando R — 1 < 0,015. Los intervalos del 95 % de 
confianza en este caso son 

- 0,068 < Aa/ao < 0,031 - 0,116 < Ame/(me)o < 0,198. (3.32) 

Los resultados se muestran en la Tabla (13.41) . Los valores de los parametros 
cosmologicos son similares a a quellos obtenidos cons i derando la variacion d e 
una sola constante a la vez ( Mosquera et al. 2008; Scoccola et al. 2008bl ). 



Ademas, los resultados son consistentes dentro de la con variacion nula de 
a y me durante la recombinacion. 

En las Fig. (I3.13p y (13.141) comparamos las degeneraciones que existen en- 
tre diferentes parametros cosmologicos y las constantes fundamentals cuan- 
do una o ambas constantes pueden variar. En ambos casos, la region permi- 
tida en el espacio de parametros es mayor cuando las dos constantes pueden 
variar. Esto es esperable dado que en ese caso, la dimension del espacio de 
parametros es mayor y por lo tanto, las incertezas en los parametros son ma- 
yores. Las correlaciones de ^ con los otros parametros cosmologicos cambia 
el signo cuando se permite variar tambien mg. Por el contrario, las correla- 
ciones de con los parametros cosmologicos no cambian el signo cuando 
se permite tambien variar a. 

Cuando solo se permite variar una constante fundamental, la correlacion 
entre esa constante y cualquier parametro cosmologico particular tiene el 
mismo signo, sin importar si la constante es a o rrie- Esto puede ser entendido 
sobre la base de que ambas constantes entran en las mismas cantidades fisicas. 
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Figura 3.13: Contornos de confianza a Icr y 2a. Lmea azul: variacion de a y 
rrif,] lmea roja: Solo variacion de a. Los parametros cosmologicos son libres 
de variar en ambos casos. 




0.94 0.96 0.98 1 1.02 1.04 1.06 1.08 1.1 1.12 



Figura 3.14: Contornos de confianza a lo" y 2cr. Lmea azul: variacion de a y 
rrie] lmea roja: Solo variacion de mg. Los parametros cosmologicos son libres 
de variar en ambos casos. 

Sin embargo, dado que la forma funcional de la dependencia es diferente, 
los valores medios del mejor ajuste para las variaciones temporales de las 
constantes fundamentals son diferentes y la distribucion de probabilidad es 
mas extendida en un caso que en el otro. De todas maneras, en los casos en 
que se permite variar solo una constante, esta prefiere un valor mas bajo que 
el valor actual. 

Finalmente, realizamos el ajuste con los datos de WMAP5 y 2dFGRS. 
Las cadenas se cortaron cuando R — 1 = 0,018. Los resultados se muestran 
en la Tabla (13.41) y en la Fig. fl3.15p . En este caso, los intervalos de confianza 
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Parainetro 


winapS 


wniap3 




0,0224 ± 0,0009 


0,0218 ± 0,0010 




0,107 ±0,008 


0,106 ±0,011 


e 




i,UOt)_Q Q29 


T 


0,086 ± 0,008 


0,090 ±0,014 


Aa/aQ 


0,003 ±0,015 


-0,023 ±0,025 




-0,017 ±0,051 


0,036 ± 0,078 


ris 


0,963 ±0,015 


0,970 ±0,019 


As 


3,052 ± 0,043 


3,054 ± 0,073 


Ho 


70 3+^'^ 

'U, -3-6,0 


70 4+''''' 



Tabla 3.4: Valores medios y errores a la para los parametros, incluyendo va- 
riacion de a y rrie. Comparacion de resultados entre los analisis con WMAP3 
y WMAP5. 



del 95 % son 



- 0,026 < Aa/ao < 0,036 - 0,115 < AmJ{me)o < 0,086. (3.33) 



(nie)o 



Son notables las fuertes degeneraciones que existen entre 

° ^ {me)o 



COM 



y G, — y n^, y finalmente entre 



Aa 



y 



Ame 

(me)o 



. Los valores obtenidos para 



h, VtcDMh . Tr-p, y ns_ concuerdan, a la, con aquellos del equipo de 



WMAP ( jSpergel et al.ll2007l ). donde no se considera variacion de a ni nie- 
Los resultados para los valores de los parametros cosmologicos son similares 
a aquellos obtenidos considerando la variacion de una sola constante a la vez. 
Los resultados son consistentes dentro de la con variacion nula de a y me 
durante la recombinacion. 

Al comparar los contornos de las distribuciones marginalizadas para Aa / 
versus Ame/(me)o, que se obtienen de los analisis con los datos de WMAP3 
y WMAP5, se distingue que los contornos no solo se hacen mas pequefios 
para los nuevos datos, sino que tambien cambian los valores medios de los 
parametros (ver Fig. ( 13.16^ ). Los valores medios de Aa/ao y de Ame/(me)o 
cambian de signo, como puede verse en la Tabla (13. 4p . 
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Figura 3.15: Distribuciones a posteriori marginalizadas, para el caso de va- 
riacion de a y mg obtenidas con los datos de CMB, incluyendo los datos de 
WMAP5 y el espectro de potencias del 2dFGRS. En la diagonal se muestran 
las distribuciones de densidad de probabilidad a posteriori para los parame- 
tros individuales indicados en el eje x. Estas distribuciones estan marginali- 
zadas sobre los otros parametros y normalizadas a 1 en el maximo. Los otros 
paneles muestran los contornos de confianza para pares de parametros. Los 
contornos representan 68 %, 95 % y 99 % de confiabilidad, desde el mas in- 
terno al mas externo, respectivamente. Analisis estadi'stico con con variacion 
conjunta de a y mg. 
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Figura 3.16: Contornos de confianza del 95%, para los parametros Aa/ao 
y Ame/(me)o, correspondientes a los ajustes realizados con WMAP3 y con 
WMAP5. 

3.4. Comparacion entre escenarios de recom- 
binacion 

Dado que los datos del fondo cosmico de radiacion son cada vez mas preci- 
sos, es necesario tener en cuenta fenomenos cada vez mas sutiles que puedan 
afectar el escenario de recombinacion. En el aiio 2008, la recombinacion del 
helio fue estudiada con sumo cuidado, como se explico en el Cap. [21 En esta 
tesis hicimos un analisis detallado de como entran las constantes fundamen- 
tales en la recombinacion del helio, y modificamos la version actualizada del 
codigo de evolucion Recfast para introducir dichas dependencias. En esta 
seccion mostramos cuales son las cotas para la variacion de a y rrie cuando 
se tiene en cuenta la recombinacion detallada del helio, y comparamos con 
las cotas obtenidas en la seccion anterior, en la cual se considero el esce- 
nario estandar de la recombinacion, descripto en ISeager et al. fl2000f ). Los 
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Analisis 




Ame/(me)o 


Variacion conjunta - PS 


0,003 ±0,015 


-0,017 ±0,051 


Variacion conjunta - NS 


0,004 ±0,015 


-0,019 ±0,049 


99999Variaci6n a - PS 


-0,002 ± 0,009 




Variacion a - NS 


-0,001 ±0,009 




Variacion mg - PS 




-0,010 ±0,030 


Variacion nif, - NS 




-0 010+°'°^^ 

U,U±U_Q Q3Q 



Tabla 3.5: Cotas para la variacion de las constantes en distintos escenarios 
de recombinacion. NS significa el nuevo escenario de recombinacion, y PS 
significa el escenario estandar. Se presentan los resultados para todos los 
analisis. 



resultados de este analisis estan publicados en IScoccola et al.l (l2008al ). 

Realizamos el analisis estadistico teniendo en cuenta la posible variacion 
de a, de me, y de a y me simultaneamente, en la epoca de la recombinacion 
en el escenario descripto en la Sec. 12. 1[ Para este analisis utilizamos los datos 
de CMB incluyendo WMAP5, y el espectro de potencias P{k) del 2dFGRS. 

Exploramos el espacio de parar aetros con 8 cadenas de Markov, y segui- 
mos el criterio de convergencia de iRaftery &: Lewis fll992f ) para detener las 
cadenas cuando R — l< 0,0180. Las cotas sobre los parametros no se ven ma- 
yormente afectadas cuando se cambia el escenario de recombinacion, como 
puede verse en la Tabla (13.51) . Los valores de los parametros cosmologicos 
practicamente no cambian entre este ajuste y el realizado teniendo en cuenta 
el escenario estandar. 

En la Fig. (13.171) . comparamos la distribucion de probabilidad para Aa/ao 
en diferentes escenarios. En la Fig. (I3.18p . hacemos lo mismo para Am-e/ (r77.e)o- 
Puede verse que la diferencia es despreciable. 

En la Fig. (I3.19P mostramos como una variacion en el valor de a en la 
epoca de recombinacion afecta la historia de ionizacion, moviendo el redsliift 
al cual ocurre la recombinacion a tiempos anteriores para mayores valores de 
a. La diferencia entre las funciones cuando se consideran los dos escenarios 
de recombinacion, para un dado valor de a, es menor que la diferencia que 
existe cuando se cambia el valor de a. Algo similar ocurre cuando se varia 
rrie. 
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-0.06 -0.04 -0.02 0.02 0.04 0.06 0.08 



A a / Uq 



Aa 

. uiiiumieiibiuiiai para 

combinacion diferentes. 



Figura 3.17: Likelihood unidimensional para para dos escenarios de re- 



3.5. Relacion fenomenologica 

Como mencionamos en la Introduccion, el ajuste de la variacion conjunta 
de a y mg se ha realizado en esta tesis de manera totalmente independiente 
de cualquier modelo que relacione dichas variaciones. Por lo tanto, la correla- 
cion encontrada entre las variaciones de estas constantes nos permite estimar 
una relacion fenomenologica entre ambas. A partir del contorno de confia- 
bihdad para Aa/ao y Ame/(me)o del ajuste mas actualizado (es decir, el 
conjunto de datos que incluye WMAP5), se puede estimar una relacion entre 
las variaciones de a y mg entre la epoca de recombinacion y la actualidad. 

Una manera de hacerlo es proponer una relacion lineal del tipo 

Arrie Aa 

acMB VocMB (3.34) 



(me)o tto 
y buscar valores apropiados para acMB Y be me 
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-0.2 -0.15 -0.1 -0.05 0.05 0.1 0.15 0.2 

A mg / (me)o 



Figura 3.18: Likelihood unidimensional para 



(me)o 



para dos escenarios de 



recombinacion diferentes. 



Para encontrar estos coeficientes podemos hacer una regresion ortogonal, 
que consiste en ajustar una elipse al contorno 2D de la Likelihood marginali- 
zada sobre todos los parametros excepto Aa/ao Y Ame/(me)o, y la relacion 
entre ambas variaciones estara dada por el semieje mayor de la recta. 

La ecuacion de una elipse rotada un angulo respecto de los ejes {x', y') 
(cuyo origen se encuentra en el punto (01,02) del sistema {x,y)) es 



Ci (x - aif + C2 (y - 02)^ + C3 (x - ai){y - 02) = 1 



(3.35) 
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Figura 3.19: Historia de ionizacion permitiendo que varie a con el tiempo. 
De izquierda a derecha, los valores de aja^ son 1.05, 1.00, y 0.95, respectiva- 
mente. Las Imeas azules (rayas) corresponden al escenario de recombinacion 
estandar, y las Imeas rojas (solidas), corresponden al escenario actualizado. 



donde 



(1 + ^2)^2 (1 



1 + ^2)^2 (1 + ^2)^2 

y 9 = tan(0). Ajustando el contorno de confiabilidad a una elipse, encontra- 
mos valores para los parametros de la misma, es decir, para a, 6, 6*, ai, y 02. 
Luego, dado que la recta que coincide con el semieje mayor de la elipse tiene 
la ecuacion 

y = 6 {x — ai) + a2 = X + {a2 — 6ai) (3.37) 
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podemos encontrar los valores acMB Y bcMB Y sus errores, los cuales calcu- 
lamos haciendo 



(^acMB = (3.38) 
^bcMB = (^a2 - o-ic^e - dc^ai (3.39) 

Encontramos que acMB = "4,2^0^45 y Bcmb = ~0,007^q^q3. Es decir, la 
relacion entre las variaciones de a y rrie es 

Arrip Aa 

- — f = -4,2 0,007. (3.40) 

(me)o ao 

Otra manera de encontrar una relacion entre Aa/aQ y Ame/(me)o es a 
traves de la regresion mmimo cuadratica. En este metodo, dada la distribu- 
cion conjunta de dos variables aleatorias X e Y, que puede expresarse como 
una distribucion de "masa" en el piano {X, Y), se determina una recta en el 
piano con la condicion de que dicha recta proporci one el nieior a juste posible 



a la masa de la distribucion. Se puede demostrar (jCramerlll972l ) que la recta 
tiene la expresion 

^ = (3.41) 

es decir 

y = 02 H [x — ai) = X + 02 Oi . (3.42) 

cTi (Ji V / 

Para encontrar los valores de p, ai, y cr2, consideramos que para una 
distribucion normal de 2 variables, la matriz de covarianza es C, tal que sus 
elementos son 



cii = al Ci2 = cov(xi,a;2) 

C21 = COv(xi,X2) C22 = CTg (3.43) 

Sea p = cov{xi, X2) / {cria2) ■ Cuando se calculan las Imeas de probabilidad 
constante (es decir, los contornos de probabilidad), se Uega a la ecuacion de 
una elipse (centrada en oi y 02): 

^^l^-2pi^±^^^ + i^^^ = l-p^ (3.44) 
a( ai (T2 

o lo que es lo mismo: 
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(xi - ai] 



- 2- 



Xi - ai X2 - a2 (x2 - 02)' 



0"T 



de donde vemos que 



0-1 



02 



(3.45) 



C3 



cr, 



Por lo tanto. 



(3.46) 



Luego, con el valor de Ci obtenemos ai, y con el valor de C2 obtenemos 
(T2. Con estos valores y los valores medios ai y 02 del ajuste anterior, de la 
Ec. f l3.42p encontramos la recta de regresion 



- — f = -2,57 0,01. 



(3.48) 



Esta relacion es diferente de la relacion encontrada en Ec. fl3.40p . Sin em- 
bargo, el signo de la pendiente es negativo en ambos casos, debido a que 
^rripj {mf.)n y Aa/an estan anticorrelacionadas. La relacion propuesta por 
Ichikawa et al.l (120061 ) . y con la cual hicieron los ajustes es 



Amg 

(?^e)o 



P- 



Aa 

ao 



(3.49) 



con p fijo. Encuentran buenos ajustes para p = 0, 1/2, 2, 4. Todos estos valores 
son positivos, y aqui se ha mostrado que, para el mejor ajuste de los datos 
(cuando se deja Aa/ao y Ame/(me)o como parametros libres), la correlacion 
entre las variaciones es negativ a, es decir, p deberia se r negativo. Por lo tanto, 
concluimos que la propuesta de Ichikawa et al.l (2006) deberia ser descartada. 
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Capitulo 4 

Modelo de Barrow &: Magueijo 



Como hemos discutimos en la Introduccion, la variacion de las constantes 
fundament ales es una prediccion de teon'as que intentan unificar las fuerzas 
de la naturaleza. En el limite de bajas energias, estas pueden ser interpretadas 
como teorias efectivas. Por otra parte, existen varios modelos fenomenologicos 
que, a partir de unas pocas hipotesis generales, describen la variacion de 
ciertas constantes fundamentals. A su vez, estos modelos pueden ser Ifmites a 
bajas energias de teorias rr ias generales. Entre ellos se encuentra el modelo de 



Barrow fc Magueijo! (120051 ) . que describe la evolucion de la masa del electron. 



Si bien, en algunos modelos, la variacion de la masa del electron aparece 
vinculada a la variacion de la constante de estructura fina, no es necesario que 
sus variaciones esten relacionadas, y diclios autores consideraron oportuno 
estudiar las consecuencias cosmologicas de un modelo simple de variacion de 
la masa del electron, independientemente del comportamiento de a. 

Aqui deduciremos soluciones analfticas aproximadas mas precisas que las 
de Barrow & Magueijo (2005). En la Seccion [4.11 resolvemos en forma mas 
rigurosa la ecuacion de Friedmann para el factor de escala, y luego mante- 
nemos todas las constantes de integracion en la resolucion de la ecuacion 
diferencial para el campo escalar. Se obtienen asi soluciones detalladas para 
la variacion de la masa del electron como funcion del tiempo, que se mues- 
tran en la Seccion IT2l Finalmente, los parametros del modelo seran acotados 
utilizando datos. 

Las diferentes observaciones astronomicas y geoffsicas que se pueden uti- 
lizar para poner cotas a la variacion de la masa del electron, corresponden a 
distintas epocas en la historia del Universo. Esto las hace realmente utiles pa- 
ra testear modelos que predigan un comportamiento temporal determinado. 
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En esta tesis se consideran distintas observaciones, las cuales se describen en 
la Seccion I4.3[ y se encuentran cotas a los parametros del modelo de Barrow 
& Magueijo utilizando estas observaciones a distintos tiempos, en la Seccion 

m 

Luego, en la Seccion S3] obtenemos limites adicionales para estos parame- 
tros, a partir de las cotas existentes para la validez del Principio de Equiva- 
lencia Debil, que se obtienen de experimentos tipo Eotvos. Comparando con 
las cotas obtenidas a partir de resultados cosmologicos, se Uega a la conclu- 
sion de que el modelo debe se r descartado. El tratam iento desarroUado en 
este Capitulo fue publicado en IScoccola et al.l (l2008bl ) . 



4.1. El formalismo 

La densidad lagrangiana del campo libre de Dirac, a part ir de la cual s e 



puede derivar la ecuacion de movimiento para el electron, es ( lSakurailll967l ): 



=Sf = —chip'yfj,{d/dXfj,)ip — meC^ip^. (4.1) 

Consideremos que la masa del electron esta controlada por un campo 
escalar "dilatonico" definido por rrie = meoexp[0], donde es el valor 
actual de la masa, que corresponde a = hoy. Para variaciones pequenas 
de nie, desarroUamos exp[0] = 1 + 0. La densidad lagrangiana total es 



^ = ~ch^-f^{d/dx^)tlj - rrieoc^^^ - meoC^^^p(j) + -d^<l)d^'(j). (4.2) 

Por lo tanto, la ecuacion de movimiento del campo se obtiene de la siguiente 
densidad lagrangiana 

=Sf = ^int + -^(j) (4-3) 

donde el termino cinetico es 

= ^d,<f>dy (4.4) 

y el termino de interaccion entre el campo y el campo de Dirac del electron 



es 



■^int = -rrieoc ipipcf). (4.5) 



62 



Este modelo tiende al modelo canonico sin variacion de la masa del electron 
cuando el valor del parametro c<j — oo. 

La ecuacion de Euler- Lagrange a resolver es la siguiente: 

d / \ 



dx^, \d{d(t)ldx^) J d(f) 

Escribimos explicitamente JSf^, teniendo en cuenta que la convencion de 
signos para la metrica es (+ — ), y que (f) es funcion solo del tiempo: 



De aqui, 

por lo tanto 
Por otra parte 



90 

Por lo tanto, la ecuacion de Euler- Lagrange queda: 



do 



(4.7) 



= (4-8) 



d {do(j)) dx^ c dt 
a=Sf \ _ d^ _u d'^(t) 



By 

= -meocVV'- (4.10) 



dXfj, \d (d(f)/dXfj,) J d(f) \9(9o0)/ d(j) 

= ^^+^eOc'^^ = (4.11) 

que se puede reescribir como: 

— = c V'V'- 4.12) 

dt^ uj 

Reemplazamos ■0V' por su promedio macroscopico, que en el Kmite no 
relativista es igual a + n^, siendo rig el numero de electrones y Up el 
niimero de positrones. Despreciando Up frente a rZe, obtenemos 
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,,2 - cW (4.13) 

Pasando a derivadas covariantes, ^ + 3^0, por lo tanto, la ecuacion 
queda 



o equivalentemente 



+ 3-0 = a-' ( <Pa^] = -^cV (4.14) 



^^3^ ^ = -^cVo- (4.15) 

/ UJ UJ 



donde hemos tenido en cuenta que a^rie = alueo y Oq = 1. 
Integrando la Ec. ( 14. 15^ obtenemos: 

<i)a^ = -Mt + A (4.16) 

con 

M = 1}^c' = P^c' (4.17) 

UJ UJ 

y A una constante de integracion. Esta es la ecuacion que debemos resol- 
ver para el campo 0, teniendo en cuenta que el factor de expansion a es 
una funcion del tiempo que debemos encontrar para el modelo cosmologico 
considerado. 



4.1.1. Contribucion de a la Ecuacion de Friedmann 

En esta tesis se deriva una ecuacion mas precisa para el campo de este 
modelo. Para ello, comencemos considerando la ecuacion de Friedmann 

/ay^c^_87rG^ (4 18) 

\aj a? 3c^ 

donde e = e^/ + + eA es la densidad de energfa total y a = da/dt. Es- 
cribimos las densidades de materia, radiacion y energfa oscura (o constante 
cosmologica) , respectivamente, en terminos de sus valores actuales: 

CAf = (^) , e/J = e/Jo (^) , = ^Ao (^) • (4.19) 
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Definimos Q = p/ pcrit donde la densidad critica es: 

Pcrit = ^yr (4.20) 

J H = a/a es la tasa de expansion de Hubble al tiempo t. 

Sea un universe espacialmente piano {k = 0). Multiplicando Ec. fl4.18p 
per Hq^, obtenemos: 



SttG 



Esta ecuacion puede escribirse como 



/ao\3 /ao\4 



(4.21) 



1 / da\'^ ^Imo , ^iJo 



Kdr / a^ 



+ + ^Ao (4.22) 



donde hemos realizado el cambio de variables: r = HqI y considerado = 1. 
Notar que r es adimensional. 

El campo <p tiene una densidad de energfa dada por 



= ^ujd,cl>d^(t) = -u (^00)^ = ^-<P\ (4.23) 
Usando la ecuacion fl4.16p . la podemos escribir como: 

Esta densidad de energfa contribuye como un termino mas en la ecuacion 
de Friedmann, que queda 



2 / M , A 

a2 Ut) ~ a3 + a4 +^^^o+ 3 a^ 



2 



Para tiempos muy remotos, los terminos mas importantes en la Ec. fl4.25p 
son los correspondientes a la radiacion y al campo 0. Para poder despreciar los 
efectos del campo en la ecuacion de Friedmann, pedimos que la contribucion 
debida a este sea mucho menor que la contribucion de la radiacion, es decir: 



( _ M ^ I A 
4:71 GUO \ HfJ'^l ^ Ho 



2 

- < ^Ro (4.26) 
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donde Tei es el tiempo r a partir del cual vale este modelo, cuando la tempe- 
ratura del universo corresponde a la energi'a de produccion de pares electron- 
positron. 

En lo que sigue, damos una estimacion del valor de Tei y del factor de 
expansion correspondiente, aei- Ambos aparecen en la condicion que impon- 
dremos al hacer el ajuste de los parametros A j M del modelo. 



Regimen de validez del modelo 

El modelo de Barrow & Magueijo tiene validez a partir del momento en 
que la temperatura del Universo tiene un valor correspondiente a la energia 
minima requerida para la formacion de pares electron-positron. En esta epo- 
ca, radiacion y materia estan en equilibrio, formando un solo fluido, de ma- 
nera tal que 

E^ = kT = 2m^(? (4.28) 

donde k es la constante de Boltzmann. Entonces, la temperatura a la cual se 
cruza el umbral de creacion de pares es 

Tel = 4.29 
k 

Por otro lado, sabemos que la temperatura decae con la expansion de manera 
inversamente proporcional 

aT = cte^T,i = ^. (4.30) 

Podemos estimar t^i a partir de la ecuacion de Friedmann para la epoca de 
dominio de la radiacion: 



a2 



(^)' = ^ ^ ada = V^dr \ = V^^l^f • (4-31) 



^Lo que habria que pedir es: 



A/ _ I _A_ 

Ho 



3 c4 a(r)2 



2 

- « f^flo (4.27) 



para todo t > t^i. 

Sin embargo, el miembro izquierdo de la desigualdad es varios ordenes de magnitud 
menor que el miembro izquierdo de la Ec. (l4.26p para los datos con los que contamos en el 
regimen de dominio de materia y radiacion, es decir BBN y CMB. 
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Suponemos que Ti = y Oi = 0. Entonces 



aei = V2Q'^Tf ^ ^ = ^ = ^ (4.32) 



Te. = //ote. = J^^;,yM ) • (4-33) 



donde hemos considerado Oq = 1. De esta mancra 

Para obtener el valor numerico de este tiempo, tomamos los siguientes valores: 

kTQ = 2,35 X 10~^ eV (temperatura de CMB) 

2me(? — 1,02 X 10^ eV (energia para la produccion de pares e~ — e^) 

Qrq/i^ = 2,47 X 10-^ 

h = //o/lOO km Mpc"^ = 0,73 

H^^ = 4,22x10^^8 (4.34) 

Esto da un valor para el tiempo cuando cesa la produccion de pares e~ — e"*" 
igual a 

iei = 1,52 X 10-^s. (4.35) 

El valor numerico de a^i es 

-10 



Uei = 2,304 X 10"'". (4.36) 



4.2. Soluciones 



Para encontrar la solucion de la ecuacion del campo primero es nece- 
sario resolver la ecuacion diferencial para el factor de expansion a{t). 

4.2.1. Solucion de la Ecuacion de Friedmann 

Buscamos la expresion para a{t) resolviendo la ecuacion de Friedmann. 
Estudiamos dos etapas del universo: 

1. regimen de dominio dc materia y radiacion 

2. regimen de dominio de materia y constante cosmologica. 
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La transicion ocurre cuando la contribucion a la ecuacion de Priedmann, 
tanto de la radiacion como del tcrmino de constante cosmologica, son iguales. 
Esto sucede cuando el factor de expansion toma el valor: 




(4.37) 



Primer regimen 



Consideremos un universo dominado por materia y radiacion. La ecuacion 
para el factor de expansion queda: 

= + ^ (4.38) 

Haciendo el cambio de variables drj — dr/a, integrando, y tomando como 
condiciones iniciales que at — 0, a — Oyrj — 0, Uegamos a que 

2 



7] 



Despejamos a{r)): 



^J^M.a + fiijo - • (4-39) 



a\r^) = ^rj- + ^rj (4.40) 

donde el supramdice R denota que estamos en la epoca dominada por materia 
y radiacion. Para encontrar r en terminos de rj, seguimos los siguientes pasos: 

dr = adri^T= f dr' = f a{ri')dr,' = ^rj^ + ^^^2 ^^^^^ 
Jo Jo 12 2 

es decir que la relacion entre r y 77 es: 

r^'^rj^ + ^n^ (4.42) 

Segundo regimen 

Consideremos ahora un universo dominado por materia y constante cos- 
mologica. La ecuacion para el factor de expansion queda: 
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Cambiando a la variable y — \ATa^7^^JWo^^^^) Uegamos a esta ecuacion: 
2 1 dy 



dT 



2 1 . n -1 
T — - — - — smh y + C . 



3 ^V^. 
Reemplazando la expresion de y y despejando a obtenemos: 

"3 



(4.44) 



aHr) 



Ao 



-^/fUoT + C' 



(4.45) 



Las constantes de integracion se eligen de manera tal de asegurar la con- 
tinuidad de a{ri) y de su derivada en el punto de ajuste tja, donde ambas 
soluciones se superponen. Por lo tanto, el punto de ajuste es aquel donde se 
cumple: 



^0 



Va + V^oVa = \i 



Ro 



Ao 



(4.46) 



Resolviendo la ecuacion cuadratica, y teniendo en cuenta que rj es siempre 
positivo, obtenemos el tja' 



Ro 



Mo 



1 + 



^Mq 4/ ^fio 
^Ro V ^Ao 



- 1 



Si Uamamos 



podemos escribir 



. _ ^Mo _ o Q-3/4 (^-1/4 



^Ro V ^Ao 



Va 



Ro 



n 



Mo 



(4.47) 



(4.48) 



(4.49) 



En la variable r, el tiempo de ajuste es: 



ta 



MoH 



Va + 



Ro 2 

Va 



2 fl 



Ta 



3/2 

Ro 



12 2 

Considerando 0^(77^) = a^(rA), ajustamos la constante C: 



2-(2-e)yrTel- (4.50) 



(4.51) 
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Resumiendo: 

En ta, se pega con de manera continua, con derivada primera 
continua. Las expresiones para el factor de expansion en cada uno de los 
regimenes son: 



^v'^ + V^^V donde r = ^r]^ 



Mo 



Ro 



sm. 



h2/3 



(4.52) 



Podemos calcular cuanto vale r hoy, tq, sabiendo que a^(ro) = 1. Asi Ue- 
jamos a que: 



2 1 

To = -- 



3^ 



Ao 



sinli ^ 



Ao 



C 



Mo 



(4.53) 



4.2.2. Soluciones para el campo cj) 
Dominio de materia y constante cosmologica 

La ecuacion para el campo es: 



/,a^ = -Mt + A. 



(4.54) 



Pasando a la variable r, y considerando a = a^(r), Uegamos a esta ecua- 



cion 



4>{ro) 
ir) 



Ao 



, Hi Qmo sinh' [I + C'] 



dr. 



(4.55) 



Integrando (ver Apendice [B]) , y teniendo en cuenta que rUe = meo exp [0] 
y por lo tanto 0o = 0, y ademas 

In ( In f + - U In ^ + ^ V ^ (4.56) 



llegamos a: 
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Amp 



M 2 I /3 \2//3 

Hi [v^^coth [C' + -y^,rj - - In (^sinh j^C' + -v^^ 

^^"^(^7' + In (7+ /IT? 



+ 3 I ln7 - 



+ 



donde hemos Uamado 



7 

Ho 3^4/0 



-coth(c" + -v^r ) + 



1 + 7' 



7 



Ao 



n 



Mo 



7 



(4.57) 



Dominio de materia y radiacion 

Partimos nuevamente de la ecuacion para el campo 

d)a^ = -Mt + A 



y tenemos en cuenta que: 



t = H,\ = H.'^r,^ 



n 



Mo 



Ro 



(4.58) 



(4.59) 



Pasando a la variable rj, y usando a = a^{ri), tendremos que integrar 



^ — dr] . 



(4.60) 



La expresion para (pA esta dada en la Ec. ( ]B.12p . La solucion de esta ecuacion 
es (ver Apendice ICl) : 
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Arrie 



iv) 



2 M 1 

3 ^Imo 



-2 In 



+ 



Ao 



+ 



A Q 



Mo 



2 



- (sinh ""^ 7 — C" 



7 



2 



4(/3r/ + l) 



VTT^ + 1 y 4/3r/ 
(e-i)v^. 3 



2 vi+Z 



1 1/2 

donde por conveniencia hemos definido (3 = -^^Imo^rJ ■ 

Finalmente, damos la expresion para la tasa de variacion de la masa del 
electron, que luego sera util para comparar con las cotas observacionales 
obtenidas a partir de relojes atomicos, en el laboratorio. Dado que 



In 



rrie 



rrie 



Por lo tanto, considerando la Ec. fl4.54p . expresamos 

rhe -Mt + A 



■nip 



a-' 



(4.61) 



(4.62) 



4.3. Datos 

En esta seccion resumimos los h'mites existentes para la variacion de la 
masa del electron, obtenidos a partir de datos observacionales correspondien- 
tes a distintas etapas cosmologicas. 



4.3.1. Nucleosmtesis primordial 

La nucleosfntesis primordial es una de las herramientas mas importantes 
para estudiar el universo temprano. El cociente barion-foton, o equivalente- 
mente la asimetria en los bariones rjB = {ns — n§)/n^ se puede determinar 
tanto a partir de calculos teoricos como de observaciones de las abundan- 
cias de los elementos quimicos livianos. Un cambio en el valor de rrie en la 
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epoca de la formacion de estos elementos respecto de su valor actual afecta 
cantidades fi'sicas tales como la suma de densidades de energi'a de electrones 
y positrones, la suma de sus presiones, y la diferencia entre las densidades 
en numero de electrones y positrones. Los cambios mas importantes en las 
abundancias primordiales debido a variacion de rrie result an del cambio en las 
tasas de decaimiento d ebil, mas que del cambio en la tasa de expansion. En 
Scoccola et al. fl2008bl ) se consideraron los datos disponibles para las abun- 
dancias de D, ^He y ^Li. El analisis estadistico indico que no hay buen ajuste 
para el conjunto total de datos. Sin embargo, se logran ajustes razonables 
si se excluye alguno de los datos. Se discuten las cotas obtenidas utilizando 
todas las abundancias y solo un subgrupo, llegando a un ajuste de 77^ con- 
sistente con el valor de WMAP solo si las abundancias consideradas son las 
de D y '^He. La cota que se obtiene en este caso es la siguiente: 



Arrif. 
(me)o 



-0,022 ±0,009. 



(4.63) 



En el caso de considerar otros pares de abundancias, el valor del mejor 
ajuste para 77^ se aleja muclio del valor obtenido por WMAP. La abundan- 
cia del ''Li no se tuvo en cuenta porq ue aun no hay conse nso acerca de las 
cotas para la misma. Se ha senalado ( jRichard et al.ll2005l ) que debe enten- 
derse mejor el transporte turbulento en las zonas radiativas de las estrellas 
para tener una estimacion c onfiable de la abundancia de ''Li. Por otra parte, 
Melendez fc Ramirea (120041 ) reanalizaron los datos disponibles para el ^Li y 
obtuvieron estimaciones ma rginalmente consistentes con WMAP. Mas aun, 
Prodanovic fc Fields! (120071 ) destacaron que la discrepancia se acentua si se 
considera la contaminacion con ^Li. Es por eso que se adopto el criterio mas 
conservador de considerar la cota que se obtiene al ajustar las abundancias 
de D y ^He. 



4.3.2. Fondo cosmico de radiacion 

El fondo cosmico de radiacion provee una de las cotas para la variacion de 
la masa del electron en el Universo temprano, durante la epoca de formacion 
del hidrogeno neutro. En el Capftulo [3] de esta Tesis, se encontraron varias 
cotas para la variacion de la masa del electron. Para este analisis utilizaremos 
la cota que se obtiene del ajuste con los datos de WMAP3 y el espectro de 
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potencias P{k) del catalogo 2dFGRS: 

- — f = -0,029 ± 0,034. (4.64) 

(me)o 

La diferencia entre esta cota y la obtenida con el conjunto de datos WMAP5 
es mmima a los efectos del ajuste de los parametros del modelo de Barrow 
& Magueijo, y no altera las conclusiones de este Capitulo. 



4.3.3. Sistemas de absorcion en quasares 



Los sistemas de absorcion en espectros de quasares constituyen "labo- 
ratorios" naturales donde buscar variaciones temporales de las constantes 
fundamentales. En particular, un metodo para poner Ifrn i tes a la variacion 
de n = ^ fue desarrollado por lVarshalovich fc Levshakovl (119931 ). Esta basa- 
do en el hecho de que las longitudes de onda de las Ifneas vibro-rotacionales 
electronicas dependen de la masa reducida de las moleculas, con distintas de- 
pendencias para distintas transiciones. De esta manera, es posible distinguir 
el corrimiento al rojo (redshift) cosmologico de una Ifnea, del corrimiento 
provocado por una variacion en fi. Para una dada transicion del hidrogeno 
molecular H2, la longitud de onda en el laboratorio (reposo), A°, se puede rela- 
cionar con la del sistema de absorcion, Aj, mediante = {1 + Zabs){^ + Ki—), 
donde Zabs es el redshift de absorcion y Ki es el coeficiente que determina 
la sensibilidad de la longitud de onda Aj a la masa del electron. Usando 
observaciones de sistemas de absorcion en quasares a muy alto redshift y 
me didas de laboratorio , varios autores han puesto cotas a la varia c ion d e 



f lPotekhin et al.l Il998l : iLevshakov et al.l l2002l : llvanchik et all l2003l . 120051 ) 



Los datos disponibles actualizados hasta el ano 2007, para cada redshift, se 
citan en la tabla 14.11 y son los que consideramos para testear el modelo de 
Barrow-Magueijo. Llamamos a este conjunto de datos QSO I. 

En el ano 2008 se actualizaron las cotas para estos sistemas de absorcion 
en quasares. Esto se logro gracias a una mejora substancial en el metodo 
de calibracion de las longitudes de onda, y tambien incorporando nuevas 
transiciones al ajuste, que antes no se tenfan en cuenta por estar situadas 
dentro de una Imea Lyman a del quasar. Los nuevos datos del conjunto QSO 
I se muestran en la Tabla 14. 2[ En la seccion 14.41 se repite el ajuste de los 
parametros del modelo usando estos datos. 



Un metodo alternativo para medir fi a alto redshift, sugerido por lFlambaum fc Kozlov 



(I2OO7I ). es usar la sensibilidad a las variaciones en /i de las transiciones in- 
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Tabla 4.1: QSO I: Redshift de absorcion, valor de con el correspondiente 
error (en unidades de 10~^), y la referenda. 



Redshift 


(me)n 


Referenda 


3.025 


-1,47 ±0,83 


Ivanchik et al. f2005) 


2.811 


-6,25 ± 13,75 


Potckhin et al. (1998) 


2.5947 


-2,11 ±1,39 


Ivanchik et al. f2005) 



Tabla 4.2: QSO I: Valores actualizados de la Tabla |4?T1 La ultima Imea corres- 
ponde al metodo de transiciones inversas del NH3 (y^y- Y errores en unida- 
des de 10"*^). 



Redshift 


(me)n 


Referenda 


3.025 


8,2 ± 7,4 


Kine et al. (2008) 


2.811 


-1,4 ±3,9 


I£ing et al. (2008) 


2.5947 


0,1 ±6,2 


Kine et al. f2008) 


0.86 


0,74 ± 1,23 


MurDhv et al. f2008) 



versas del amoniaco a 24 GHz. Un cambio en la frecuencia debido a una 
variacion de se puede discernir del redshift cosmologico comparandolas con 
transiciones cuya sensibilidad a esta cantidad fisica sea menor. Para compa- 
rar, se utilizan transiciones rotacionales de moleculas tales como CO, HCO"^ 
y HCN. Por ahora se conoce un solo sistema de absorcion en quasares con 
absorcion en NH3. Agregamos la cota que se obtiene a partir de este sistema 
al conjunto QSO I. La misma se muestra en la ultima entrada de la Tabla 

m 

Otro metodo para acotar la variacion de las constantes fundamentales 
esta basado en la comparacion entre la transicion hiperfina de 21 cm, en ab- 
sorcion, para el hidrogeno {va) y una transicion resonante optica {vb)-, tambien 
del hidrogeno. El c ociente ^ es prop orcional a x = a^Qp"^ donde Qp es el 
factor g del proton (ITzanavaris et al.ll2007l ). Por lo tanto, un cambio en esta 
cantidad resultara en una diferencia en el redshift medido a partir de la Imea 
de 21 cm y el espectro medido a partir de Ifneas de absorcion en el optico, da- 
da po£^_^^^^^^Dado que estamos trabajando en el contexto del modelo 
de lBarrow fc Magueijd (120051 ). la linica constante fundamental que puede va- 
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riares rUe. La Tabla muestra las cotas presentadas por iTzanavaris et al. 



(120071 ) combinando las medidas de los redsliifts obtenidos del optico y del 
radio. A este conjunto de datos lo Uamaremos QSO II. Este metodo tiene 
el inconveniente de que es dificil determinar si ambas llneas (radio y optico) 
fueron originadas en el mismo sistema de absorcion. Por lo tanto, una dife- 
rencia en la velocidad de las nubes de absorcion puede esconder una variacion 
de X. 



A me 



Tabla 4.3: QSO II: Redshift de absorcion y valor de ^'"j con su correspon- 
diente error (en unidades de 10~^), obtenidos de comparar redsliifts de radio 
y de Ifneas moleculares. 



Redshift 




0.24 


1,21 ± 2,10 


0.31 


-0,61 ±4,27 


0.40 


3,22 ±3,15 


0.52 


-2,95 ± 1,05 


0.52 


0,26 ±3,67 



Redshift 




1.78 


-2,59 ±0,90 


1.94 


3,30 ± 0,44 


2.04 


5,20 ± 2,76 


2.35 


-2,54 ± 1,82 



4.3.4. Cotas a partir de datos geofisicos 



La vida media de emisores de particulas (3 de vida media larga ha sido de- 
terminada tanto en medidas en laboratorio, como comparando con la edad de 
meteoritos, la cual se encuentra del analisis del decaimiento de particulas a. 
Las cotas mas fuertes para la variacion de la vida media. A, viene de la com- 
paracion del decaimiento del ^^^Re durante la formacio n del sistema solar y 



el pre sente, flOlive et al.ll2004f ): ^ = (-0.016 ± 0.016). ISisterna fc Vucetich 



(Il990f ) derivaron una relacion entre el corrimiento en la vida media y una 



posible variacion de las constantes fundamentals. Como aquf estamos con- 
siderando solo la variacion de me, la cota que consideramos se obtiene de la 
siguiente relacion: 



AA 

X 



Arue 
'(me)c 



(4.65) 



donde a = —600 para ^^'Re. 
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4.3.5. Cotas de laboratorio 



La comparacion entre frecuencias de diferentes transiciones atomicas a lo 
largo del tiempo es litil para poner cotas fuertes a la tasa de variacion de 
las constantes fundament ales en la actualidad. En particular, la frecuencia 
hiperfina del Cesio se puede aproximar por uqs — 5'cs^«^-Ry-Pcs(«), donde 
gcs es el factor nuclear g, Ry es la constante de Rydberg expresada como 
frecuencia y -Fcs(a) es una funcion adimensional de a y no depende de nie al 
menos a primer orden, mientras que las frecuencias de tra nsiciones opticas 
se pueden expresar como Vopt — R yFja). Varios autores (IBize et al.l 12003 : 
Fischer et al.ll2004l : IPeik et al.ll2004l ) han medido diferentes transiciones opti- 
cas y las han comparado con la frecuencia de la separacion hiperfina del nivel 
fundamental en ^^^Cs neutro. Estas medidas se pueden usar para poner limi- 
tes a la variacion de rhg/ (^e)o- En la Tabla 14.41 se listan cotas obtenidas 
de diferentes experimentos. Estos datos seran comparados con la prediccion 
para me/ (me)^ que se deduce de la ecuacion fl4.54p : 



-Mt + A 



^ Ho 



[m, 



(4.66) 



Tabla 4.4: Relojes atomicos comparados, valor de con su correspondien- 

te error (en unidades de 10~^^), intervalo de tiempo para el cual se midio la 
diferencia, y la referenda. 



Frecuencias 




At[yr] 


Referenda 


Hg"*" and Cs 


0,2 ± 7,0 


5 


Fischer et al. f2004) 


Yb+ and Cs 


1,2 ±4,4 


2.8 


Peik et al. (2004) 


Hg+ and Cs 


0±7 


2 


Bize et al. f2003) 



4.4. Ajuste de los parametros del modelo 

El modelo que estamos estudiando tiene dos parametros libres, cuyos 
valores se pueden estimar buscando el mejor ajuste para el conjunto de datos 
descripto en la seccion anterior. Ajustaremos las cantidades adimensionales 
M/Hq y A/Hq, por el metodo de maxima verosimilitud. Para ello, hacemos 
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una grilla de valores para M/Hq j A/Hqj evaluamos el en cada uno de los 
puntos de la grilla. Para asegurar que el valor del campo sea pequefio para 
todo tiempo, lo cual debe ser cierto ya que hemos aproximado exp[0] ~ 1 + 0, 
agregamos a la expresion del un nuevo termino que tenga en cuenta el 
vinculo dado por la ecuacion fl4.26p . Ademas, tomamos una expresion elevada 
al cuadrado porque el siempre debe ser positivo. La expresion para el 
queda: 



X 



+ 



1 47r Gcj I H, 



Ho 



^Ro 3 c4 



''el 



(4.67) 



Existe una relacion entre uj y M/Hq. Recordemos la Ec. (I4.17p : 



M = 

UJ 



(4.68) 



Reescribimos esta relacion buscando una expresion para peo- Sea N}, el numero 
total de bariones, A^h el numero de nucleos de atomos de Hidrogeno, Ane el 
numero de nucleos de atomos de Helio, y Ne el numero total de electrones. 
Se tiene que, dado que el universo es neutro: 



ATe = ATh + 2ArHe 



(4.69) 



y dado que cada micleo de hidrogeno aporta un barion y cada nucleo de helio 
aporta cuatro bariones. 



ATfe = A^H + 4ArHe = N, + 2Nne. 



(4.70) 



Ahora, sea /ne la fraccion de bariones en forma de nucleos de He con respecto 
al numero total de bariones, es decir: /ne = 4A'He/Afe (ya que cada He tiene 
4 bariones). Entonces 

Ne = Ni,- 2Nue = Nb- 2^Nb = AT, (^1 - . (4.71) 

Por otra parte, la densidad (en masa) de electrones es la densidad numerica 
multiplicada por la masa del electron. Por ende, se tiene que 



M, _ N,m, _ AT, (1 - if) 



V 



V 



Trip 



1 - ^) pJ-^ (4.-2) 
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donde en la ultima igualdad se tuvo en cuenta que 

Pb = —TT- - —TT^ ^ -rT = — • 4.73 

V V V nip 

Escribimos la densidad (en masa) de bariones en funcion de la densidad 
critica de energia: 

Pb = c-^nbPc- (4.74) 

Por lo tanto, 

PeO = 1 TT I pb = 1 - ^ ^bpcC (4.75) 



y dado que pc = se tiene 



TUp \ l J rrtp 



que 



PeO = 1 - ^ — ^b — -7T 4.76 
2 J nip Stt G 



por lo tanto, 



PeO 4 PeO 4 1 / i /hs \ me 3 C 1 

a; = — c = — ^ — ^ =1 — i'b — — ^ 4.77 

M Hi (m\ \ 2 J nip 8n Gfj^^ ^ ' 



o lo que es lo mismo: 



^ = - f 1 - ^1 ^^b^. (4.78) 
Stt V 2 J nip /^jf ^ ^ ^ 



^0 



Por lo tanto, la expresion para el Que vamos a minimizar es la siguiente: 



X 



E 



Am. 



(ti 



1 - 



He 



Ho 



aei 



(4.79) 

La suma es sobre la totalidad de los datos, y tj se refiere al tiempo correspon- 
diente a cada uno de ellos. Este primer analisis estadistico fue realizado antes 
del afio 2008 y por ende, se utilizaron las cotas para Ame/rrie obtenida con 
los quasares de la Tabla HTTl El resultado del analisis estadistico muestra que 
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no se logra buen ajuste para el conjunto completo de datos. Sin embargo, 
se puede encontrar un buen ajuste excluyendo el quasar a 2; = 1,94 de la 
Tabla I4.3[ Con el objetivo de estudiar como afecta cada conjunto de datos 
a los valores obtenidos para los parametros del modelo, repetimos el analisis 
estadistico, excluyendo de a un grupo de datos por vez. Encontramos que 
la cota que proviene de nucleosmtesis es crucial para determinar el valor de 
A/Hq. Por otra parte, el grupo de datos de la Tabla Wl] es importante para 
determinar el valor de M/Hq. Los resultados se muestran en la Tabla l¥]5| 
donde hemos llamado QSO I al conjunto de datos de la Tabla HTTl y QSO II 
al de la Tabla H73l Estas conclusiones no se modifican cuando se utilizan las 
cotas dadas en la Tabla 



datos 




A/Ho 


\^ - 

N-2 


todos los datos 


' ''-"-'-2,02 




1.14 


todos — nucleo 


' ''-"-'-2,02 


'-''-'-'-'-13,3 


1.22 


todos - CMB 






1.16 


todos — decaim j3 


-7 40+^'^^ 


3,60_i5o 


1.13 


todos - QSO I 




3,60_i59 


1.29 


todos - QSO II 


— ^ 70+^'^^ 
J? ' 'J-2,44 


3,60_i49 


0.30 


todos — relojes 




3,60_i5o 


1.44 



Tabla 4.5: Parametros para el modelo que da el mejor ajuste de los datos. El 
valor para el parametro M/Hq esta en unidades de 10~^, y el valor para el 
parametro A/Hq esta en unidades de 10~^^. 



Dado que el modelo esta descripto en terminos de la constante de acopla- 
miento u, derivamos su valor para el mejor ajuste, a partir de los resultados 
anteriores. Partimos de la siguiente relacion: 



^ = ± f 1 _ AeU, f 4 ) (4.80) 

Stt V 2 J nip V^o 



y usamos los siguientes valores numericos: 
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G = 6,67 X lO-^cm^g-^s-^ 

^ = ££!!i = 0,0418 

" (0,73)2 

^ = 1836,15. 



(4.81) 



Falta estimar el valor de /ne- Para ello, hacemos las siguientes considera- 
ciones: 



^^b /He 

_ 4iVHe _ 4 _ 1 

iVjj 

Por otra parte, obscrvacionalmente se determina que el cociente entre la 
masa en He y la masa en H es 

y de la tabla periodica, se determina que 

Il^Si = = 3,9712. (4.84) 

mn 1,0079 ' ^ ' 

Entonces, /ne = 0,1946788, y por ende: 

^ _ /He ^ 0,9027. (4.85) 

Por lo tanto, dado que del ajuste de los parametros del modelo se encuen- 
tra que 

M --.+2,10^,^-6 



-2 = -7,30:^;i^ X lO-'' (4.86) 



el valor para ujG/c^ es 
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= -0,336tSl (4.87) 



donde hemos encontrado los errores de uG/c^ propagando los errores de 
M/Hq, de la siguiente manera: 



Esto lo hemos podido hacer dado que el intervalo a la alrededor del mejor 
ajuste para el parametro M/Hq no contiene al cero. 

La cantidad uoG/c^ es adimensional. Para comparar con el resultado de 
Barrow & Magueijo, debemos usar cantidades adimensionales ya que estas no 
dependen de las unidades. Barrow & Magueijo utilizan unidades naturales, 
por eso presentan sus resultados como G\uj\ > 0,2. En unidades usuales, este 
resultado se traduce en Glul/c^ > 0,2, que es el que usamos para comparar 
con nuestros resultados. 

El contorno de confianza que resulta del analisis estadistico puede verse 
en la figura lCTl El corte que se observa en el mismo se debe a que la condicion 
impuesta en el diverge para M = 0. 

En el afio 2008, fueron actualizados algunos de los datos de la seccion 
14.3.3^ que se muestran en la Tabla W72[ Adicionalmente, se ha publicado una 
cota para la variacion de la masa de electron obtenida con un nuevo metodo 
que utiliza las Ifneas de las transiciones inversas del NH3. Hemos repetido 
el ajuste teniendo en cuenta dichas actualizaciones. El analisis muestra que 
no se obtiene un buen ajuste al considerar el conjunto completo de datos. 
Sin embargo, se logra un buen aiust e cuando se descartan t res sistemas de 



M 



absorcion en quasares del grupo de iTzanavaris et al.l (120071 ) (el que se en 



cuentra a z = 1,94, que habia sido descartado tambien en el ajuste anterior, 
el que esta a z = 1,78, y la entrada 4 de la Tabla 14.31 de la izquierda, con 
z = 0,52). La nueva cota para la variacion de la masa del electron obtenida 
con los datos de WMAP5 no va a cambiar el ajuste de M/Hq debido a que 
las barras de error para las cotas que provienen de CMB (y BBN) son mucho 
mayores que las correspondientes a los demas datos, como puede verse en la 
Fig. dH). 

El resultado del nuevo ajuste, realizado con 16 datos, se presenta en la 
Tabla 1121 Observar el cambio en el valor de la cota para el parametro M/ Hq. 
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-2e-13 I 1 1 1 1 1 1 1 1 — 

-1.4e-05 -1.2e-05 -1e-05 -8e-06 -6e-06 -4e-06 -2e-06 2e-06 



Figura 4.1: Contorno de confiabilidad para los parametros del modelo de 
Barrow & Magueijo. 



Este valor es mas consistente con cero que el valor obtenido en el ajuste 
anterior (ver Ec. fl4.86p ). Lo que cambio entre un ajuste y el otro es el grupo 
de datos QSO I. El nuevo conjunto de datos es menos favorable para la 
variacion de me. Esto se traduce en un valor mas consistente con cero para 
uno de los parametros que controla la amplitud de la funcion (j)(t), M/Hq, 
responsable de la variacion de mg. 





A/Ho 


- 

Af-2 




'-'5'-' "-"-1,32 


0.87 



Tabla 4.6: Parametros del mejor ajuste considerando los datos actualizados 
al 2008. Las unidades de M/H^ son 10"^ y de A/Hq son 10^^^ 



83 



A partir de la cota para M/Hq, encontramos la cota para Gu/c 



(4.89) 



El contorno de confianza que resulta de este analisis estadistico puede verse 
en la figura I4.2[ Nuevamente, el corte que se observa es debido a que la 
condicion impuesta en el diverge para M = 0. 



4e-13 




2e-13 



-2e-13 



-2e-06 



Figura 4.2: Contorno de confiabilidad para los parametros del modelo de 
Barrow & Magueijo, segiin los datos actualizados. 



4.5. Principio de Equivalencia Debil 

De manera complementaria, se pueden obtener cotas para el parametro u 
del modelo a partir de los limites impuestos por los experimentos que testean 
la validez del Principio de Equivalencia Debil. Para obtener una expresion 
que podamos comparar con dichos datos, consideremos las variaciones espa- 
ciales del campo escalar 0. Dado que estudiaremos las implicancias de dichas 
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0.01 


-0.01 
-0.02 
-0.03 
-0.04 
-0.05 
-0.06 
-0.07 



1 — 1 — I — 1 — 1 — 1 — 1 — I — 1 — 1 — 1 — 1 — |- 

Regimen: rad + mat 

Regimen: rad + cte cosmol. 

— iHi datos quasars + 

BBN + CMB x 



_i_ 




1e-15 



1e-10 



1e-05 



100000 



1e+10 



Figura 4.3: Solucion del modelo de Barrow & Magueijo correspondiente a los 
parametros del mejor ajuste. Se presentan tambien los datos utilizados para 
el ajuste, con sus correspondientes barras de error. 

variaciones a nivel local, despreciaremos las variaciones temporales de 0. A 
partir de la ecuacion de Euler-Lagrange: 



d 



dXf, \d{d(j)/dXf,) 
y teniendo en cuenta que 



(4.90) 



(4.91) 



llegamos a la siguiente ecuacion: 



rrieC 

□ = Ue 

UJ 



(4.92) 



85 



Cuando despreciamos las variaciones temporales, 00 = — V^0, y por lo tanto, 
la ecuacion se transforma en: 

2 2 2 

VV = = -pe = -Cp 4.93 

UJ UJ UJ 

donde ( = Pe/ P- 

Para encontrar la solucion de esta ecuacion, usaremos el metodo de la 
Funcion de Green, y consideraremos que toda la masa (M) esta concentrada 
en el centre de masa del cuerpo, de manera tal que p = M6^{r). Dado que 

V^G(x, x') = 6(x - x') G(x, x') = --^^ — r (4-94) 

4:71 \X — X'\ 

la solucion de la ecuacion fl4.93p es 

1 1 

= (M- (4.95) 

Att UJ r ^ ^ 

Consideremos el caso de una partfcula de prueba en presencia de un cuer- 
po mas masivo como, por ejemplo, el sol o la tierra. En este caso, ademas de 
la fuerza gravitatoria, la partfcula experimenta una fuerza debida a las varia- 
ciones espaciales del campo escalar generado por los electrones presentes 
en el sol (o en la tierra). 

Supongamos que la partfcula de prueba tiene una masa total y que la 
masa en electrones es m^. Planteamos la segunda Ley de Newton: 

^ F = -mtV$ - mlc^Wcj) = rritd (4.96) 

de donde se deduce la siguiente expresion para la aceleracion que sufre la 
partfcula de prueba: 



-V<l> -c'Vcj) = -V$ - Ctc'V0 (4.97) 

nit 



donde $ es el potencial newtonianoH 



^Hemos considerado que la fuerza inducida por (j> sobre la materia viene dada por 
mfC^Wcj) ya que la densidad lagrangiana de inter accion es 

^mt = -meC^tpipc/) = -meC^n{r)(f){f) = -TOeC^(5(r - fo)0(f) 

d^xJf^t = -meC^Hro) (4.98) 

y por lo tanto, el potencial de interaccion en tq, inducido por 0, es —nieC^ (j){ro) . 
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El campo gravitatorio se deriva de la funcion potencial $ = —GM/r. De 
est a manera: 

V$ = — (4.99) 

Por otra parte, el campo escalar que generan los electrones en el sol 
(ec. I4.95p . Ueva a la siguiente expresion: 



1 c2 1 

Vc/) = CsM-f 

Air uj 



(4.100) 



En conclusion, la expresion para la aceleracion que sufre una particula de 
prueba en el campo gravitatorio del sol, y en el caso de que el campo (j) tenga 
variaciones espaciales, es: 



GM 



(4.101) 



Vemos que la aceleracion depende de la relacion entre la cantidad de electro- 
nes y la cantidad de bariones, es decir, depende de la composicion qui'mica 
de los cuerpos. 



4.5.1. Experimentos de Eotvos 

Dados dos cuerpos A j B, que caen con aceleraciones a a y Sb respecti- 
vamente, se define el parametro de Eotvos de la siguiente manera: 

r^{A,B)^^^^^^ (4.102) 

En nuestro caso de interes, cuando la aceleracion tiene la expresion de la 
Ec. f l4.10ip . el parametro 1] toma la forma: 



^ = l^T^C. ^ , ^ ^^Cs (Ca - Cb) (4.103) 

AttGuj 1 + ^ ((^ - (s) 47r Gu 

Este parametro se ha medido en distintos experimentos cuyos resultados 
y referencias se detallan en la Tabla 14. 7[ Por lo tanto, se puede obtener una 
cota sobre el parametro c^/Guj. Sin embargo, necesitamos saber el valor del 
parametro C, de los elementos con los cuales se realiza el experimento. 
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Sea n la densidad numerica de atomos de la substancia monoatomica 
analizada, dc mimcro atomico Z, y peso atomico P. Recordamos que Z re- 
presenta el numero de protones en el niicleo (dado que el atomo es elcctrica- 
mente neutro, esto equivale al numero de electrones), y P es el promedio de 
la masa de los isotopos de dicho elemento, pesados segiin la abundancia con 
la cual aparecen en la tierra, y en unidades de masa atomica (uma), que se 
define como la doceava parte de la masa de un atomo de ^^C. Por lo tanto, 
la expresion para el parametro ( es: 

^ p n P ^ ' 

Dado que 1 uma = 1,661 x 10"^'' g = 931,57 MeV, y que la masa del 
electron es igual a 0,511 MeV, el parametro C se escribe como: 

^ Z X 0,511 

donde P se escribe en umas, como aparece en la tabla periodica de los ele- 
ment os. 

Si en lugar de ser monoatomica, la substancia esta compuesta de una 
mezcla de distintos elementos, el parametro C es: 

C ^ (4.106, 

donde Ni cs el numero de atomos de cada elcmcnto. Si tcncmos una cierta 
cantidad de masa Mj del elemento i, y esta cantidad es un porcentaje pi de 
la masa total M de la substancia, tenemos que: 

N.= ^ = '^ (4.107) 



Por lo tanto, el parametro ( queda: 

' Pi 7 



C = — — = 2^ -B-Zime (4.108) 



o lo que es lo mismo: 



Pi Z, 0,511 
Pi 931,57 

donde Pj se expresa en umas. 



(4.109) 
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Haremos un ajuste del parametro c^/Guj con los datos de la Tabla HTl 



Materiales 


r]{A,B) X 10" 


Referencia 


en $ del: 


Al-Au 


1,3 lb 1,0 


Roll et al. ri964) 


sol 


Al-Pt 


0,0 ±0,1 


Braffinskii & Panov (1972) 


sol 


Cu-W 


0,0 ±4,0 


Keiser & Faller fl982) 


sol 


Be-Al 


-0,02 ±0,28 


Su et al. ri994) 


tierra 


Be-Cu 


-0,19 ±0,25 


Su et al. ri994l 


tierra 


Si/Al-Cu 


0,51 ±0,67 


Su et al. ri994) 


sol 



Tabla 4.7: Resultados de experimentos de Eotvos. 



Como mencionamos anteriormente, es necesario calcular el parametro C, 
para los elementos utilizados en dichos experimentos. Los resultados aparecen 
en la Tabla En el experimento realizado con Si/Al-Cu, la muestra de 
Si/Al tiene una masa de 10,0432 g (62,7% Si - 37,3% Al). Para el caso del 
Sol y de la Tierra, necesitamos la composicion quimica detallada, la cual se 
muestra en la Tabla 14. . 



elemento 


Z 


P 


Cx 10^ 


Al 


13 


26.98 


2.643 


Au 


79 


196.97 


2.200 


Pt 


78 


195.08 


2.193 


Cu 


29 


63.55 


2.503 


W 


74 


183.84 


2.208 


Be 


4 


9.01 


2.435 


Si 


14 


28.09 


2.734 


Si-Al 






2.700 


Sol 






4.704 


Tierra 






2.568 



Tabla 4.8: Valores del parametro C,. 



Para encontrar una estimacion del parametro /Guj, hacemos una grilla 



■^La composicion quimica del Sol y de la Tierra, fueron tomadas respectivamente de 
http:/ / astronomy, swin. edu. au/ cosmos /C/ Chemical+ Composition 
http://seds.lpl.arizona.edu/nineplanets/nineplanets/earth.html . 
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Element o 


Pi (solar) 


Pi (terrestre) 


Hidrogeno 


73% 




Helio 


25% 




Oxigeno 


0.80 % 


29.5 % 


Carbono 


0.36 % 




Hierro 


0.16 % 


34.6 % 


Neon 


0.12% 




Niquel 




2.4% 


Nitrogeno 


0.09 % 




Silicio 


0.07% 


15.2% 


Magnesio 


0.05 % 


12.7% 


Sulfur o 


0.04% 


1.9% 


Titanio 




0.05 % 


Otros combinados 


0.04% 





Tabla 4.9: Composicion quimica del Sol y de la Tierra. Los valores son por- 
centajes en masa. 



de valores de este parametro y calculamos el x^, cuya expresion es 

. 2 



X 



E 



erf 



(4.110) 



donde rji y cTj son los datos y errores de la Tabla 14.71 Al minimizar esta 
cantidad, encontramos el mejor ajuste para el parametro c^/Guj. A Icr, el 
resultado es 



(7,5 ±66,6) X 10" 



(4.111) 



Como este intervalo contiene al cero, no podemos liacer una propagacion de 
errores para encontrar una cota para el parametro que nos interesa, que es 
Guj/c^. Sin embargo, podemos dar una cota inferior para el mismo. Efectiva- 
mente, dado que 

(4.112) 



Gu 



< 74,1 X 10" 



llegamos a que el limite impuesto por el Principio de Equivalencia Debil es 

Guj 



> 1349. 



(4.113) 
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4.6. Validez del modelo 



En el paper original (IBarrow fc Magueiidl2005l ) se hacen algunas aproxi 



maciones acerca de la evolucion del factor de escala, por lo que la solucion 
del modelo se simplifica. En esta tesis se presentan soluciones mas precisas 
de la evolucion del factor de escala, y en consecuencia, se obtienen soluciones 
mas exactas para el campo escalar que induce variaciones en la masa del 
electron. Otra mejora en las soluciones proviene de no despreciar la prime- 
ra constante de integracion, la cual es de suma importancia en el Universo 
temprano. Efectivamente, una vez que se Integra la ecuacion (14.541) se puede 
escribir 

= -M - ^ = -M {t + 8,47 x lO^r) (4.114) 

donde se ban usado los valores del mejor ajuste de los parametros M / y 
A/Hq, y el valor h = 0,73. Notar que el segundo termino del miembro dere- 
cho domina en el Universo temprano, en particular durante la nucleosfntesis 
primo rdial. 



En lBarrow fc Magueijd (120051 ) se presenta una cota de G\uj\/c'^ > 0,2. En 
dicho trabajo, este resultado se obtiene a partir de la precision de las cotas 
para la variacion de mg provenientes de los quasares, pero no se consideran las 
medidas individuales, con sus valores medios y errores. Ademas la expresion 
que relaciona Ame/(me)o y el parametro G\uj\/c'^ es muy simplificada. En 
cambio, en esta tesis se ban utilizado todos los datos disponibles que ponen 
limites a la variacion de la masa del electron, en distinto s tiempos cosmologi- 



cos. H ay una diferencia entre los analisis, dado que en iBarrow &: Magueijo 



(120051 ) se considera como cota el valor Ame/(me)o < 10~ , mientras que en 
esta tesis se ha considerado cada dato con su correspondiente error. Al ob- 
tener cotas a partir de experimentos que testean la validez del Principio de 
Equivalencia Debil, ellos obtienen una cota G|ci;|/c^ > 10^, comparable con 
la obtenida en esta tesis a partir de la misma consideracion. Sin embargo, 
al comparar con la cota obtenida al restringir la variacion temporal de la 
masa del electron a lo largo de la historia del Universo con los datos exis- 
tentes hasta el momento, se obtiene que 1,023 < Guj/c^ < 4,368 (intervalo 
de confianza a 1 a). (Antes de la actualizacion de datos del 2008, habfamos 
obtenido —0,615 < Goj/c^ < —0,045 a 3 a). De esta comparacion U e gamos 



a la conclusion de que el modelo propuesto en IBarrow fc Magueijd (120051 ) 
para la variacion temporal de la masa del electron es inconsistente y debe ser 
descartado. 
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Un ultimo comentario acerca del signo de uj. En el ajuste con datos ante- 
riores a 2008, a 2cr, el valor de uj es negativo. Esto no debe resultar sorpren- 
dente. De hecho, terminos cineticos negativos en el Lagrangiano ya ban sido 
considerados por otros autores, por ejemplo, en modelos de k-e ssence con 
una componente de energia fantasma (pbantom) ( Caldwelil I2OO2I ) . De todas 
maneras, con los datos actualizados esto no sucede. 
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Capitulo 5 
Radios Planetarios 



Consideremos un sistema fisico en el laboratorio, libre de la accion de la 
gravedad. Desde el punto de vista de la termodinamica, el estado del sistema 
se puede describir con una funcion potencial energta interna U que depen- 



variables extensivas. Las derivadas primeras del potencial definen las varia- 
bles intensivas temperatura T, presion P, y potencial qmmico Quedan 
establecidas, entonces, las Uamadas ecuaciones de estado. En particular 



Cuando sometemos el sistema a la accion de la gravedad, cada particula 
del "gas" experimenta la fuerza gravitatoria, que se suma a las fuerzas inter- 
nas de interaccion con las otras particulas. Sin embargo, la dependencia de la 
energia interna U con las otras variables termodinamicas no cambia cuando 
el sistema esta en un campo gravitacional y las ecuaciones de estado son las 
mismas. 

Estudiemos ahora un gas autogravitante como sistema termodinamico. 
Supongamos que el sistema tiene simetn'a esferica, esta en equilibrio hi- 
drostatico, y ademas no hay fuentes ni sumideros de masa, es decir, la masa 
se conserva. Para encontrar su estructura (cantidades termodinamicas en 
funcion de la coordenada radial) debemos resolver el siguiente sistema de 

^Una variable intensiva es una cantidad fisica cuyo valor no depende de la masa de la 
substancia. 



de de la entropfa S, el volumen V, y el niimero de particulas N. Estas son 




(5.1) 
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ecuaciones: 



dP Gm{r)p{r) 

= 7^ — ^^-^^ 

dm{r) 



dr 



A7rp{ry. (5.3) 



Es necesario, ademas, utilizar la ecuacion de estado que relaciona la presion 
con la densidad. Esta se obtiene para el sistema particular de estudio, en el 
laboratorio, cuando esta libre de la accion de la gravedad. 

El objetivo de este capitulo es encontrar cual es la dependencia de las 
cantidades termodinamicas con las constantcs fundamcntales, y como se ve 
afcctada la estructura de un planeta si las constantes fundamcntales toman 
un valor diferente del actual. En particular, analizaremos como cambia el ra- 
dio del planeta. Estudiaremos luego la estructura interna de tres cuerpos del 
sistema solar: La Tierra, el planeta Mcrcurio, y la Luna, y buscaremos una 
relacion teorica cntre la variacion del radio y la variacion de las constantes 
fundamcntales. Para estos sistemas contamos con cotas observacionales para 
posibles variaciones de sus radios entre la epoca de formacion del sistema 
solar y la actualidad. Con los Kmites observacionales, y las expresiones teori- 
cas obtenidas, se dara una cota para la variacion "local" de las constantes 
fundamcntales en los ultimos 4.000 millones de aiios. 



5.1. Deduccion teorica 

Partimos del hamiltoniano de un sistema general de electrones e iones, 
donde despreciamos la gravedad pero consideramos la interaccion electro- 
magnetica entre las particulas. Identificando el hamiltoniano con la energia 
interna, deducimos la dependencia con las constantes fundamcntales de todas 
las cantidades termodinamicas. 

Luego, consideramos un planeta como sistema autogravitante, y estudia- 
mos como la gravedad infiuye en que la estructura del mismo no sea invariante 
frente a un cambio de escala. Asimismo, analizamos como la variacion de las 
constantes fundamcntales afecta la estructura del planeta, actuando a traves 
de una modificacion efectiva de la constante de gravitacion de Newton. 

Finalmente, deducimos una expresion para la variacion del radio del pla- 
neta cuando las constantes fundamcntales varian en el tiempo. 
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5.1.1. Hamiltoniano y cambios de escala 

Sistema sin gravedad 

Consideremos el hamiltoniano (en unidades cgs) de un sistema cualquiera 
de electrones y nucleos, tal como moleculas o solidos, en ausencia de gravedad 

^ = E^ + -El^--E%^ (5.4) 

en donde Rij = |rj — rj|. El segundo termino corresponde a la repulsion 
electron-electron, y el tercero a la interaccion electron-micleo. Hemos despre- 
ciado la energia cinetica de los nucleos, que es muy pequeiia comparada con 
la de los electronesJl 

Si en el hamiltoniano introducimos un factor de escala para las longitudes 

r = /?x P = y//? (5.5) 

el hamiltoniano toma la forma 



2 



(5.6) 



Si introducimos la represent acion de coordenadas y = —iKV, podemos 
elegir (3 de modo tal que los coeficientes de la energia cinetica y la energfa 
potential sean iguales 

IT^ P = — = = as (5.7) 



donde es el radio de Bohr. Con este factor de escala, el hamiltoniano toma 
la forma 

H = —Ha (5.8) 



^Es la aproximacion dc Born-Oppcnhcimer. Las correcciones por movimicnto de los 
nucleos son del ordcn 0{me/M). 
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en donde el hamiltoniano adimensional, independiente de las constantes fun- 
damentales, es: 



(5.9) 



Las ecuaciones (15.51) y (15.81) muestran que en condiciones adiabaticas y 
en ausencia de gravedad la variacion de las constantes cambia las longitudes 
y las energias con los factores de escala 

r = aBX E = — Ea- (5.10) 

as 

Como consecuencia, el volumen V cambiara como 



V = alVA. 



(5.11) 



Reescaleo de variables termodinamicas 

Para obtener como cambian otras cantidades termodinamicas, recordemos 
la definicion de la Entalpfa 7i: 



'}i = U + PV = n{S,P,N). 
Por tener dimensiones de energia, su cambio frente al reescaleo es 



(5.12) 



Ti. = — "Ha- 

as 



Por otra parte, 



p_dn_ e^mUdVA _ ^BHa _ 

~ dV ~ aB OVa dV ~ a% OVa ~ ^' 



"B 



La compresibilidad adiabatica variara como 

1 



(5.13) 



(5.14) 



(5.15) 



Supongamos que las constantes fundamentals dependen del tiempo. En 
consecuencia, el factor de escala qb tambien tendra una dependencia tempo- 
ral 



r(t) 



aB(t)x 
aB(to)x 



(5.16) 
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donde to el valor del tiempo hoy. Sin embargo, x no depende del tiempo, 
porque la dependencia temporal esta en las constantes fundamentals. Por 
lo tanto, al hacer el cociente queda 

r = /3ro (5.17) 



donde hemos redefinido j3 = aB{t)/aBo, con a^o = ^^(^o)- Tambien defini- 

2 

mos por conveniencia 7 = Los submdices representan las cantidades 

termodinamicas cuando el valor de las constantes fundamentals es el actual. 
De este modo, obtenemos las relaciones de escala 

r = f3ro V = f3'Vo P = p-'po ..^ 

Las relaciones termodinamicas se dcducen a partir del hamiltoniano del 
sistema en ausencia de gravedad, por lo tanto, estas relaciones de escala 
siguen valiendo aiin cuando estudiemos el sistema incluyendo la fuerza gra- 
vitatoria. 



Sistema autogravitante 

Supongamos ahora que tenemos en cuenta los efectos de la gravedad en 
el sistema. El hamiltoniano sera 



H = yi^ + ly^-ly^-Gy^^ (5.19) 

donde los indices i,j hacen referenda a los electrones, a hace referenda a 
los niicleos, y a, /3 hacen referenda a cualquicr particula. En la ultima suma, 
los terminos que mas contribuyen son los que ticnen el doble producto de 
las masas de los niicleos, debido a que la masa de cualquier niicleo es mucho 
mayor que la masa del electron. Por ello, el hamiltoniano se puede aproximar 
como 

donde nuevamente a y 6 hacen referenda a los niicleos del sistema. 
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En este caso, el cambio de variables 

r = /?x P = y//3 (5.21) 

no nos permitira expresar al hamiltoniano como un factor de escala por un 
hamiltoniano adimensional, independiente de las constantes fundament ales. 

Sin embargo, las relaciones de escala para las cantidades termodinamicas 
se mantienen, y las usaremos en las secciones siguientes, dado que nos intere- 
san las ecuaciones de estado, que se calculan a partir del hamiltoniano libre 
de gravedad. Por ultimo destacamos que el subindice no hara referenda a 
las cantidades "hoy", sino a las cantidades reescaleadas con un factor que 
tiene en cuenta el valor de las constantes al tiempo t, y el valor actual. Esto 
sera explicado mejor en las secciones subsiguientes. 



5.1.2. Ecuacion de equilibrio hidrostatico 

Retomemos el estudio de un sistema autogravitante, con simetn'a esferica, 
en equilibrio hidrostatico, en el cual la masa se conserva. Veamos como se 
modifica la ecuacion de equilibrio hidrostatico cuando reescaleamos las can- 
tidades termodinamicas. Supongamos que las constantes fundamentales a y 
rrie pueden tener valores arbitrarios, diferentes de sus valores actuales. Usare- 
mos la hipotesis adiabdtica, que establece que en presencia de una variacion 
adiabatica de las constantes fundamentales, la forma de las ecuaciones no 
cambia. 

Con esta condicion, la ecuacion de equilibrio hidrostatico toma la forma 
habitual: 

= -G—-^f-^ (5.22 

dr 

en donde 

m{r) = Att p{r')r'^dr' (5.23) 







En estas ecuaciones tanto a como me pueden variar lentamente con el tiem- 

poEl 

El sistema esta en equilibrio, tiene un radio total R y una estructura 
radial en las cantidades termodinamicas dada por las funciones 



"Lentamente" significa que la escala de tiempo de variacion es mucho mas grande que 
las escalas ti'picas en el sistema solar. 
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P{r) Presion 

p(r) Densidad 

d{r) = r Distancia desde el centre hasta cierta capa 

Ks{r) Compresibilidad adiabatica (5-24) 



donde r es la coordenada radial que va desde a i?. 

Supongamos que r = asx, es decir, medimos r en unidades de radios de 
Bohr aB- Si hacemos un cambio de escala en el sistemfl Y ahora medimos 
en unidades de a^o, estamos cambiando la parametrizacion, y = asox. 
El cambio de escala queda definido por r = ^?"o = P'^o- Las variables 
termodinamicas cambian segiin 

Pi^) = ^Po{ro) p{r) = (3-^Po{ro) K = ^Ko{ro). (5.25) 

P 1 

Cuando reescaleamos las longitudes, d[R) = Rq, tal que R — I3Rq. Sin 
embargo, esto no significa que el nuevo radio del planeta sea Rq. Esto se debe 
a que el hamiltoniano con gravedad no es invariante de escala, y por lo tanto, 
la configuracion de equilibrio del planeta puede ser diferente si las constantes 
fundamentals tienen un valor distinto del actual. Esto implica que la nueva 
configuracion de equilibrio es tal que el nuevo radio es i?o, y la estructura del 
planeta es 



P(ro) Presion 

p(ro) Densidad 

c?(ro) 7^ ro Distancia desde el centro hasta cierta capa 

'^s'(^o) Compresibilidad adiabatica (5.26) 



donde tq es la coordenada radial que va desde a -Rq- 

Veamos como cambia la ecuacion de equilibrio hidrostatico cuando hace- 
mos este cambio de escala. Vamos a hacer la aproximacion de que la masa se 
conserve^, es decir la masa m en un punto r medida en el sistema de coor- 
denadas S es la misma que la masa mo en el punto tq medida en el sistema 

* Similar a la expansion del Universo. 

^Dado que la masa del electron es mucho menor que la masa del proton y de bariones 
mas pesados, el cambio en la masa por el reescaleo es despreciable. 
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de coordenadas Sq, donde r = /?ro. Es decir, 



m(r) = mo(ro). (5.27) 
Asimismo, relacionamos la densidad en uno y otro sistema de coordenadas 



m(r) 1 m(r) 1 mo(ro) 1 / x 
^ = ^ /Q 3 = m A/'i 3 = 75^ . /q 3 = T^^PoO^o • 5.2^ 



Reemplazando en la ecuacion de equilibrio hidrostatico 

y simplificando queda 

dPgiro) _ G mo{ro)po{ro) _ ^ mo(ro)po(ro) 
dro 7/? rg ' 

Vemos que la ecuacion de equilibrio hidrostatico mantiene su forma si 
reescaleamos la constante de gravitacion de Newton de la siguiente manera 

Ge/ = = ^ = G (^) . (5.31) 

Esta Gef depende del valor de la constante de estructura fina, pero no 
depende del valor de la masa del electron. Asi, la estructura de equilibrio del 
sistema, solucion de la Ec. fl5.30p . sera diferente si a varia, mientras que si 
nie varia solo se vera un cambio en el tamaiio del sistema debido al reescaleo, 
pero ningiin cambio intrinseco. Interpretemos esto ffsicamente: 



Cambio de mg 

El cambio en la masa del planeta por un cambio en rrie es despreciable. 
Bajo esta consideracion, G no cambia, y tampoco lo liace el radio, mas alia de 
un reescaleo. El radio del planeta es independiente de su masa. 



Cambio de a 

Si a varia, cambia la fuerza electromagnetica. El sistema esta menos li- 
gado si a es mas pequeiia, y dado que la repulsion entre los micleos da la 
compresibilidad del material, este es mas compresible. El radio del planeta 
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es menor, debido a que la gravedad se ve resaltada; es el mismo efecto que 
se producin'a si G fuese mayor. Es decir, se obtiene una G efectiva, tal que 
Gef > G. 

5.1.3. Cambios en el radio 

Queremos determinar como se modifica el radio de un planeta si las cons- 
tantes fundamentales varian. Es decir, queremos comparar Ro (el radio del 
planeta si las constantes fundamentales tienen el valor actual), con Rq, que 
es el radio del planeta si las constantes fundamentales tienen otro valor, pero 
medido en las mismas unidades que R^. Por lo tanto, la cantidad de interes 
sera ^(M) = Rq — Rq (ambas en unidades de a^o) y no -Rq — -R (que estan 
en distintas unidades). 

Para resolver la ecuacion (15.301) conviene usar la masa como variable 
independiente. Esto tiene la ventaja adicional que la masa es invariante en 
el tiempo 

dmo(ro) = po{ro)'i7crldro. (5.32) 

Sustituyendo en (15.301) e integrando entre una capa (donde la masa acumu- 
lada desde el centro sea m) y la superficie 



° Gef f^'^ modmo 

47r L, rn 



Po{m) = - I dPoinio) = I (5.33) 

donde hemos usado la condicion de contorno P = en la superficie. Esta 
es una solucion formal, dado que ro(mo) es desconocido. La funcion Presion 
reescaleada es Po^m). 

Definamos la siguiente funcion 

1 /"^^^ modrriQ 
A' 





^Para entender esto, hacemos la siguiente analogi'a con la expansion del Universo y el 
tamafio de los ciimulos de galaxias: 

El tamafio ffsico (no comovil) de un ciimulo de galaxias crece con la expansion del 
Universo. Supongamos que queremos comparar el tamafio tipico de un ciimulo actual con 
el tamafio de los primeros ciimulos del Universo, para determinar si hay algiin tipo de 
evolucion intrfnseca. Al mirar los ciimulos lejanos, los estamos viendo como eran en el 
pasado, cuando el factor de expansion era a{t), menor que el actual cq. Para comparar 
el tamafio en aquella epoca con el tamafio hoy, tenemos que reescalear el tamafio del 
ciimulo lejano, de manera tal de simular que el Universo ya se expandio, y luego comparar 
con el tamafio hoy. En otras palabras, se debe medir en las mismas unidades para poder 
comparar. 
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es decir 



Po{m) = G,ff{m). 



(5.35) 



^Que representan las cantidades con submdice 0? 

Si hubiesemos partido de un planeta de masa M^, con distribucion de 
presion Po(^o) Y de densidad Pq^tq), y hubiesemos reescaleado con /? tal que 
solo varie me, liubieramos llegado a la distribucion P(r) y p(r) en equilibrio 
hidrostatico, y al radio R, tal que R = (3Rq. En este caso, las cantidades con 
subfndice representan la estructura del planeta hoy. 

En cambio, si van'a a no podemos partir de {P,p,R), hacer el mismo 
reescaleo y llegar a {Pq,pq,Ro) en equilibrio hidrostatico. Las cantidades 
con subfndice solo son un reescaleo, pero no representan la estructura del 
planeta hoy. 

Al resolver la ecuacion de equilibrio hidrostatico reescaleada, tenemos la 
presion reescaleada: 

Po(m) = Geffim). (5.36) 

Si a varia, Gef 7^ G, mientras que si solo rUe varia, Gg/ = G. Ademas, si 
Gef = G, Pq coincide con la presion actual Pq- En cambio, si Gef 7^ G, Pq es 
la presion reescaleada, pero no coincide con la distribucion de presion hoy. 

Supongamos que tenemos un cambio en Pq dado por un cambio en Gef 
(punto a punto en la coordenada masa). En este caso, 

APo AGe/ . -n Pq — Pq G — Gef 

—— ~ —-^ sigmfica — ^ ~ — ^. (5.37) 

rQ Kj Pq ^ 

donde hemos despreciado el carnbio en f{m) por un cambio en Gef, de manera 
analoga al tratamiento de IWill (Il98lh . 



Vale aclarar que en cada punto las relaciones termodinamicas para las 
distribuciones Pq y Pq son las mismas, por lo que en ambos casos vale 

= -V.dPo- ^'-''^ 

Si la presion reescaleada cambia por variacion de las constantes funda- 
mentals, produciendo un APq, el cambio en el volumen sera AVq, que sera la 
diferencia entre el elemento de volumen hoy y el elemento del volumen en 
coordenadas reescaleadas. 

Sea la transformacion 
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— # — * 

foi'iTT') = ro{m) + ^(m) 



(5.39) 



que mide el desplazamiento de una capa de masa por variacion de las cons- 
tantes fundament ales. Entonces ^(M^) es la diferencia entre el radio real hoy, 
y el radio reescaleado en el pasado, si las constantes toman otro valor. 

Sea la transformacion T : W C — > dada por la Ec. (I5.39p . Escrita 
en componentes es 



Xo = Xo+ ^1 

yo = yo + 6 

= + ^3. 



(5.40) 



El Jacobiano de T mide como se distorsionan los voliimenes por la trans- 
formacion. Asi, Vq se transforma en Vq = JVq, con 



d{xo,yo,zo) 


= det 


{xo, 2/0, ^0) 





dxo 



dxn^ 



dxo 



V 



Ml l _|_ ^2 9|3 

dyo dyo dyo 

O^l OK2 _|_ <7|3 

920 dzo 



dzo / 



J= -}--^^-^--> =det ^ 1 + ^ ^ ^l + VQ.f+0(2) 

(5.41) 

donde hemos considerado que la transformacion es diferencial. De esta ma- 
nera, 

AVo Vo - Vo 



(5.42) 



Vo Vo 

donde AVq significa la diferencia entre el volumen real y el reescaleado (a una 
dada masa acumulada m) y Vq significa el volumen reescaleado. Suponiendo 
que ^ es radial, y tomando la divergencia en coordenadas esfericas. 



De 



Uegamos a que 



Ko 



ro dro 



1 dVo 
VodPo 



APn 



1 AV^o 

Ko Vq 



1 1 djrlO^ 
Korg dro 



(5.43) 



(5.44) 



(5.45) 
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por lo que 



De aqui 



Entonces 



APn 



Aa 



e(ro) 



AGef 1 



1 

/to (iro 



(5.46) 

(5.47) 
(5.48) 



En su estudio acerca de consecue ncias obser vables de posibles violaciones al 
Principio de Equivalencia Fuerte, IWilll (Il98ll ) presenta una expresion equiva- 
lente a esta, aunque en otro contexto y deducida a partir de consideraciones 
diferentes. 

Escribiendo ( 15.48^ en funcion de la masa, y evaluando en Mt, llegamos a 



e(MT 



AG 



1 



Ro 



(5.49) 



donde Rq es el radio reescaleado, que lo aproximamos por el radio actual i?o, 
Poir') es la distribucion real de presion hoy, y ko('^o) la reescaleada, pero 
que la aproximamos por la real hoy. La justificacion para aproximar el radio 
es la siguiente 



Po(Mt) = Po(Mt) - ^(Mt) = Ro{Mt) U - f^ffjM ~ Po(Mt) (5.50) 

porque en la Ec. (15.49^ ya estamos a orden 0{—^). Asimismo 
111 1 

(5.51) 



\ MMt)) 

al orden en el que estamos trabajando. 

Por otra parte, justificamos que ko('"o) aproximadamente la distribu- 
cion real por lo siguiente. Las relaciones termodinamicas se deducen a partir 
del hamiltoniano sin gravedad. Al reescalear de manera de expresar las can- 
tidades termodinamicas en funcion de los que valen con las constantes fun- 
damentals actuales, encontramos las mismas relaciones termodinamicas que 
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encontrariamos hoy. Para ver que la funcion es la misma (porque Rq ^ Ri, 



Ko(ro) = KQ^r^erd) + ^ ) ■ {r^erd - t') (5.52) 



verd 



y i^yerd — f') ~ por lo cual lo despreciamos. 

Por lo tanto, el cambio en el radio del planeta en unidades de qbo sera 



Ai? = ^Mt) = —7^^, / Poir'ohoir'oK'dr',. (5.53) 
^ Rq Jo 



5.2. Modelado de los planetas 

La variacion temporal de las constantes fundamentales produce una al- 
teracion en el valor de la constante de Newton efectiva Gef- Esto, a su vez, 
provoca cambios en la estructura interna de los planetas. 

Ajustaremos el valor de AGef/G utilizando las cotas observacionales para 
la variacion de lo s radio s de la Tierra, Mercurio y la Luna, presentados por 



McElhinny et al.l (119781 ). Luego, haremos una estimacion para Aa/ao entre 
la epoca de formacion del sistema solar y la actualidad. 

Sin embargo, para poder relacionar AGej / G con AR es necesario modelar 
la estructura interna de los planetas mencionados anteriormente. 



5.2.1. Estructura de la Tierra 

Las ondas sismicas proveen una manera de explorar el interior de la Tie- 
rra, y su comportamiento predictivo hace posible obtener modelos de alta 
resolucion de las propiedades internas de nuestro planeta. 

Consideremos que la Tierra tiene simetn'a esferica, y que sus propiedades 
termodinamicas en funcion del ra dio est an bien descriptas por el Preliminary 
Reference Earth Model fPREMl dLav. T. and Wallace. T.lflQQsh (ver Tabla 



lA.ip . Los parametros de diclio modelo son la densidad p, las velocidades de 
ondas sismicas P y S", el coeficiente de atenuacion de shear Q^, el modulo de 
bulk adiabatico o incompresibilidad Ks, la rigidez /x, la presion y la gravedad. 

El modulo de bulk o incompresibilidad K se define como la resistencia de 
un material al cambio de volumen, cuando esta sujeto a presion. Si es el 
cambio de volumen por un cambio de presion AP, tenemos 
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es decir que K = donde k es la compresibilidad. 

Para estimar el valor de la integral involucrada en la expresion de Ai? 
para el caso de la Tierra, dividimos al planeta en capas y, llamando rj al 
radio de cada capa hacemos la siguiente aproximacion: 

rR N 

/ P{r')Ks{r')r''^dr' ~ V P{ri)Ks{ri)r1{ri - ri_i). (5.55) 
Usando esta aproximacion en la Ec. fl5.53p . obtenemos la relacion: 

= -288,09 km. (5.56) 



AG, 



G 



Considerando que el radio actual de la Tierra es i?o = 6371 km (segiin el 
modelo PREM), se obtiene 

0,04522 ( ^\ . (5.57) 



Ro \ G 

5.2.2. Estructura de Mercurio 

Para calcular la cota sobre el radio de Mercurio, necesitamos conocer la 
presion y la compresibi lidad en funcion de l a coor denada radial. Usamos los 



modelos calcu lados por [Harder fc Schubertl (120011 ) . de los cuales, elegimos el 



presentado en ISpohn et al.l (120011 ) . Este es un modelo con simetn'a esferica, 
de dos capas: un micleo de hierro, de densidad pc = 8gcm"'^ y un manto 
de silicatos, con una densidad de pm = 3,35 gcm~^. El radio de Mercurio es 
R = 2439 km, y el radio del micleo es Rc = 1860 km. 

Calculos para la presion 

Para encontrar la presion en funcion del radio, suponemos que el planeta 
esta en equilibrio termodinamico, es decir que se satisface 



dr 
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entonces 



r dP = P(r) - P(R) = Pir) = ~ f ^Il^^yiP^yidr' (5.59) 
Jr Jr r'^ 

donde hemos considerado que P{R) = 0. Teniendo en cuenta que 



m(r) = An p{x)x dx (5.60) 
Jo 

y Uamando 

g(r) = ^ (5.61) 

escribimos 

rR 

P{r) = I g{r')p{r')dr'. (5.62) 

J r 

Buscamos la expresion para la presion dentro y fuera del niicleo. 

Dentro del niicleo 

Para < r < separamos la integral fl5.62p en dos partes 

Rc rR 



P{r) = / g{r')p{r')dr' + / g{r')p{r')dr' . (5.63) 

Jr J Rc 

En la primera integral, r < es decir p(r') = p^? por lo tanto, la resolvemos 
de la siguiente manera 

Rc j-Rc ^ f-r' 

g{r')p{r')dr' = AnGpc J j p{x)x^dxdr'. (5.64) 
Ahora < x < r' < i?c, y por ende p{x) = pc- De esta manera, 

[ \{r')p{r')dr' = A'kGpI ! \ I x^dxdr' = ATrGpl^Rl - r^). (5.65) 

Jr Jr ''^ Jo 6 

Resolvemos ahora la segunda integral, donde Rc < r' < R, j por lo tanto, 
p{r') = pm 

R I r' 

g{r')p{r')dr' = ATtGpm / I p{x)x^dxdr' . (5.66) 
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Rc ^ Jo 



Separamos la integral en la variable x de la siguiente manera 



Rc rr' E)3 _ td3\ 



p{x)x dx = Pc X dx + p,n X dx = Pc^ + pm 

'0 Jo JRc 

Reemplazando en la ecuacion ( 15.66p . 



(5.67) 



g{r )p(r )dr = AnG 



(5.68) 



Por lo tanto, la presion dentro del niicleo, en funcion del radio, es 



pW(r) = AnG 



(5.69) 



Fuera del niicleo 

En el manto, Rc < r < R, por lo tanto 

"R rR 



P{r) = I g{r')p{r')dr' = ATrGpm [ ^ [ p{x)x^dxdr'. (5.70) 

Jr Jr Jo 



Aqm, Rc < r < r' < R, por lo tanto separamos 

f p{x)x^dx = Pc [ x^dx + pmj x^dx = pc^+Pm— — "^""^ (5-71) 

Jo Jo J Rc "J "J 

En consecuencia, la presion en el manto, en funcion del radio es 



pM(r) =47rG 



Pm i^Pc Pm) 



R r 6 



(5.72) 



Calculos para la compresibilidad 

Para la incompresi bilidad, usamos los valore s del modulo de bulk de la 
Tabla I del trabajo de [Harder fc Schubertl (120011 ) y consideramos que 



dK 

K{P) =Ko + 



(5.73) 
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Para el manto, = 120GPa, y = 4^25. Para el niicleo, 

^(core) ^ i27GPa, y |f = 2,20. Es decir que 

j^im) _ ^20 GPa + 4,25P('") (GPa) (5.74) 
K^'\P) = 127GPa + 2,20P(")(GPa) (5.75) 

(5.76) 

Recordamos que la incompresibilidad K es la inversa de la compresibilidad 
Finalmente, buscamos la expresion para 



AR 


1 


AG, A 




Gef J 






1 







R 



Rc j-R 



(5.77) 

la cual se puede integrar semi-analiticamente usando el programa Mathema- 
tica. Para determinar en que unidades de presion quedara el resultado de la 
integral, analizamos las unidades de las cantidades que la componen 

G = (6,67428 ± 0,00067) x 10"^^ kg~^ s'^ 
[r] = km 

[p] = gcm-l (5.78) 

Teniendo en cuenta que Igcni"^ = lO^^kgkm"^, y que IPa = lkgm~^s~^, 
Uegamos a que las unidades de P(r) seran: 

[P(r)] = 10^2 kgkm"^ s'^ = lO^Pa = IGPa. (5.79) 

Por lo tanto, si medimos las distancias en km, el resultado de la Ec. fl5.77p 
estara tambien en km. 

AP 

= -48,91km. (5.80) 



AG, 



G 

Considerando que el radio actual de Mercurio es Rq = 2439 km, se obtiene 
AP „ fAGe 



Rq 



-0,02005 • (5-81) 



109 



5.2.3. Estructura de la Luna 



El modelo de la Luna que vamos a considerar esta descripto en lKuskov fc Kronrod 
(Il998l ). Se trata de un modelo con simetria esferica, que consta de zonas di- 



ferenciadas: corteza, manto y un posible nucleo. En este estudio despreciare- 
mos el nucleo, dado que aun no se ha podido confirmar su existencia, radio 
y composicion. Por ot ra parte, hay evidencia s de que si existiera, deberia ser 
pequeno (< 460 km) (jWieczorek et al.ll2006l ). Por lo tanto, consideraremos 
que las caracterfsticas del manto inferior se extienden hast a el centro de la 
luna. 

Las expresiones para la densidad y compresibilidad en funcion de la coor- 
denada radial se obtendran a partir de las caracterfsticas de la estructura 
interna. Por otra parte, la presion en el interior de la Luna se puede describir 
con la formula 

r 



Pn 



1 - 



R 



(5.82) 



donde la presion central es Pq 
y r es la coordenada radial. 



47,lkbail!l, el radio medio es R = 1738 km, 



Corteza 

La corteza esta compuesta de una mezcla de Al-Ca, con una densidad 
promedio de 3 g cm~^ y se extiende desde la superficie hasta una profundidad 
de 58 km. En la interfaz corteza- manto, la densidad esta eri el siguiente rango: 
3,24 < Pm-cr < 3,32 gcm~^. En lKuskov fc KronrodI (119981 ) se consideran dos 
modelos, segiin el valor de la densidad en la interfaz: 



Modelo I 
Modelo II 



Pm-cr = 3,24 gem 3 
r = 3,32 gcm~^. 



P 



Para hallar la compresibilidad en la corteza, suponemos que la densidad 
tiene un crecimiento lineal con la profundidad h, es decir 



s dp 

P[h) = Pm-cr + -TT 

an 



{h hjYi—cr) 



''Recordemos las unidades: 

Ikbar = lO^bar = lO^dy/cm^ = lO^g cm~^seg~^ 



(5.84) 



(5.83) 
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donde h es la variable que caracteriza la profundidad y /i^ 



58 km es la 



profundidad de la interfaz corteza-manto. La densidad promedio es 



P 



h 



m—cr Jo 



+ 



dp 
dh 



{h ^m— cr) 



dh 



2 dh J , 



(5.85) 



de aqui deducimos que 



dp 
dh 



hr. 



{,Pm—cr P) ) 



por lo que podemos escribir la expresion para la densidad como 



P{h) = pm-cr + 



hr. 



Por otra parte, 



1 dV 



1 dp 



(^Pm~cr p) ■ i^h hm—cr) 



1 dp dh 1 dp 1 



VdPJs pdPJs pdhdP pdh[^) 



(5.86) 



(5.87) 



(5.88) 



y derivando la expresion (15.821) obtenemos la compresibilidad en funcion de 
la profundidad 



Ks{h) 



dp 



R 



(5.89) 



Pih) dhj,_ 2Po(l-4) 
En funcion de la coordenada radial, la compresibilidad queda expresada por 

1 dp\ i?2 



(5.90) 



p(r) dh J , 2Por 
Reemplazando los valores numericos queda (para 1680 < r < 1738 [km]): 



265,38 



(3,24 + ^(1680 - r)) r kbar 

353,83 1 
(3,32 + ^(1680 -r))r kb^ 



para el Modelo I 



paraelModeloII. (5.91) 
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Manto 



El manto lunar consiste en varias regiones de velocidades si'smicas cons- 
tantes, separadas por fronteras discontmuas. En las Tablas 15.11 (para el Mo- 
delo I) y 15.21 (para el Modelo II), se muestra la estructura interna de la luna, 
dando las propiedades de las distintas capas que conforman el manto, in- 
cluyendo los valores de las velocidades medias compresional (Vp) y de shear 
[Vs). Las distintas capas del manto se indican con: MS I (manto superior 
I), MS II (manto superior II), M int. (manto intermedio), y M inf. (manto 
inferior). 

Usaremos la siguiente relacion entre las velocidades si'smicas Vp y Vs, 
el modulo de bulk adiabatico Ks y la densidad p para sacar valores de la 
incompresibilidad para las distintas capas del manto: 



Si p se mide en gem ^, mientras que Vp, y en kms ^, las unidades de Ks 
son 10 kbar. Recordar que la compresibilidad ks es la inversa del modulo de 
bulk Ks. 

Finalmente calculamos 



Separamos la integral en 5 integrales, con ks constante en cada una de ellas 



Ks 



v'p - \vl 



(5.92) 



P 




(5.93) 




(5.94) 
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Por lo tanto, el resultado para el Modelo 1 es: 



AG, 



-83,622 km 



(5.95) 



mientras que para el Modelo 11, el resultado es: 

Ai? 

= -78,498 km. 



AG, 



G 



(5.96) 



Consider ando que el radio actual de la Luna es Rq = 1738 km, se obtiene 
para el Modelo 1 



AR 

Rn 



-0,04811 



AG, 



G 



(5.97) 



y para el Modelo 11 


AR 

Rq 


= -0,04517 






(5.98) 


region 


Vs 


Vp 


P 




Prof, inicial 


Prof, final 




[kms^-^] 


[kms^-^] 


[g cm^^] 


[kbar^-*^] 


[km] 


[km] 


corteza 






eq. 05.871) 


eq. 05.891) 





58 


MS 1 


4.493 


7.674 


3.239 


0.00097 


58 


270 


MS 11 


4.464 


7.674 


3.260 


0.00095 


270 


400 


M int. 


4.305 


7.568 


3.367 


0.00091 


400 


800 


M inf. 


4.549 


8.223 


3.389 


0.00074 


800 


nucleo 



Tabla 5.1: Estructura para el Modelo 1. 

5.3. Ajuste y resultados 

Para realizar el ajuste, utilizaremos los limites p ara la variacion de los 
radios planetarios dados en iMcElhinny et al.l (119781 ). Para la Tierra, estas 
cotas fueron obtenidas a partir de datos paleomagneticos con un metodo 
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region 


Vs 
[kms^-^] 


Vp 

[kms^"^] 


P 

[g cm^'^] 


[kbar^-*^] 


rroi. micial 
[km] 


rroi. nnai 
[km] 


corteza 






eq.fl5.87|) 


eq. (15.891) 





58 


MS I 


4.511 


7.748 


3.319 


0.00092 


58 


270 


MS II 


4.461 


7.772 


3.319 


0.00089 


270 


400 


M int. 


4.294 


7.526 


3.375 


0.00092 


400 


800 


M inf. 


4.532 


8.202 


3.376 


0.00074 


800 


nucleo 



Tabla 5.2: Estructura para el Modelo II. 

que compara las posiciones de regiones que presentan evidencias del campo 
magnetico terrestre correspondientes al mismo tiempo geologico, y estima 
cual sen'a el radio de la Tierra en esa epoca teniendo en cuenta el posible 
movimiento de placas tectonicas. La cota para la Luna se obtiene estudiando 
su superficie y determinando la edad de sus "mares" y "tierras altas" . Asi se 
estimo que la superficie de la Luna adquirio su forma actual hace 3,9 x 
10^ alios, cuando tuvieron lugar los ultimos eventos de mayores cambios 
de estructura. Finalmente, en las escarpas de la superficie de Mercurio hay 
evidencias de una contraccion del manto que ocurrio hace mas de 4 x 10^ 
anos. En la Tabla 15.31 se detallan los limites obtenidos para la variacion del 
radio de estos planetas. 



Planeta 


Ra/ Rq 


edad (10^ anos) 


Tierra 


1,020 ±0,028 


0,4 


Mercurio 


1,0000 ±0,0004 


3,5 


Luna 


1,0000 ±0,0006 


3,9 



Tabla 5.3: Paleoradios Ra dados en funcion del radio actual de la Tierra, 
Mercurio, y la Luna. 

Usaremos estos datos para realizar un ajuste de AGef/G y buscar una 
cota a la variacion de la constante de estructura fina a, entre el tiempo de 
formacion de los planetas y hoy. Supondremos que a tuvo un valor constante 
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en el pasado, aunque no necesariamente igual a su valor actual. Recordemos 



AR = ^{Mt) = Ro - Ro ^^ = i-S^. (5.99) 

Ro Rq 

siendo Rn/Rn el cocient e entre el paleoradio y el radio actual del planeta. En 



McElhinny et al.l (119781 ) se explica que un cambio homologo en el planeta no 
seria detectado por las observaciones, dado que escalea todas las dimensiones 
lineales de igual manera. Es decir, las cotas para la variacion del radio 
estan dadas para diferencias que excedan un simple cambio de es- 
cala. Es por eso que las cotas dadas en la Tabla 15.31 se comparan 
con Rq/Rq de la expresion ( I5.99p . 

Realizamos el ajuste para AGef/G. En las secciones anteriores obtuvimos 
valores iVj para el cociente entre AR/Rq y AGef/G en el caso de la Tierra, 
Mercurio y la Luna. Estos valores se resumen en la Tabla [5^ De esta manera 



Planeta 


Tierra 


Mercurio 


Luna (modelo I) 


Luna (modelo II) 




-0,04522 


-0,02005 


-0,04811 


-0,04517 



Tabla 5.4: Valores del cociente entre AR/Rq y AG^f/G para la Tierra, Mer- 
curio, y la Luna. 

Minimizamos el definido por 



X 



i=l 



1 _ AT, . 



(5.101) 



El valor de la cota que se obtiene para AGef/G, ya sea que se utilice el 
modelo I o II para la Luna, es 



AG, 



ef 



G 



0,000 ±0,011. 
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(5.102) 



Buscamos la cota para la variacion de a. Para ello, consideramos que 



G-G, 
G~ 



ef 



= 1 - 




(5.103) 



G 



Obtenemos el siguiente ajuste 



Aq; a 

= 1 = 0,000 ±0,011. 



(5.104) 



ao ao 



El limite encontrado resulta relevante por tratarse de una cota local, e 
independiente de otros resultados a tiempos similares. Si bien esta cota es 
mas debil que la obtenida a partir del reactor nuclear natural que opero en 
Oklo hace 1,8 x 10^ anos, la cual tiene una precision de 10~^, o de los de- 
caimientos cuya precision es de 10~^, esta cota resulta interesante por ser 
obtenida a partir de consideraciones muy generales, mientras que las otras 
aqui mencionadas se obtuvieron haciendo fuertes suposiciones en los modelos 
utilizados. 
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Capitulo 6 

Resumen y Conclusiones 



Esta tesis constituye un analisis de la variacion de la masa del electron 
nie, y de la constante de estructura fina a, en distintas etapas de la evolucion 
del Universe. 

Una parte importante de este trabajo se focaliza en la epoca de forma- 
cion del hidrogeno neutro, que se conoce como recombinacion del Universe. 
Hemos estudiado el escenario estandar de la recombinacion, y modificado 
las ecuaciones para las fracciones de ionizacion, para hacer posible que las 
constantes en estudio tomen valores arbitrarios. Sin embargo, mantenemos 
la forma de las ecuaciones, como propone la hipotesis adiabatica cuando las 
variaciones son pequefias. Analizamos las dependencias de las cantidades fisi- 
cas con las constantes fundamentals, incluyendo las cantidades relevantes en 
la recombinacion del helio, y otros fenomenos que modifican la historia de 
recombinacion, tales como las transiciones desde los estados tripletes del ni- 
vel n = 2 al n = 1 en el hidrogeno, o el efecto del hidrogeno neutro sobre 
la recombinacion del helio. Destacamos que las dependencias en este nue- 
vo escenario, que se de tallan en el Capitulo [21 es una de las contribuciones 
originales de esta tesis ( Scoccola et al. 2008al ). Modificamos las versiones del 
codigo Recfast que resuelven las ecuaciones de ionizacion en el escenario 
estandar y en el escenario actualizado, para poder utilizarlas luego para hacer 
los ajustes con los datos. 

Utilizando datos del fondo cosmico de radiacion, informacion acerca de 
Hq obtenida a partir del HST Key Project, y el espectro de potencias de las 
galaxias del catalogo 2dFGRS, hemos encontrado limites observacionales a 
la variacion de a y rrie entre la recombinacion y la actualidad. Realizamos 
analisis estadisticos, donde ajustamos el valor de un conjunto de parametros 
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cosmologicos, junto con la variacion de las constantes. Analizamos la varia- 
cion de cada una de ellas por separado, y el caso en que ambas constantes 
varian simultaneamente. Realizamos distintos ajustes utilizando diferentes 
conjuntos de datos, y comparamos nuestros resultados con los existentes en 
la literatura. En esta tesis se presentan resultados usando un prior de Hq que 
fue calculado a partir de objetos astronomicos cercanos. Estos analisis son los 
linicos en la literatura donde se considera la posibilidad de que el resultado 
final para el valor de Hq presentado por el HST Key Project este afectado por 
la variacion de las constantes, debido a que algunos de los objetos utilizados 
para medir Hq se encuentran a alto redshift, donde las constantes pudieron 
tener un valor diferente del actual. 

Se presentan cotas actualizadas con los datos de fondo cosmico de radia- 
cion incluyendo los datos de WMAP5 y el espectro de potencias del 2dFGRS, 
para la variacion de a, la variacion de rrie y la variacion conjunta. Ademas, se 
estimo una relacion fenomenologica entre la variacion de las constantes, que 
fue encontrada a posteriori, en el analisis donde se variaron ambas constantes 
de manera indepe n dient e. Comparamos nuestros analisis con los realizados 



en llchikawa et al.l (120061 ) . donde se propone que las variaciones estan rela- 
cionadas a partir de una ley de potencias, motivada en el contexto de una 
teon'a de cuerdas heteroticas, y realizan el ajuste presuponiendo dicha rela- 
cion. Dado que los valores del exponente de la ley de potencias utilizada por 
estos autores son positivos y que del ajuste independiente que realizamos en 
esta tesis estimamos una correl acion negativa entre la s variaciones de a y rrie, 
concluimos que la propuesta de llchikawa et al.l (l2006n de beria ser descartada. 
Los resultados estan publicados en iLandau et al.l (120081 ). 

Se presenta la linica cota existente para la variacion de rrie entre el des- 
acople de la materia y la radiacion y la epoca actual, y esta fue obtenida a 
partir de los datos mas actualizados del fondo cosmico de radiacion junto con 
el espectro de potencias del catalogo 2dFGRS. La cota enco ntrada con los 



datos d e WMAP3 y el espectro del 2dFGRS fue publicada en IScoccola et al. 



(l2008bl ). antes de que salieran publicados los datos de WMAP5. 

Los datos actuales comprenden mediciones muy precisas de las aniso- 
tropias del fondo cosmico de radiacion. Sin embargo, no se ha podido medir 
con precision la polarizacion. En los proximos aiios, una gran cantidad de 
experimentos (algunos en marcha, otros previstos) permitiran obtener me- 
diciones precisas de la polarizacion del fo ndo cosmi co de radiacion. Algunos 



de ellos son CMBP ol (iFraisse et al.ll2008r). QUIET (ISamtleben 



(iNorth et al.l 120081 ). BICEP (iTakahashi et al.ll2008l ). y Planck (lDupacll2007l ) 



2008), Clover 
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Pero a pesar de que las misiones futuras tendran mayor precision, aiin asi no 
permitiran ajustar la variacion de las constantes fundamentals con una pre- 
cision mejor que 10 

La cotas obtenidas observacionalmente para la variacion de las constan- 
tes fundamentales constituyen una de las pocas herramientas con las cuales 
testear teorias que intentan unificar las cuatro interacciones de la naturale- 
za, como hemos mencionado en la Introduccion. En esta tesis, analizamos 
el modelo fenomenologico de B arrow- Magueijo, que propone una variacion 
de la masa del electron inducida por cambios de un campo escalar en el 
espacio-tiempo. Mejoramos las soluciones originalmente presentadas por di- 
chos autores, teniendo en cuenta de manera mas detallada la evolucion del 
factor de expansion del Universo, y haciendo menos aproximaciones que en 
el trabajo original. Estimamos los parametros del modelo usando las cotas 
para la variacion de la masa del electron que se obtienen en esta tesis a par- 
tir de los datos del fondo cosmico de radiacion, y otros datos astronomicos 
y de laboratorio correspondientes a distintos tiempos cosmologicos. Por otra 
parte, utilizando cotas observacionales obtenidas en experimentos de Eotvos 
que testean la validez del Principio de Equivalencia Debil, encontramos un 
Ifmite inferior para el parametro que aparece en el Lagrangiano propuesto en 
este modelo, que viene de la variacion en la masa del electron inducida por 
cambios espaciales del campo escalar. Comparando esta ultima con la cota 
obtenida con los datos a distintos tiempos, cuando estudiamos evolucion tem- 
poral del campo, se llega a la conclusion de que el modelo e s inconsistente 



ydebe ser descartado. Estos resultados fueron publicados en IScoccola et al. 
fl2008bh . 

Por ultimo, estudiamos la variacion de las constantes fundamentales entre 
la actualidad y la epoca de formacion de los planetas del Sistema Solar. A 
partir de consideracion muy generales, analizamos como se modificaria el 
radio de un planeta por variacion de las constantes a y vtif.. Explicamos por 
que la variacion de la masa del electron no induce cambios en el radio que 
sean factibles de ser observados. Obtenemos una cota para la variacion de la 
constante de estructura fina a partir de Ifmites sobre la variacion del radio 
de la Tierra, la Luna y Mercurio, desde la epoca de los liltimos cambios 
estructurales mayores en la superficie de los mismos. Si bien es una cota 
debil, es importante notar que se trata de un resultado local, e independiente 
de otras medidas. Ademas, fue obtenido a partir de consideraciones muy 
generales. 
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Apendice A 
Modelo PREM 



Radio 


Prof. 


Dcnsidad 


Vp 


Vs 


Q, 


^ adiah 




Presion 


Gravedad 


(km) 


(km) 


{g/cm?) 


(km/s) 


(km/s) 


(kbar) 


(kbar) 


(kbar) 


(cm/s^) 





6371 


13.08 


11.26 


3.66 


85 


14253 


1761 


3638.5 





200 


6171 


13.07 


11.25 


3.66 


85 


14231 


1755 


3628.9 


73.1 


400 


5971 


13.05 


11.23 


3.65 


85 


14164 


1739 


3600.3 


146 


600 


5771 


13.01 


11.2 


3.62 


85 


14053 


1713 


3552.7 


218.6 


800 


5571 


12.94 


11.16 


3.59 


85 


13898 


1676 


3486.6 


290.6 


1000 


5371 


12.87 


11.1 


3.55 


85 


13701 


1630 


3402.3 


362 


1200 


5171 


12.77 


11.03 


3.51 


85 


13462 


1574 


3300.4 


432.5 


1221.5 


5149.5 


12.76 


11.02 


3.5 


85 


13434 


1567 


3288.5 


440 


1221.5 


5149.5 


12.16 


10.35 








13047 





3288.5 


440 


1400 


4971 


12.06 


10.24 








12679 





3187.4 


494.1 


1600 


4771 


11.94 


10.12 








12242 





3061.4 


555.4 


1800 


4571 


11.8 


9.98 








11775 





2922.2 


616.6 


2000 


4371 


11.65 


9.83 








11273 





2770.4 


677.1 


2200 


4171 


11.48 


9.66 








10735 





2606.8 


736.4 


2400 


3971 


11.29 


9.48 








10158 





2432.4 


794.2 


2600 


3771 


11.08 


9.27 








9542 





2248.4 


850.2 


2800 


3571 


10.85 


9.05 








8889 





2055.9 


904.1 


3000 


3371 


10.6 


8.79 








8202 





1856.4 


955.7 


3200 


3171 


10.32 


8.51 








7484 





1651.2 


1004.6 



Tabla A.l: Parametros del modelo PREM 
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Radio 


Prof. 


Densidad 


Vp 


Vs 


Q, 


K-adiab 




Presion 


Gravedad 


(km) 


(km) 


(g/cm^) 


(km/s) 


(km/s) 


(kbar) 


(kbar) 


(kbar) 


(cm/s^) 


3400 


2971 


10.02 


8.19 








6743 





1441.9 


1050.6 


3480 


2891 


9.9 


8.06 








6441 





1357.5 


1068.2 


3480 


2891 


5.56 


13.71 


7.26 


312 


6556 


2938 


1357.5 


1068.2 


3600 


2771 


5.5 


13.68 


7.26 


312 


6440 


2907 


1287 


1052 


3800 


2571 


5.4 


13.47 


7.18 


312 


6095 


2794 


1173.4 


1030.9 


4000 


2371 


5.3 


13.24 


7.09 


312 


5744 


2675 


1063.8 


1015.8 


4200 


2171 


5.2 


13.01 


7.01 


312 


5409 


2559 


957.6 


1005.3 


4400 


1971 


5.1 


12.78 


6.91 


312 


5085 


2445 


854.3 


998.5 


4600 


1771 


5 


12.54 


6.82 


312 


4766 


2331 


753.5 


994.7 


4800 


1571 


4.89 


12.29 


6.72 


312 


4448 


2215 


655.2 


993.1 


5000 


1371 


4.78 


12.02 


6.61 


312 


4128 


2098 


558.9 


993.2 


5200 


1171 


4.67 


11.73 


6.5 


312 


3803 


1979 


464.8 


994.6 


5400 


971 


4.56 


11.41 


6.37 


312 


3471 


1856 


372.8 


996.9 


5600 


771 


4.44 


11.06 


6.24 


312 


3133 


1730 


282.9 


998.8 


5650 


721 


4.41 


10.91 


6.09 


312 


3067 


1639 


260.7 


1000.6 


5701 


670 


4.38 


10.75 


5.94 


312 


2999 


1548 


238.3 


1001.4 


5701 


670 


3.99 


10.26 


5.57 


143 


2556 


1239 


238.3 


1001.4 


5771 


600 


3.97 


10.15 


5.51 


143 


2489 


1210 


210.4 


1000.3 


5871 


500 


3.84 


9.64 


5.22 


143 


2181 


1051 


171.3 


998.8 


5921 


450 


3.78 


9.38 


5.07 


143 


2037 


977 


152.2 


997.9 


5971 


400 


3.72 


9.13 


4.93 


143 


1899 


906 


133.5 


996.8 


5971 


400 


3.54 


8.9 


4.76 


143 


1735 


8.6 


133.5 


996.8 


6061 


310 


3.48 


8.73 


4.7 


143 


1630 


773 


102 


993.6 



Tabla A. 2: Parametros del modelo PREM (cent.) 
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Radio 


Prof. 


Densidad 


Vp 


Vs 


Q, 


Kadiab 




Presion 


Gravedad 


(km) 


(km) 


(g/cm^) 


(km/s) 


(km/s) 




(kbar) 


(kbar) 


(kbar) 


(cm/s^) 


6106 


265 


3.46 


8.64 


4.67 


143 


1579 


757 


86.4 


992 


6151 


220 


3.43 


8.55 


4.64 


143 


1529 


741 


71.1 


990.4 


6151 


220 


3.35 


7.98 


4.41 


80 


1270 


656 


71.1 


990.4 


6186 


185 


3.36 


8.01 


4.43 


80 


1278 


660 


59.4 


989.1 


6221 


150 


3.36 


8.03 


4.44 


80 


1287 


665 


47.8 


987.8 


6256 


115 


3.37 


8.05 


4.45 


80 


1295 


669 


36.1 


986.6 


6291 


80 


3.37 


8.07 


4.46 


80 


1303 


674 


24.5 


985.5 


6291 


80 


3.37 


8.07 


4.46 


600 


1303 


674 


24.5 


985.5 


6311 


60 


3.37 


8.08 


4.47 


600 


1307 


677 


17.8 


984.9 


6331 


40 


3.37 


8.1 


4.48 


600 


1311 


680 


11.2 


984.3 


6346.6 


24.4 


3.38 


8.11 


4.49 


600 


1315 


682 


6 


983.9 


6346.6 


24.4 


2.9 


6.8 


3.9 


600 


753 


441 


6 


983.9 


6356 


15 


2.9 


6.8 


3.9 


600 


753 


441 


3.3 


983.3 


6356 


15 


2.6 


5.8 


3.2 


600 


520 


266 


3.3 


983.3 


6368 


3 


2.6 


5.8 


3.2 


600 


520 


266 


0.3 


982.2 


6368 


3 


1.02 


1.45 








21 





0.2 


982.2 


6371 





1.02 


1.45 








21 








981.5 



Tabla A.3: Parametros del modelo PREM (cent.) 
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Apendice B 
Integrales I 



Tenemos que resolver est a ecuacion diferencial: 

es decir, algo de la forma: 

a + a; , f dx f x dx 

-dx = a — — 27 r + / ~. — 27 7 (-^ 



sinh (c + bx) J sinh (c + bx) J sinh (c + bx) 



I II 



Sea 

y = c + bx ^ dy = b dx ^ dx = - dy 

b 



T f dx a f dy ci <^ / , \ 

1 — a 27 7 — T 2~ ~ ~T cotny = — - cotn [c + bx) 

J sinh (c + bx) b J sinh y b b 

(B 

ya que 

d(— cothy) d / cosh|/\ sinh^ y — cosh^ y 1 
dy dy \ sinh y J sinh^ y sinh^ y 

Por otra parte, 
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I si 



X dx 



sinh^ (c + bx) IP' J sinh^ y IP' 



1 [ {y- c)dy 1 



y dy c 

+ 7^ cotn y 



y dy 
sinh^ 1/ 



— c 



Calculamos ahora 



si 



1/ c?l/ 



sinh y 



sinh^ y 6^ 



-|/ coth y + coth 



donde hemos integrado por partes. Ademas: 

coshy dy 



cothy dy 



sinhy 



In (sinh y) 



Entonces, 



y dy , 1 / . 1 \ 

2~ — ~y coth y + In (smJi y) 

sinh y 



Por lo tanto: 



dy 



sinh y 

coth 3/ 



:b.5) 



:b.6) 



(B.7) 



(B.8) 



1 c 
II = — (-y coth y + In (sinh y)) + —coth y 

I ^ _|_ ) 1 c 

= — coth (c + fcx) + — In (sinh (c + bx)) + — coth (c + bx) 

JjZ Jj2 JjZ 



En definitiva: 

+ x 



!— 

J sinh^( 



(B.9) 



C^; -|- a) 1 

dx = --^^ — — - coth (c + fex) + — In (sinh (c + (B.IO) 



(c + 6x) b 62 

Reemplazando las cantidades originates, Uegamos a esta expresion: 
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X 



coth (c' + ^ v^^o) + ^ ('^^^ + ^a/^^o)) 



(B.ll^ 



Reemplazando el valor de tq, reordenando terminos y llamando 7 = 
a/ a/q , la solucion para la epoca de dominio de materia mas constante 
cosmologica es: 



M 



Va^r coth (^C + ^ v^^^ 



+ 



-coth(C" + -v^r 



C + In f 7 + Vl + 72 



7 



+ 



Esta solucion vale desde el momento de ajuste, (donde los dos regfme- 
nes se superponen), hasta hoy, Tq. Evaluamos la Ec. (IB.lip en ta- 



^0 + 



M n 



Ao 



2 



Hp A 
M 



3^ 



Ao 



3\/^Ao 



+ 777^ — In 7 



-l/2x 



M 



7 



91] 



Ao 



(B.12) 
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Apendice C 
Integrales II 



La ecuacion a integral es: 



't'A 



— drj 



VA 



HoV'- 



,2 / "Mq 



Por lo tanto, integramos la siguiente integral indefinida: 



M 
2 Ho 



7] 



V + V^) + A 



Hon 

M 



2 I 



+ 



2 Hi ^ 
A 



6 ./n 



«0 



4 



RO 



=77 + 1 



dr] 



dr) 



Han 



4 ^/n 



dr] 



Sea (3 = Qmo/V^o, a ^ (3/6 y b ^ /3/4 {a = 2 6/3), 



7 = 



M 



at] + 1 A f df] 



Primer o resolvemos: 



(C.l) 



(C.2) 



(C.3) 
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ar] + 1 



drj — — 



br] + b/a 



a 
b 
a 



{bri + 1)' 
dfj f b 



{bri + 1) 
In [hr] + 1) + 



+ 



dr] 



- 1 



dfj 



a — b 



62 {brj + 1) 



-In (677 + 1) 



4 



36 (677 + 1) 



21n( -/?77+l 



1 



{jPv + 1) 



(C.4) 



La otra integral que se debe resolver es: 



/ 



7?2 (677 + 1)^ J V J J {hq + l) J (hq + lf 

= -261n77 - - + 261n(677 + l) -6- (C.5) 

'77 ' (677 + 1) 



Por lo tanto, 



M 



+ 



6 6//2v^ 
A 6 



2 In (677 + 1) 



1 



—2 In 77 h 2 In (677 

brj 



(677 + 1) 

1)- 



+ 



(677 + 1). 



(C.6) 
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Reemplanzando las cantidades originales, la solucion es: 



M 



+ 



3 HqQmo V 



(21n(6r7 + 1) 



1 



[br] + 1) 



+ 



2 M 

3 Hq^Imo 
A 



-2^^^-^ + 21n(6, + l)-^^ 



21n(6r/^ + l) 
1 



+ 



-2hir]A - 



ipr]A + l), 
h21n(6r/^ + 1) 



(C.7) 



con 



1 



1 VI 



Mo 



4^ 



Ro 



Considerando que es continuo, y reemplazando 0^ (de la Ec. ( IB. 121) ) 
en la ec. ( ]C.7p . el campo para la epoca dominada por materia y radiacion 
es: 

2M 1 / 1 



kR 



+ - 



3 Hq VImo V 
1 A VImo 



2 \n{br] + 1) 



(br] + 1) 



+ - 



2 M 1 



1 



+ 
1 



+ 



21n(6r/A + l) 



(^ + ?/(OV?TT-?i„,r"=)- 



3 i^o 

(sinh-i(7) - sinh-i(r^/') + /(O) ^^^^^ + ^ ln(7) 



A Q 



Mo 



4 '""--5^ 



- ln{br]A + 1) 



1 2v^TT 

4(6r/^ + 1) ^ 3 



donde 7 se definio en el Apendice [B] y 



2 + (e-2)v/r+i 



(C.9) 
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Ademas: 

sinh"^ x = ]n(x + Vx^ + l) (C.IO) 
A es la constante de integracion, y 

M = (C.ll) 
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